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PREFACE 


This  publication1  presents  results  of  a  seminar  organized  by  and 
held  at  DCEC  on  26  November  1979  for  the  purpose  of  exchanging  infor¬ 
mation  among  Industry  and  government  representatives  on  the  subject 
of  frequency  selective  fading  and  its  effects  on  digital  LOS  radio 
performance.  The  seminar  was  prompted  by  DCEC  interest  in  frequency 
selective  fading  effects  on  LOS  radio  and  by  similar  Interest  shown 
by  other  organizations  as  indicated  in  the  available  literature.  It 
is  hoped  that  results  of  this  seminar  will  lead  to  further  exchange  of 
analytical  and  test  results  in  the  future. 

Each  seminar  speaker  provided  a  copy  of  the  materials  used  in  his 
presentation,  and  these  copies  are  presented  herein.2  In  some  cases, 
these  are  only  the  visuals  that  were  used  with  extemporaneous  speeches. 

In  other  cases,  copies  of  supplemented  materials  prepared  by  authors 
than  the  speakers  are  included  because  these  materials  were  a  key  element 
in  the  presentation. 

Further  questions  regarding  the  proceedings  of  the  seminar  may  be 
directed  to  Dr.  David  R.  Smith,  (COM}  703-437-2316,  (VON)  364-2316. 


1  This  publication  Is  the  second  in  a  series  of  DCEC  reports  on  frequency 
selective  fading  and  its  effects  on  DCS  transmission  system  performance. 

2 

For  those  presentations  which  used  previously  published  articles,  per¬ 
mission  of  the  author(s)  has  been  obtained  and  copyright  information 
has  been  provided  in  accordance  with  requirements  of  the  IEEE  and 
Bell  Telephone  Laboratories. 
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1.  ASSESSMENT  OF  SELECTIVE  FADING  EFFECTS  ON  DIGITAL  LOS  LINKS 


David  R.  Smith 

Defense  Communications  Engineering  Center 
Res ton,  VA 


NOTE:  The  materials  provided  in  this  section  are  copies  of 
the  visuals  used  with  the  extemporaneous  remarks  given 
by  Dr.  Smith. 
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ASSESSMENT  OF  FREQUENCY  SELECTIVE  FADING  ON  DIGITAL  MICROWAVE  LOS  LINKS 


BACKGROUND 


POWER  RATIO  OF  DIRECT  TO  REFRACTED  COMPONENT  (Pj/P^ 


Refractive 
Layers  or 
Irregularities 


PROPAGATION  STUDIES/MEASUREMENTS 


IV.  TEST  AND  EVALUATION  EFFORTS  (CONTINUED 


BER  AND  EACH  DIVERSITY  CHANNEL 


WSCSBING  PAOS  BUNK-NQff  Jli>£D 


2.  MULTIPATH  DEGRADATION  OF  DIGITAL  MICROWAVE  LINK  AVAILABILITY 


John  K.  Webb 

The  MITRE  Corporation 

Bedford,  MA 


NOTE:  The  material  provided  in  this  section  is  a  text  of  the 
remarks  given  by  Mr.  Webb. 
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The  MITRE  Corporation  provide.-;  technical  support  to  the  Air  Force  AFT  El. 
Program  Office  which  has  Program  Management  responsibility  for  the  DEO  system. 

DEB  will  introduce  all-digital  transmission  to  the  European  DC3.  Stage  I 

£ 

now  largely  completed,  uses  3-level  partial  response  12.$  Mb/s  with  FM  radios. 
Stages  2  thru  4  will  install  117  links  of  up  to  26.3  HB/s  v.'ith  the  Army 
DRAMA  radios  using  QPRS  at  4  and  8  GHz  for  2  bits  per  Hz  bandwidth  efficiency. 
The  DRAMA  radio  has  an  alternate  1.0  bit/Hz  QP3K  capability. 

In  July  78,  I  was  asked  to  examine  the  performance  that  should  be  expected 
on  the  longer  DEB  links  with  higher  rate  QPRS  modulation.  Concern  was  suggested 
by  several  BTL  and  BNR  papers  presented  at  ICC-78.  They  showed  that  outages 
due  to  multipath  degradation  of  signal  quality  could  exceed  outages  due  to  signal 
power  fading  below  the  receiver  threshold.  Barnett's  work  measured  considerable 
outage  due  to  multipath  degradation  alone  on  a  26  mile  link  with  a  margin  so 
high  that  fading  below  threshold  would  not  happen.  Papers  by  Prabhu,  Greenstein, 
Rummler  and  Emshwiller  showed  that  the  outages  observed  by  Barnett  were  consis¬ 
tent  with  the  intersymbol  interference  to  be  expected  with  the  time  dispersive 
multipath  delays  observed  during  the  test  period. 

({2])  Jakes'  paper  offered  a  statistical  model  for  cumulative  outage  prediction 
that  was  based  in  several  sources.  It  embodies: 

.  Ruthroff's  atmospheric  delay  model  which  assumes  atmospheric 
conditions  that  allow  delay  to  increase  as  the  third  power 
of  distance. 

.  Coefficients  for  mean  values  of  delay  from  Prabhu  and  Greensteir,. 

.  and  coefficients  that  fit  the  model  expression  to  samples  of 
outage  data  measured  on  several  test  links  and  measured  multipath 
tolerance  for  several  modulation  methods. 
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c  .  .  none  o  f  which 

: I . h . j v •  v.  data.  i.'c-ver- 
iy  o?r-l  v.jj  to  the  lowsr 


LNR  provided  a  basic 
thir  check  cn  the  Jakes 
t?  of  the  relative  multi- 


The  Jakes  model  fits  outage  measurements  f .  . 
v/ere  accompanied  by  reported  or  measured  collateral  ; 
theless,  it  is  the  only  model  we  have  that  can  he  dir,' 
data  rates  and  the  longer  distances  of  the  DEB  system. 

The  work  reported  by  Barber,  .  nderson  and  Fa  tel 
for  including  some  form  of  diversity  improvement,  and  c 
model  at  a  33  mi 1 e  distance.  Also  it  provides  an  estr. 
path  tolerance  of  QPRS  and  4  or  8  PSK.  BNR  found  GPRS  to  he  nearly  the  same  as 
8  PSK;  later  private  communication  with  Jakes  suggests  bilk's  conclusion  to  be 
based  on  linear  amplitude  and  slope  effects  and  that  QPKS  should  perform  more 
like  QPSK  with  multipath  as  it  exists  in  nature. 

((3^  This  graph  shows  the  Jakes  model  calculation  of  r,cn-divc-’~sity  cumulative 
outage  for  26.8  Mb/s  DEB  links  at  4  and  8  GHz  as  a  function  of  distance.  The 
curve  shows  outage  time  growing  exponentially  with  distance.  The  lower  curve 
shows  the  effect  of  diversity  improvement  in  the  amount  measured  by  BNP  in 
their  tests.  However,  the  DRAMA  radios  use  hitless  baseband  switched  diversity, 
not  coherent  IF  combining  and  adaptive  equalization  as  in  the  CNR  radios. 
Application  of  BNR's  diversity  improvement  to  a  system  using  a  different  diversity 
technique  is  of  course  not  supportable;  rather  wo  need  to  refine  our  understand¬ 
ing  of  how  a  switching  diversity  system  would  work  in  terms  of  correlation  of 
multipath  degradation  in  the  diversity  pair.  The  inadequacy  of  the  cumulative 
prediction  is  further  highlighted  by  the  DCS-allocated  LOS  performance  require¬ 
ment;  it  is  expressed  in  terms  of  probabilities  of  outage  of  several  durations, 
not  as  cumulative  outage.  The  results  did  show  that  multi  path  degradation  may 
significantly  add  to  the  outages  expected  from  signal  p"..er  faring  alone.  The 
concern  was  sufficient  to  develop  an  extended  mole,  ti.rf  ieu  a  proposal  tor 
tests  to  refine  uncertainties  that  wore  inherent  in  tl,-  . .ore  aicut  the 

tests  later,  but  first  the  ox  tended  model. 
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The  model  was  extended ,  largely  by  layered  hypotheses,  to  provide  a  nuasi-ira' 
basis  for  accommodating  diversity  features.  The  Jakes  model  is  used  initially 
to  calculate  cumulative  non-diversity  multipath  outages;  the?  cumulative  out.  eg 
is  allocated  to  each  fade  depth  accord! rg  to  this  population  reported  by 
Barnett.  It  is  a  histogram  showing  the  percentage  of  the  total  multi  path  outage 
associated  with  each  fade  depth.  Here  we  assume  that  fades  of  say  20  dB  will 
contribute  nine  percent  of  the  cumulative  outage  for  either  Bi'L  or  DOD.  Now  each 
fade  depth  has  an  average  duration  and  a  distribution •(Ts^  Here  we  used  relation 
ships  suggested  by  Vigants  for  signal  fading  with  the  implicit  assumption  that  in 
fades  were  multipath  outages  occur,  the  duration  of  the  outage  will  correspond 
with  duration  of  the  fade.  Thus  we  can  sum  the  outage  time  contribution  of  each 
fade  depth  above  threshold  for  a  specified  fade  duration  (t).  For  simultaneous 
fades  of  both  channels  below  threshold,  we  used  conventional  expressions  from 
Vigants--there  we  don't  care  about  multipath  degradation  because  both  channels 
have  inadequate  signal  level.  For  conditions  where  the  protected  channel  fades 
below  threshold,  we  put  in  a  provision  to  add  the  outage  contribution  due  to 
multipath  degradation  of  the  protection  channel.  Lack  of  data  so  far  prevents 
us  from  exercising  this  feature  of  the  model.  Outage  reductions  due  to  switch¬ 
ing  diversity  protection  are  included  by  multiplying  the  outage  contribution  of 
each  fade  depth  by  the  correlation  coefficient  of  multipath  outage  on  a  diversity 
pair  for  each  fade  depth--that  data  of  course  is  not  available  from  measurements 
reported  in  the  literature. 

(fsj)  The  non-diversity  outages  summed  by  this  extended  model  are  shewn  in  this 
table  for  a  5C  mile  link  with  a  fade  margin  of  32  dB.  Mere  we  see  that  outages 
of  all  durations  are  predominantly  due  to  multipath.  The  38, £00  figure  is  from 
Jake's  model  but  allocated  for  fade  depths  less  chan  3?  o.Z  according  to  Barnetts 


2-4 


histogram.  Fo  the  right  we  ere  outage  Lima  decreasing  with  outage  duration 
as  Vigants  average  duration  statistics  are  apoliod  to  each  fade  depth. 
Below,  a  line  of  small  numhe. i  shu./s  that  oucages  due  to  the  signal  fading 
below  threshold  tend  toward  insignificance  beside  multipath  caused  outages. 
Also,  we  see  that  outage  ti:;_  decreases  slowly  as  outage  duration  increases. 
This  is  because  fade  durations  are  long  at  the  shallower  fade  depths  at 
which  much  of  the  multipath  degradation  occurs. 

pit  of  the  same  numbers  computed  over  a  range 
of  distances.  The  same  32  dB  margin  applies  in  all  cases.  Again  we  see  the 
same  alarming  growth  of  outages  with  distance  for  all  outage  durations,  pri¬ 
marily  because  the  Jakes  and  Kuthroff  models  allow  outages  to  grow  as  the  third 
power  of  distance.  Observations  of  multipath  delay  on  longer  links  do  not 
support  that  such  growth  really  happens;  therefore  we  should  expect  these 

curves  to  truncate  in  some  fashion  with  a  trend  in  evidence  somewhere  between 

1 1 

30  and  40  miles.  But  if  we  erect  an  ordinate  at  say  40  miles  we  see  over  10' 
seconds  of  .2  second  outages  that  have  to  be  dealt  with  by  the  diversity 
switch. 


& 


The  next  figure  is  a  gv 


So  far,  vie  have  seen  hew  one  model  works  without  outage  reductions  or 
improvements  available  with  the  hitless  diversity  switch.  For  a  diversity 
evaluation,  fading  correlation  coefficients  were  taken  from  a  sample  of  data 
for  30-foot  vertical  an  Ur, no  separation,  the  coefficients  approached  unity  in 
the  shallower  fade  region  we  are  concerned  about.  As  expected,  that  run  showed 
negligible  outage  reduction.  We  could  see  in  the  data  why  switching  diversity 
works  in  analog  systems  —  they  fa!1  below  threshold  in  deep  fades  and  deep  fades 
become  quite  uncorrelatn.j  ir.  Urne.  Thus  a  low  probability  of  simultaneous 
fading  below  threshold  a:--*  r;  U-i-rr  of  outage.  Mow  wo  can  conclude 


it  conditions 


Li  pa  to  dear; 


in  shallow  fade: 


a  di.crsiby  pair 


si  tv  improvement  in  the 
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digital  radio  that  does  not  otherwise  ameliorate  multipath  by  adapting  the 
receiver  tc  the  corrupted  waveform. 

This  modeling  activity  then  suggests  several  features  for  a  link  test 
to  refine  our  ability  to  predict  the  performance  of  longer  digital  radio 
links.  Two  of  the  more  important  features  are  distance  dependence  and 
correlation  of  diversity  degradation  as  a  function  of  antenna  installation 
criteria.  (f8j)0n  distance  dependence,  there  are  several  reasons  to  suggest 
that  the  curves  of  outage  versus  distance  increase  more  slowly  with  distance 
than  the  model  shows.  The  10  ns  of  delay  predicted  by  Ruthoff's  model  and 
used  ir.  the  Jakes  model  is  in  agreement  with  measurements  on  the  26  mile 
Georgia  path;  however,  that  value  has  to  increase  as  the  cube  of  distance 
to  support  use  of  the  model  at  greater  distances.  The  atmospheric  delay 
observed  by  channel  probes  is  not  known  to  have  ever  exceeded  30  ns  on  any 
path  of  up  to  65  miles.  When  the  BNR  system  parameters  are  plugged  into  the 
Jakes  model,  the  outage  calculated  at  32  miles  is  3  times  greater  than  was 
measured  in  BNR's  non-diversity  test.  These  facts  suggest  that  realistic 
delay  values  for  distances  longer  than  30  miles  will  increase  more  slowly 
than  Ruthroff's  model .  This  seems  reasonable  because  longer  paths  involve 
higher  terminal  elevations;  accordingly  a  smaller  fraction  of  a  long  path 
involves  the  lower  elevation  atmosphere  where  many  anomolous  refractive 
conditions  arc  known  to  occur.  If  the  multipath  degradation  is  due  to  lower 
atmosphere  anomolies  than  performance  should  not  degrade  as  severely  with 
path  length.  Tests  on  long  paths  are  required  to  resolve  this  uncertainty. 

Tne  diversity  question  is  one  of  how  the  correlation  of  multipath 
degradation  on  a  diversity  pair  is  rented  to  antenna  installation,  if 
correlation  of  fade  intervals  s  closely  tied  ;o  intervals  cf  mul  tipatli 
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lie; :  irii,i',  then  30  feet  or  less  cf  spacing  appears  Inadeouate;  however,  ITS 

o(':erva':  ior:s  a.  Ft.  iiucju  of  a  100  foot  spaced  antenna  pair  showed  very  iKiCOi - 
rente.:  multipach  outages.  Certainly,  measurements  of  multipath  degradation 
need  to  to  made  with  antenna  spacing  as  a  variable  The  tests  should  establish 
ante:  na  installation  criteria  that  enable  realization  of  diversity  improvement 


wi ti •  the  switching  technique. 

This  figure  shows  on  a  data-rate/di stance  coordinate  system  where  the  test 
data  is  in  relation  to  DEB.  DEB  is  shown  as  a  variable  distance  with  numbers 
besiuc  the  ordinate  to  show  the  population  of  links  in  each  decade  cf  distance. 
There  is  not  much  near  DEB.  The  BTL  and  BNR  points  are  in  the  lower  right  corner. 
The  ETL  point  might  move  in  the  future  toward  higher  data  rates  but  surely  not 
to  longer  distances.  The  Haw-tel  data  at  90  miles  is  an  exception--but  it  didn't 
work.  Haw-tel  observed  whole  days  of  outage  which  correlated  with  an  elevated 
inversion  layer.  The  Army  tested  a  link  similar  to  DEB,  12.9  Mb/s,  over  33  miles 
at  Ft.  Huachuca.  Fading  was  not  observed  which  is  typical  of  microwave  tests 
in  the  high  dry  climate.  Therefore,  none  of  the  other  existing  work  provides 
data  .on  a  system  like  DEB  in  the  European  location.  We  have  recommended  testing 
the  DEB  equipment  on  a  long  path  in  a  climate  approximating  central  Europe.  Two 
paths  have  been  recommended  in  the  Massachusetts/New  Hampshire  region,  one  45 
miles,  the  other  68  miles.  A  test  thru  a  summer  fading  season  cn  either  of  these 
path.;;  should  greatly  refine  our  ability  to  predict  DEB  digital  radio  performance 


in  Europe. 

©  AS  part  of  our  recommendations  to  do  live  link  testing  considerable  thought 
has  been  given  to  location  dependence  in  terms  of  atmospheric  phenomenology. 

The  first  question  we  addressed  is  "Why  not  do  the  testing  it  Ft.  Huachuca  where 
y  :.os  an  established  test  facility?  First  the  ava  ';<■>  atmospheric 
rt  statistics  for  Germany  and  a  high  dry  cl  imatic  ai  ca  were  compared— 

*•  i".  -  i :  in  oir-ck,  El  Paso,  Texas  is  red  and  El  Paso  is  climatically 
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similar  to  Ft.  Huachuca.  Summer  ir.  solid  lines  and  winter  is  dashed  lints. 

The  red  lines  (El  Paso)  show  little  winter  to  summer  difference  arid  are  r.ct 
very  different  from  the  black  dashed  Vine-winter  in  Germany.  The  refractive 
index  extremes,  both  sub  and  super  refraction  shown  by  the  solid  black  line 
which  is  the  summer  fading  season  in  Germany,  is  not  in  evidence  in  the 
El  Paso  or  Ft.  Huachuca  climate.  It  is  r,o  suprise  that  microwave  links  work 
extremely  well  at  Ft.  Huachuca  and  if  the  deployment  area  was  Algeria  or 
Afganistan,  Ft.  Huachuca  would  be  a  representative  test  area.  It  is  for  these 
reasons  that  we  have  recommended  New  England  as  a  test  location;  it  is  more 
representative  of  the  central  Europe  climate  and  terrain  than  are  Arizona, 
Washington,  D.C.,  or  Eglin  AFB  Florida. 

(0JThe  proposed  test  would  operate  a  space  diversity  digital  radio  link  in 
parallel  with  a  channel  probe  such  as  the  ITS  apparatus.  The  radio  link  would 
be  instrumented  to  measure  digital  error  rate  and  signal  level  in  each  tenth 
second  interval.  Pre-processing  would  discard  data  for  uninteresting  intervals 
while  the  interesting  data  would  be  saved  on  tape  for  later  analysis.  The 
channel  probe  would  measure  delays  on  the  path  fur  correlation  with  errored 
intervals.  A  proposed  site  is  a  45  mile  path  from  Waltham,  MA  to  Pack  Monadnock, 
N.H.  A  NOAA  weather  radar  looks  down  on  the  path  from  Worcester,  MA;  this  aids 
in  identifying  outages  caused  by  precipitation.  The  tests  should  be  accompanied 
by  radio  sonde  flights  during  periods  of  observed  multipath  degradation.  Such 
a  test  would  add  a  great  deal  to  the  body  of  knowledge  as  to  the  interrelation¬ 
ship  of  the  atmosphere,  delays  in  the  medium  ana  digital  radio  performance. 

These  recommendations  for  tests  were  made  in  December  1978.  The  results 
should  be  important  to  both  the  Army  who  is  procuring  the  DRAMA  radios  and  who 
«'h-.  'il cl  speak  for  the-ir  1  inVli  tiers  as.  well  as  agencies  who  apply  the  radio  in 
s /stems .  The.  tft-’fs  have  not  happened  fo-  r.l:  i'-'ast  two  reasons:  -he  radios 
bsV.  j  ci  e  all  oca  lad  c-.»  teals  at  ft.  Iiuu'.Hj.'?  or«l  to  the  production 
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installation  schedule.  Existin'!  program  direction  ,-euld  have  to  be  changed 
to  allocate  radios  to  the  proposed  tests  and  provide  funding  to  a  designated 
testing  agency. 

k  there  has  been  a  report  of  a  n ewer 
test  by  Mr.  Giuffrida  of  BTL  at  ICC  79.  The  test  methods  were  similar  to 
those  reported  earlier  by  Barnett,  the  equipment  and  location  were  the  same 
and  probably  much  of  the  instrumentation.  The  tests  used  three  periods  during 
a  fading  season,  each  greater  than  a  month,  to  measure  cumulative  outage  with 
combinations  of  switching,  coherent  IF  combining  and  adaptive  equalization. 

This  graph  plots  together  Giuffrida 's  separately  plotted  test  results.  The 
system  that  provided  the  lowest  outage  was  the  lowest  black  dashed  line,  it  is 
functionally  the  same  as  the  system  that  BNR  selected  for  production  after 
their  1977  tests.  The  highest  red  line  is  functionally  like  the  DRAMA  radio; 
space  diversity  with  baseband  switching.  There  are  many  questions  that  can 
be  asked  by  those  who  would  choose  to  discount  the  significance  of  this  data; 
nevertheless  it  is  a  comparison  of  six  diversity  techniques  that  suggests  the 
technique  embodied  in  the  DRAMA  radio  to  provide  among  the  poorest  diversity 
improvement  of  those  techniques  tested.  One  can  say  the  test  might  have  been 
done  wrong,  or  that  adaptive  equalization  is  a  requirement  for  either  switching 
diversity  or  IF  combining,  or  the  data  can  support  that  tests  and  work  are  needed 
to  determine  how  to  obtain  significant  diversity  improvement  by  such  means  as 
antenna  installation  criteria  where  radios  are  constrained  to  dependence  on 
switching  diversity. 


(12 n  Since  the  publication  of  this  wor 


DIGITAL  RADIO  MULTIPATH  DEGRADATION  MECHANISM 


WE  NUMBER  OF  SECONDS  OF  OUTAGE  (BER  >  10  ')  PER  WORST-FADING  MONTH: 


DISTANCE -STATUTE  MILES 


CUMULATIVE 


2- 


OUTAGE  VERSUS  OUTAGE  DURATION  FOR  ALL  FADE  DEPTHS 


CUMULATIVE  OUTAGE 


AimavuvAv 
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LOG  BER  IN  ONE  SECOND  INTERVALS 


AMPLITUDE  DISTORTION  MEASUREMENTS  ON  LOS  LINKS 


Larry  Hause 

Institute  of  Telecommunications  Sciences 
National  Telecommunications  and  Information 
Administration 
Boulder,  CO 


NOTE:  The  materials  provided  in  this  section  are  copies  of  the 
visuals  used  with  the  extemporaneous  remarks  given  by 
Mr.  Hause. 
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AMPLITUDE  DISTORTION  MEASUREMENTS 
ON  LOS  LINKS 

1.  REVIEW  PREVIOUS  ITS  WORK  AT  MT.  CORNA,  ITALY. 

2.  ENUMERATE  THE  PURPOSES  OF  THE  PROPOSED 
AMPLITUDE  DISTORTION  MEASUREMENTS. 

3.  DESCRIBE  WHAT  MEASUREMENTS  WILL  BE  MADE. 

A.  DESCRIBE  PROPOSED  DATA  ANALYSIS. 
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MEASUREMENTS  THAT  ITS  HAS  BEEN  MAKING  ON 
THE  CONVERGING  PATHS  AT  MT.  CORN A 

1.  1  YEAR  OF  RECEIVED  SIGNAL  LEVEL  DATA. 


2.  DIVIDE  THE  DATA  INTO  5  CATEGORIES. 


3.  ONE  15  GHz  CHANNEL  ON  EACH  PATH. 


A.  TWO  8  GHz  CHANNELS  ON  EACH  PATH  EXCEPT 
PAGANELLA  -  CORNA. 


5.  OBTAIN  FADING  DEPTH  DISTRIBUTIONS  FOR  EACH 
CATEGORY. 


6.  SHOW  THE  PERIODS  THAT  EACH  FADING  MODE 
CATEGORY  IS  ACTIVE. 


RADIO  RADIO 


.  Mt.  Venda  to  Mt.  Corna  8  GHz  primary 
1979. 


Path  Length  (GCD)  =  90.24  km 
Carrier  Frequency  =  8.32  GHz 
Horizontal  Polarization 


Primary- 


Fraction  of  Time,  P,  Below  Level 

Terrain-  Type  -  average  Vet  Radons  Loss  per  Ant.  =  1. 7 


Avg.  Annual  Temp.  =  27  °C 


Mean  Annual  Rainfall  =  800  mm 
No.  of  Rainy  Daye/yr.  =111 
Thunderstorm  Ratio  B  =  0. 15 


Multipath  and  rain  loss  variability  for  the 
path  from  Mt.  Venda  tower  to  Mt.  Corna  tower 
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ITS  MEASUREMENTS  AT  CEGGIA 

1.  PATHS  CONVERGING  AT  CEGGIA  ARE  MONITORED. 

2.  THESE  PATHS  ARE  MONITORED  ONLY  ON  8  GHz. 

3.  THE  CEGGIA  PATHS  DO  NOT  MEET  THE  DESIRED 
CLEARANCE  CRITERIA  FOR  LOS  PATHS. 
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WEATHER  STATISTICS  GATHERED  FOR  ITS 
BY  AFCS 


1.  DATA  FROM  SEVEN  WEATHER  STATIONS. 


2.  DATA  FROM  MORE  STATIONS  ARE  AVAILABLE  BUT  NOT 
APPLICABLE. 


3.  THE  DATA  IS  PRIMARILY  TEMPERATURE  AND  RAIN 
STATISTICS. 
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PURPOSES  OF  THE  AMPLITUDE 
DISTORTION  MEASUREMENTS 


1.  CAN  WE  DETECT  THE  DISTORTION? 

2.  HOW  OFTEN  DOES  IT  OCCUR? 

3.  WHAT  IS  THE  RANGE  OF  SEVERITY? 

A,  WHAT  EFFECT  DOES  THE  DISTORTION  HAVE 
ON  THE  SYSTEM? 
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THE  SYSTEM  AND  TYPES  OF  PATHS 


1.  ITS  WILL  MONITOR  THE  DEB  SYSTEM  IN  NORTH  CENTRAL 
ITALY. 


2,  IT  IS  A  12  Mb/s  SYSTEM  WITH  A  1A  MHz  BANDWIDTH 
AND  8  GHz  CARRIERS. 


3.  IT  USES  THREE  LEVEL  PARTIAL  RESPONSE 
MODULATION. 


A.  THERE  ARE  THREE  LINE-OF-SIGHT  PATHS  WHICH 
TERMINATE  AT  THE  RECEIVING  AND  RECORDING 
SITE. 


5.  TWO  OF  THE  PATHS  ARE  ABOUT  90  km  AND  ONE  IS 
ABOUT  130  km  LONG. 


DATA  GATHERING 


1.  THE  SPECTRUM  FROM  BOTH  THE  PRIMARY  AND  DIVERSITY 
RADIOS  WILL  BE  MONITORED, 


2.  THE  SIGNALS  WILL  BE  MONITORED  PRIMARILY 
DURING  PERIODS  OF  MULTIPATH  FADING. 


3.  THE  SIGNALS  WILL  BE  MONITORED  OVER  A  PERIOD 
OF  TWO  MONTHS  ON  VARIOUS  OF  THE  THREE  PATHS, 


4.  THE  IF  SPECTRUM  WILL  BE  MONITORED  AT  A  TEST 
POINT  AFTER  AGC  ACTION. 


RECORDER  CHANNELS 


1.  TIME  CODE. 

2.  SPECTRUM  ANALYZER  SAW  TOOTH. 

3.  SPECTRUM  WITH  FREQUENCY  MARKER. 

4.  RECEIVED  SIGNAL  LEVEL. 

5.  IDENTITY  OF  RECEIVER  ON  LINE. 

6.  FRAME  VIOLATION  EVENTS. 


7.  FORMAT  VIOLATION  EVENTS. 


DATA  ANALYSIS 


1.  ITS  WILL  OBTAIN  SPECTRAL  DENSITY  GAIN 
DIFFERENCE  VALUES  IN  Db  FOR  CORRESPONDING 
FREQUENCIES. 


2.  A  STANDARD  SPECTRUM  OBTAINED  DURING  A 
QUIET  RSL  PERIOD  WILL  BE  USED  FOR 
COMPARISON. 


3.  ALL  SIGNAL  DIGITIZING  THAT  IS  REQUIRED 
WILL  BE  DONE  AT  ITS  LABS  IN  BOULDER,  CO. 


A.  DATA  SEARCHES  WILL  BE  MADE  ON  BOTH  THE  BASIS 
OF  FRAME/FORMAT  EVENTS  AND  MAGNITUDE  OF 
SPECTRAL  DENSITY  DIFFERENCES. 
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Spectrum 


T  ime 
Code 


Analog  Magnetic 
Tape  Recorder 
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DATA  CLASSIFICATION  SPECIFICATIONS 

The  purpose  of  these  specifications  is  to  categorize  the 
RSL  data  within  each  hour  in  terms  of  the  various  types  of 
mechanisms  that  are  operating  on  the  radio  path  or  possible 
adverse  equipment  situations  that  prevailed  within  the  hour. 
These  specifications  are  used  to  prepare  a  data  analysis  log 
so  that  data  contaminated  by  equipment  malfunction  can  be 
removed  and  so  that  data  having  the  characteristics  typical  of 
certain  fading  mechanisms  are  properly  classified  and  can  be 
checked  against  meteorological  data  and  site  operations  log 
entries.  The  log  is  needed  in  order  that  signal  level  time 
distributions  can  be  prepared  for  the  various  propagation 
mechanisms  of  interest. 

Six  types  of  received  signal  level  recording,  hourly 
period  classifications  will  be  used: 

1.  (Equipment  Malfunction,  EM) 

Communications  or  test  equipment  malfunction  such  that 
data  taken  on  that  particular  channel  should  not  be 
used. 

2.  (Undisturbed  Condition,  UC) 

Essentially  undisturbed  signal  level  period  during 
which  the  signal  does  not  fade  even  momentarily  more 
than  10  dB  below  the  long  term  median. 

3.  (Multipath  with  High  Median,  MHM) 

A  period  during  which  multipath  fading  occurs  (greater 
than  10  dB  below  the  median)  but  the  hourly  median 
signal  level  does  not  fade  more  than  6  dB  below  the 
long  term  median  signal  level. 

4.  (Multi-path  Fading  with  Power  Fading,  M&P) 

A  period  during  which  multipath  fading  occurs 
(greater  than  10  dB  below  the  hourly  median)  but  the 
hourly  median  is  also  depressed  by  more  than  6  dB 
below  the  long  term  median. 

5.  (Rain  Attenuation  Fading,  RA) 

A  period  of  rain  attenuation  fading  (greater  than 
10  dB  below  the  long  term  median  at  some  time  during 
the  hour) .  Rain  attenuation  fading  leaves  a 
distinctive  trace  on  the  chart  recorder  and  also 
produces  highly  correlated  recordings  between  the 
primary  and  space  diversity  channels. 

6.  (Power  Fading  without  Rain  or  Multipath,  PF) 

A  period  of  depressed  median  signal  level  (greater 
than  6  dB)  but  without  multipath  or  rain  attenuation 
characteristics. 
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4.  OVERVIEW  OF  MULTIPATH  MEASUREMENTS  USING  A  PN  PROBE 


Bob  Hubbard 

Institute  of  Telecommunications  Sciences 
Boulder,  CO 


NOTE:  The  materials  provided  in  this  section  are  copies 
of  the  visuals  used  with  the  extemporaneous  remarks 
given  by  Dr.  Hubbard 
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THE  WIENER  RELATIONSHIP 


1.  FREQUENCY  TRANSFER  FUNCTION 

*y(«)  =  |H0“)  I  2  ^(W) 


*  («)  =  HOW)  ^(w) 


2.  IMPULSE  RESPONSE  FUNCTION 

Rxy(T)  =  h<T>  ®  Rxx(T) 

Since  h(t)  =  F  |H(jcc)l 

R»(T)  * F  {*»<“»} 

V  - F  {%<“'} 

f/  1  denotes  a  Fourier  transform. 


Where 


Figuat  IS.  Maximum  length  ihi^t  ae.gi.4tta  itqutnce  ehaaacteaiA tic* 


RESPONSE  AMPLITUDE 


L*J  p<rm*r  L*J  Lof 

Att.  (^1  Meter  I  ^1  Amplifle 


RECEIVER  FOR  8.  6  GH*  CHANNEL  PROBE 


Holeokolo  \  Wouno  100 


S 


geometries. 


Figure.  17, 


TIME  DELAY  IN  NANOSECONDS 

Impulse.  A.e-&pon4&  between  Ha.Zzaka.Za  and  Kona,  Hawaii, 
MaAC.fr  /S,  1974  at  2  izcond  intzx.val& . 
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197 


Fiqure  6.  Sample  of  the  fading  records  over  the  Keahole-Haleakala  path  on 
24  March  1974. 


lUgust  1974 


5 


■ - STANDARD  ATMOSPHERE  -  i - 


(a)  Refractivity  profile  measured  at 
1700  CST  on  18  November,  1976. 


- MI>.J  ' 

(c)  Impulse  response  during  a  15  dB  fade. 

Figure  12.  Refractivity  profile  and  power  impulse  responses 

recorded  during  a  ducting  propagation  conditon  along 
the  Gulf  Coast  of  Florida. 


i 
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(b)  Time-lapse  photograph  of  2 Q  s  period  at  1256  PST. 


Figure  10(a).  The  average  power  impulse  response  measured  between 

1255  and  1300  PST,  31  August  1978,  at  a  PN  clock  rate 
of  50  MHz. 


(b)  Time-lapse  photograph  of  20  s  period  at  1304  PST. 

Figure  11(a).  The  average  power  impulse  response  measured  between 

1300  and  1305  PST,  31  August  1978,  at  a  PN  clock  rate 
of  50  MHz. 
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(a)  Time-lapse  response  measured  at  1303  PST,  corresponding  to  the 
fade  in  Ch.  1  and  the  loss  of  synchronization  in  the  data  stream 
shown  in  Figure  12 (b) . 


(b)  Time-lapse  response  measured  at  1314  PST,  where  the  data  stream 
resynchronized. 


(c)  Time-lapse  response  measured  at  1310  PST,  where  the  data  stream 
resynchronized  for  a  few  seconds. 

Figure  14.  Time-lapse  photographs  of  the  power  impulse  response 
measured  during  the  fading  period  in  Ch.  1  beginning 
at  1308  PST,  31  August  1978. 


(a)  Dynamic  .multipath  measured  in  Ch.  2  (high  antenna)  at  1317  PST. 


(b)  A  flat-fade  response  observed  in  Ch.  1  (low  antenna)  at  1331  PST. 


(c)  Another  example  showing  a  power  shift  between  components  in  Ch.  1, 
with  only  minor  change  in  the  Ch.  2  response  at  1346  PST. 


Figure  17.  Time-lapse  photographs  of  the  power 
measured  between  1317  and  1346  PST, 


impulse  response 
31  August  1978. 
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5.  PATH  TESTING  OF  LONG,  OVER  WATER  DIGITAL  RADIO  PATHS 


Jim  Weblemoe 

Pacific  Missile  Test  Center 
Pt.  Mugu,  CA 


NOTE:  The  materials  provided  in  this  section  are  copies 
of  the  visuals  used  with  the  extemporaneous  remarks 
given  by  Mr.  Weblemoe. 
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OFFSHORE  MICROWAVE  PATHS 


031 v 


AD-A087  793  DEFENSE  COMMUNICATIONS  ENGINEERING  CENTER  RESTON  VA  F/G  17/2.1 

PROCEEDINGS  OF  SEMINAR  ON  FREOUENCY  SELECTIVE  FADING  AND  ITS  EF— ETC<U) 
JAN  80  0  R  SMITH 

UNCLASSIFIED  pCEC-TN-2-QO _ SBIE-AD-E100  375 _ NL _ 


Flyura  5(a).  Rafractiva  indtx  prof 11a  for  Ft.  Mugu  at  1058  VST,  31  August  1978 


6.  MULTIPATH  FADING  CHANNEL  MODEL 


Tom  Gluffrida 

Bell  Telephone  Laboratories 
Holmdel ,  NJ 


NOTE:  The  materials  provided  in  this  section  are  copies 
of  two  Bell  System  Technical  Journal  articles  used 
as  the  basis  for  the  extemporaneous  remarks  given 
by  Dr.  Giuffrida. 
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Copyngbt  ©  1979  Amarican  Talaphooa  and  Taiagraph  Company 
The  Bell  System  Technical  Journal 
VoL  56.  No.  5,  May-Juna  1979 
Pruutdi a 


A  New  Selective  Fading  Model:  Application  to 
Propagation  Data 

By  W.  D.  RUMMLER 
(Manuscript  received  September  21 , 1 978) 

Channel  transmission  models  for  use  in  estimating  the  performance 
of  radio  systems  on  line-of-sight  paths  at  6  GHz  are  explored.  The 
basis  for  this  study  is  the  simple  three-ray  multipath  fade,  which 
provides  a  channel  transfer  function  of  the  form  H(u)  -  afl  -  b  exp 
-j(i c  —  uo)r],  where  a  is  the  scale  parameter,  bis  a  shape  parameter, 
r  is  the  delay  difference  in  the  channel,  and  <*>  is  the  (radian) 
frequency  of  the  fade  minimum.  This  model  is  indistinguishable  from 
an  ideal  channel  model,  within  the  accuracy  of  existing  measure¬ 
ments.  The  propagation  data  that  confirm  the  model  were  obtained 
in  summer  1977  from  a  26.4-mile  hop  near  Atlanta,  Georgia.  The 
received  power  at  24  sample  frequencies  spaced  at  1.1  MHz  and 
centered  on  60343  MHz  was  continuously  monitored  and  recorded 
during  periods  of  anomalous  behavior.  The  model  is  applied  to 
estimating  the  statistics  of  the  channel  delay  difference,  r.  The  aver¬ 
age  delay  difference  giving  rise  to  significant  selectivity  in  the  channel 
is  between  5  and  9  ns.  The  distribution  of  delay  difference  is  obtained 
for  delay  differences  greater  than  10  ns.  The  channel  is  found  to  have 
more  than  3  dB  of  selectivity  (difference  between  maximum  and 
minimum  attenuation  in  band)  due  to  delay  differences  greater  than 
20  ns  for  more  than  70  seconds  in  a  heavy  fading  month.  (This  is 
comparable  to  the  time  the  channel  attenuation  of  a  single  frequency 
exceeds  40  dB.)  The  three-path  model  requires  further  simplification 
for  narrowband  channel  application.  For  a  channel  with  30  MHz 
bandwidth,  a  model  with  fixed  delay  of  6.3  ns  provides  a  sufficiently 
accurate  representation  of  all  observed  channel  conditions.  The  re¬ 
sulting  nonphysical  model  is  used  to  statistically  characterize  the 
condition  of  the  fading  channel  The  statistics  of  the  parameters  of 
the  fixed  delay  model  are  almost  independent  and  of  relatively  simple 
form.  The  distribution  of  the  shape  parameter  b  is  of  the  form  (1  - 
bfr*.  The  distribution  of  a  is  lognormal  For  b  >  0.5,  the  mean  and 
standard  deviation  of  -20  (log  a)  are  25  and  5  dB,  respectively;  the 
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mean  decreases  to  15  dB  for  smaller  values  of  b.  The  probability 
density  function  of  wo  is  uniform  at  two  levels;  measuring  wo  from  the 
center  of  the  band,  the  magnitude  of  wot  is  five  times  as  likely  to  be 
less  than  n/2  than  to  be  greater.  A  companion  paper  describes  the 
use  of  this  model  for  determining  the  bit  error  rate  statistics  of  a 
digital  radio  system  on  the  modeled  path. 

I.  INTRODUCTION 

Performance  prediction  of  a  digital  radio  system  on  a  line-of-sight 
microwave  channel  requires  an  accurate  statistical  model  of  the  chan¬ 
nel  Because  different  digital  radio  systems  may  have  different  sensi¬ 
tivities  to  the  various  channel  impairments,  the  model  must  be  com¬ 
plete  to  the  extent  that  it  must  be  capable  of  duplicating  the  amplitude 
and  phase  (at  least  approximately)  of  all  observed  channel  conditions. 
To  facilitate  laboratory  measurements  and  computer  simulations  for 
calculating  outage,  the  model  should  be  realizable  as  a  practical  test 
circuit  and  should  have  as  few  parameters  as  possible.  Most  important, 
the  parameters  should  be  statistically  well  behaved. 

Two  types  of  models  have  been  generally  considered  for  line-of-sight 
microwave  radio  channels:  power  series  type  models1-3  and  multipath 
models.*"*  A  power  series  model  will  require  a  few  terms  only  if  the 
channel  is  a  multipath  medium  with  a  small  spread  of  delays  relative 
to  the  reciprocal  bandwidth  of  the  channel3  This  implies  that  one 
must  understand  the  channel  as  a  multipath  medium  to  understand 
the  behavior  of  a  power  series  model  Hence,  we  have  limited  our 
characterization  efforts  to  multipath  models. 

The  basis  for  this  study  is  the  simple  three-ray  multipath  fade.7  If 
the  fading  in  a  channel  can  be  characterized  by  a  simple  three-path 
model,  the  channel  will  (as  shown  in  Section  Q)  have  a  voltage  transfer 
function  of  the  form 


H(w)  -  a[l  -  be*'"’-’-""].  (1) 

where  the  real  positive  parameters  a  and  b  control  the  scale  and  shape 
of  the  fade,  respectively,  r  is  the  delay  difference  in  the  channel  and 
wo  is  the  radian  frequency  of  the  fade  minimum.  The  plus  and  minus 
signs  in  the  exponent  correspond,  respectively,  to  the  channel  being  in 
a  nonminimum  phase  or  minimum  phase  state.  Note  that,  with  appro¬ 
priate  choices  of  parameters,  this  model  can  be  reduced  to  a  two-path 
model  or  a  scaled  two-path  model  etc. 

It  has  been  shown  previously,7  and  is  illustrated  in  Section  Q,  that 
the  simple  three-path  fade  overspecifies  the  channel  transfer  function 
if  the  delay  is  less  than  H  B,  where  B  is  the  observation  bandwidth. 
The  critical  value  of  r  for  a  30- MHz  channel  is  about  5.3  ns.  which  is 
comparable  to  the  mean  delay  in  the  channel  As  a  consequence,  unless 
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the  channel  response  can  be  determined  to  an  accuracy  on  the  order 
of  0.001  d£,  a  unique  set  of  parameters  a,  b,  r,  and  f0  cannot  be 
determined  for  more  than  half  the  faded  channel  conditions  encoun¬ 
tered.  To  avoid  this  problem,  one  must  suppress  or  fix  one  of  the 
model  parameters.  Section  Q  shows  that  the  delay,  r,  is  the  only 
parameter  which,  when  fixed,  produces  a  reasonable  model. 

While  a  model  with  a  fixed  delay  may  appear  to  be  a  strange  choice, 
it  has  ail  the  required  characteristics  for  modeling  the  channel  transfer 
function.  Figure  1  shows  the  amplitude  of  the  channel  transfer  function 
of  eq.  (1)  on  a  power  scale  and  on  a  decibel  scale  for  r  <■>  6.31  ns.  With 
r  fixed,  the  response  minimum  is  shifted  with  respect  to  frequency  by 
varying  /».  Varying  a  changes  the  overall  level  and  b  changes  the 


Fig.  1— Channel  mod*!  function.  Jf(u> m  a[l  -  bnpl-j  2  f,)r)],torr  »6.3  m. 
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“shapeliness.”  If  the  minimum  is  within  the  30-MHz  bandwidth  of  a 
channel,  the  fixed  delay  model  can  generate  notches  with  a  wide  range 
of  levels  and  notch  widths.  With  the  minimum  out  of  band,  it  can 
generate  a  wide  range  of  combinations  of  levels,  slopes,  and  curvatures 
within  the  channel  bandwidth.  Section  VI  shows  that  the  model 
versatility,  with  r  chosen  to  be  6.31  ns,  is  sufficient  to  characterize  a 
30-MHz  channel  in  the  6-GHz  common  carrier  band. 

Section  II  provides  a  brief  discussion  of  the  simple  three-path  fade. 
A  comparative  discussion  of  the  relative  merits  of  the  different  possible 
simplifications  of  this  model  leads  to  the  choice  of  the  fixed  delay 
model. 

The  data  used  for  detailed  evaluation  of  models  were  obtained  from 
a  6-GHz  experiment  in  Palmetto,  Georgia,  in  June  1977.  The  radio 
channel  was  equipped  with  a  general  trade  78-Mbit/s,  8-PSK  digital 
radio  system,  and  the  received  spectrum  was  monitored  with  a  set  of 
24  filters  with  band  widths  of  200  kHz  spaced  at  a  1.1 -MHz  separation 
across  this  channel  During  fading  activity,  the  received  power  of  each 
of  these  frequencies  was  measured  five  times  each  second,  or  once 
every  2  seconds,  depending  on  how  rapidly  the  channel  was  changing; 
sampled  power,  quantized  in  1-dB  steps,  was  recorded  by  the  midas 
system.*  The  data  base  used  for  this  study  consists  of  approximately 
25,000  scans  representing  8400  seconds  of  fading  activity;  about  8700 
scans  were  recorded  during  periods  when  the  equipment  was  indicating 
errors.  These  data  represent  about  60  percent  of  the  fading  activity  of 
a  heavy  fading  month;  therefore,  the  derived  statistics  must  be  viewed 
as  provisional  and  subject  to  some  modification  as  additional  data  are 
processed.  At  the  very  least,  the  data  base  is  sufficiently  large  to 
indicate  what  can  happen  on  the  channel  and  to  form  a  basis  for 
choosing  and  validating  a  model 

As  described  in  Section  m,  the  model  parameters  were  estimated 
for  each  scan  by  fitting  the  magnitude  squared  of  the  transfer  charac¬ 
teristic  [eq.  (1)]  to  the  observed  channel  shape  as  characterized  by  the 
power  received  at  the  sampling  frequencies.  Phase  is  subsequently 
derived  by  assuming  the  channel  is  minimum  phase.  Problems  are 
encountered  in  realizing  a  minimum-phase  solution  because  of  quan¬ 
tization  noise  and  the  presence  of  certain  channel  shapes  caused  by 
large  delays.  The  procedure  for  handling  these  difficulties  is  described. 

The  statistics  of  the  parameters  of  the  fixed  delay  model  are  dis¬ 
cussed  in  Section  IV.  Equations  providing  an  idealized  description  of 
the  statistics  of  the  parameters  of  the  model  are  also  given  here. 

In  Section  V,  the  determination  of  the  delay  difference  present  in 
the  channel  is  considered.  In  the  first  subsection,  it  is  demonstrated 


*  Multiple  Input  Diu  Acquisition  System,  constructed  by  G.  A.  Zimmerman:  see 
Ref.  1. 
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that,  during  the  observed  period  of  fading  activity,  the  average  delay 
is  9  ns.  A  lower  bound  on  the  distribution  of  delay  difference  for  large 
delays  is  developed  in  the  second  subsection.  A  third  subsection 
provides  an  example  of  a  channel  scan  that  can  best  be  approximated 
by  a  three-path  fade  with  a  delay  difference  of  26  ns.  Fades  with  at 
least  this  delay  and  with  a  more  moderate  amount  of  shape  (2  dB  or 
more)  were  encountered  for  about  60  seconds  of  the  data  base  studied. 
Thus,  one  might  expect  26-ns  delays  to  be  present  during  about  100 
seconds  of  a  heavy  fading  month. 

The  presence  of  such  large  apparent  delays  raises  questions  as  to 
the  accuracy  with  which  the  fixed  delay  model  represents  the  channel. 
These  questions  are  addressed  in  Section  VI  where  the  statistics  of  the 
errors  in  modeling  scan  fits  are  described.  The  errors  are  small  and  do 
not  compromise  the  usefulness  of  the  modeL 

Results  and  conclusions  are  briefly  summarized  in  Section  VTI. 


II.  CHOICE  OF  MODEL 

In  this  section,  we  provide  a  brief  description  of  the  simple  three- 
path  model  and  show  why  it  cannot  be  used  to  estimate  delays  when 
the  delay  bandwidth  product  is  less  than  <4.  In  a  comparative  discus¬ 
sion,  we  show  why  the  fixed  delay  model  is  the  only  simplification  of 
the  model  that  is  manageable. 

Z1  Simpim  thrm+umth  mtxM 

Consider  a  channel  characterized  by  three  paths  or  rays.  The  am¬ 
plitude  of  the  signal  on  each  of  these  three  paths,  as  seen  by  the 
receiver,  is  1,  at,  and  a*.  The  second  and  third  paths  are  delayed  with 
respect  to  the  first  by  rt  and  r2  seconds,  respectively,  where  t2  >  n.  We 
define  the  simple  three-path  model  by  requiring  the  delay  between  the 
first  two  paths  to  be  sufficiently  small,  i.e., 

(w?  -  Wl)Tl  <K  l,  (2) 

where  u?  and  «t  are  the  highest  and  lowest  (radian)  frequencies  in  the 
band.  The  complex  voltage  transfer  function  of  the  channel  at  a 
frequency  u  may  be  illustrated  with  a  phasor  diagram.  Figure  2a  shows 
the  phasor  diagrams  for  wi  and  superimposed.  By  designating  the 
amplitude  of  the  (vector)  sum  of  the  first  two  paths  by  a;  the  angle  of 
the  sum  by  <p  ■»  wot  -  n,  where  r  is  equal  to  r*.  the  delay  difference  in 
the  channel;  and  the  amplitude  of  the  third  ray  by  ab,  we  obtain  the 
simplified  diagram  in  Fig  2b.* 


*  Not*  that,  if  the  third  amplitude  is  greater  than  the  sum  of  the  first  two.  wt 
iaterchani*  the  assignments  of  amplitudes  a  and  ab  and  obtain  a  nonminimum  phase 
fade. 
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midband  and  at  both  edges.  It  has  been  shown  elsewhere  that  fades 
matched  in  this  way  will  be  identical  in  band  to  within  a  few  tenths  of 
a  decibel  at  most,  and  will  have  almost  identical  envelope  delay 
distortion.  Given  noisy  quantized  measurements  of  |  H( w>  I  over  the 
band,  it  is  impossible  to  distinguish  between  such  fades  unless  we  fix 
one  of  the  four  parameters.  Let  us  consider  each  of  the  four  possibili¬ 
ties. 


2.2  Pseudo  two-path  tad  a 

If  one  fixes  the  amplitude,  a,  at  unity,  the  simple  three-path  fade 
reduces  to  a  two-path  fade  with  independent  control  of  the  frequency 
of  the  minimum  of  the  response.  The  difficulty  with  this  model,  as 
may  be  seen  by  referring  to  Fig.  2b,  is  that  it  can  provide  in-band 
minima  only  for  |H(w)  |  <  1  and  maxima  in-band  only  for  \H{u)  \  >  1. 
In  other  words,  the  model  cannot  match  an  in-band  maximum  at  an 
arbitrary  fade  leveL  In  addition,  it  was  found  that  during  approximately 
half  the  periods  when  the  radio  equipment  was  indicating  errors,  the 
channel  could  not  be  well  modeled  with  a  pseudo  two-path  model. 

2.3  Scaled  two-path  fade 

If  one  fixes  the  phase.  <b  -  uo-r  -  ir,  in  the  simple  three-path  model 
at  0,  the  fade  reduces  to  a  scaled  two*path  fade.  (For  a  two-path  fade, 
we  require  the  additional  condition  a  *  I.)  This  is  the  most  physically 
desirable  of  the  reduced  three-path  models  because  it  may  be  derived 
without  recourse  to  the  three-path  formalism.  Unfortunately,  it  is 
mathematically  intractable,  particularly  when  dealing  with  amplitude 
data  only.  In  fitting  the  model  to  a  given  channel  shape  (in  the  manner 
described  in  Section  HI  for  the  fixed  delay  model),  one  obtains  a 
function  of  a,  b,  and  r  that  must  be  minimized  to  obtain  the  best  fit. 
Because  of  the  ur  term  in  the  exponent  of  the  model,  this  function  has 
a  local  minimum  in  every  interval  of  r  of  length  0.17  ns,  the  reciprocal 
of  6  GHz.  Since  the  possible  range  of  ~  extends  to  about  30  ns,  one 
may  have  to  perform  hundreds  of  minimizations  to  find  the  best  fit  to 
a  single  channel  scan.  Even  then  this  “best  fit”  may  have  no  minimum 
phase  realization,  and  there  is  no  known  procedure  that  leads  to  one. 

2.4  Fixed  b  modal 

If  one  fixes  the  amplitude  b  in  the  simple  three-path  model,  the 
resulting  reduced  model  has  all  the  mathematical  difficulties  of  the 
scaled  two-path  model  and  no  satisfactory  physical  interpretation. 

2.3  Fixed  delay  model 

It  is  demonstrated  in  the  remainder  of  this  paper  that  the  fixed 
delay  model  described  in  Section  I  is  useful  and  effective  in  character¬ 
izing  the  channel. 
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III.  ESTIMATION  PROCEDURES 

This  section  describes  how  the  model  parameters  are  estimated  from 
the  channel  scans  and  how  realizability  difficulties  are  surmounted. 

3.1  Parameter  animation 

The  channel  data  consist  of  a  set  of  25,000  scans  of  the  channel 
power  spectrum.  Each  scan  consists  of  a  power  measurement  at  each 
of  24  frequencies  at  1.1-MHz  spacing  across  the  channel.  (Actually, 
only  23  frequencies  are  used  since  the  19th  was  inoperative  during  this 
test  period).  The  power  measurements  are  recorded  in  decibels,  and 
each  must  be  referenced  to  the  average  power  level  of  that  frequency 
at  mid-day.  With  proper  conversion  and  calibration,  the  basic  data 
characterizing  a  scan  are  a  set  of  power  ratios.  We  designate  the  power 
ratio  at  nth  frequency  by  V,,  where 

Yn  -  y(w„)  n  »  1,  2,-  •  •,  24.  (3) 

We  wish  to  model  the  channel  with  a  voltage  transfer  function  of 
the  form  given  in  eq.  (1),  which  we  repeat  here  for  convenience 

HM  -  a(l  -  6eiy,u-*,'T],  (1) 

Thus  our  estimate  of  Yn  will  be 

Yn  “  I  H (w* )  | 2  •*  a  -  y3  Ccs(w«  -  Uo)t,  (4) 

where 

a  ■  a2(l  +  b 2) 

0  -  2 a2b.  (5) 

For  convenience,  we  measure  frequency  in  the  units  of  the  frequency 
separation  of  the  power  measurements.  Thus, 

-  2 itfn  -  2ffn(l.l  x  10s)  n  -  1,  2,  3,  •  •  •,  24.  (6) 

If  we  choose 

1 

T  “  N(l.l  x  10")’ 

then 


-  2jt  (8) 

N 

For  the  fixed  delay  model,  we  choose  iV  -  144  which  gives  a  model  r 
of  6.31  ns.  Thus,  the  in-band  frequencies  correspond  to  n  values 
between  1  and  24,  and  the  channel  transfer  function  given  by  the 
model  is  periodic  for  n  modulo  144,  corresponding  to  a  frequency  shift 
of  144  x  1.1  x  10®  -  158.4  MHz. 
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The  weighted  mean-square  error  between  the  estimated  and  ob¬ 
served  power  is  given  by 


24 


S  cay,  -  Y a2 


l  c* 

n»l 


where  the  summation  skips  n  —  19  as  described  above,  and  where  C„ 
is  a  weighting  applied  to  the  measurement  at  frequency  .  Since  the 
original  data,  from  which  the  YVs  were  derived,  were  uniformly  quan¬ 
tized  on  a  logarithmic  scale,  it  is  appropriate  to  use  a  weighting  that  is 
approximately  logarithmic.  Hence,  we  use  the  weighting  function 


G. 


(10) 


A  number  of  different  weighting  functions  were  tested,  but  the  one 
given  by  (10)  is,  generally,  the  most  satisfactory. 

Estimates  of  a,  b,  and  f0  may  be  obtained  by  minimizing  the  weighted 
mean-square  error,  E.  It  is  shown  in  the  appendix  that  one  may  obtain 
closed  form  estimators  for  a,  /?,  and  A>  by  substituting  eq.  (4)  into  (9) 
and  minimizing  E,  first  with  respect  to  a.  then  with  respect  to  ft  (or 
vice  versa),  and  last  with  respect  to  f0.  In  the  resulting  scheme,  the 
estimator  of  f0,  the  frequency  of  the  model  minimum,  is  a  function  of 
data  only.  The  estimators  of  a  and  13  are  functions  of  the  estimated  fo 
and  the  data.* 

After  estimates  of  a  and  (3  have  been  calculated,  the  parameters  a 
and  b  of  the  model  are  obtained  by  inverting  the  relationships  given 
by  eq.  (5). 


(11) 

(12) 


It  is  clear  from  (11)  and  (12)  that  we  can  realize  the  channel  shape 
with  the  model  only  if  a  2  /?.  This  is  to  be  expected.  Since  |  H(w)|2  is 
a  power  transfer  function,  it  must  be  positive  for  all  frequencies,  which 
is  possible  only  if  a  2  j3  [see  eq.  (4)].  Thus,  the  condition  a  2  /3  allows 
us  to  obtain  a  minimum  (or  nonminimum)  phase  transfer  function 
whose  magnitude  squared  is  the  minimum  weighted  mean-square  error 
fit  to  the  observed  power  transfer  response  of  the  channel. 


*  For  mathematical  simplicity,  we  actually  usa  an  estimator  for  $  conditioned  on  U. 
a.  and  the  data. 
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3.2  Application  of  estimators 

If  the  procedure  described  above  is  3trictiy  applied  to  the  set  of 
25,000  scans  in  the  data  base,  one  finds  that  about  35  percent  of  the 
scans  cannot  be  modeled  with  real  values  of  a  and  b.  A  study  of  these 
problem  scans  revealed  that  the  estimator  for  f0,  the  frequency  of  the 
fade  minimum,  was  biased  for  two  types  of  scans.  One  type  is  a  scan 
with  little  shape,  dominated  by  quantization  noise;  the  other  is  a 
selective  channel  shape  having  a  steep  slope  across  the  band.  Both 
types  of  scan  are  illustrated  in  Fig.  4.  The  scan  in  Fig.  4,  which  is 
almost  flat,  was  fabricated  to  illustrate  the  severity  of  the  quantization 
problem.  The  other  scan  is  typical  of  the  more  shapely  troublesome 
scans. 

To  obtain  a  good  realizable  fit  to  such  channel  shapes  requires 
degrading  the  quality  of  the  fit;  that  is,  moving  the  parameters  away 
from  the  values  that  minimize  the  fit  error,  eq.  (9).  Given  the  form  of 
the  estimation  scheme,  this  is  easily  accomplished  by  moving  the 
frequency  of  the  fade  minimum,  fa,  away  from  its  original  “optimum” 
value  and  reoptimizing  the  remaining  parameters  to  obtain  values  of 
a  and  b  that  are  optimum  for  the  new  value  of  fa.  Figures  5  and  6 
illustrate  the  results  of  such  a  quasi-optimization  regarding  /„  as  a  free 
parameter.  They  show  the  fit  error  E  and  the  values  of  a  and  b  as  fa  is 
varied  from  its  original  optimum  value.  Figure  5  corresponds  to  the 
flat  fade  in  Fig.  4  and  Fig.  6  to  the  sloped  fade. 

The  shapes  of  the  curves  in  Figs.  5  and  6  are  typical  of  those 
obtained  when  the  channel  has  no  minimum  in  band.  The  weighted 
error  in  the  fit,  E,  is  not  very  sensitive  to  the  estimate  of  f„,  the 


s 


Fig.  5— Locua  of  weighud  fit  error  and  modal  parameters  with  U  aa  a  has  variable  for 
Oat  fad*  in  Fig.  4. 


frequency  of  the  modeled  fade  minimum.  The  minimum  of  £  is  broad 
and  flat,  due  to  quantization  and  instrumentation  noise  in  the  channel. 
The  variation  of  the  parameter  a  with  f„  is  also  typically  very  gradual 
The  salient  features  of  the  variation  of  b  with  respect  to  fo  are  clearly 
seen  in  Fig.  6,  and  are  also  present  and  labeled  in  Fig.  5.  As  f0  is  varied 
from  its  original  optimum  value,  b  varies  from  a  value  of  1  to  a  value 
of  0  in  a  sideways  3-curve  with  two  stationary  points,  a  maximum  and 
a  minimum.  Extensive  simulations  with  known  channel  characteristics 
indicate  chat  a  good  choice  of  parameters  is  the  set  corresponding  to 
the  point  where  b  is  locally  minimized.  To  illustrate  this  point,  assume 
that  the  channel  shape  is  that  given  by  the  model  with  6.3-ns  delay, 
fo  at  18.5  x  1.1  MHz,  a  «■  0.04,  and  b  -  0.7.  One  can  construct  a  plot 
similar  to  Figs.  5  and  6  for  this  simulated  fade,  with  the  result  shown 
in  Fig.  7.  The  curves  in  this  figure  illustrate  the  results  cited  above,  in 
that  the  true  value  of  /<■  occurs  near  a  minimum  value  of  b.  A  better 
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*  1.1  MHz 


FREQUENCY  OF  MOOEL  MINIMUM  MEASURED  FROM  MIOSAND 

Pig.  6 — Locus  of  weighted  fit  error  and  model  parameters  with  f„  hi  free  variable  for 
typical  scan  in  Fig.  4. 

choice  for  the  case  shown  and  for  others  that  have  been  simulated 
would  be  “on  the  shoulder”  between  the  minimum  and  b  *■  1;  however, 
such  a  criterion  is  difficult  to  quantify. 

To  summarize,  if  the  standard  routine  does  not  provide  a  realizable 
fit  to  a  scan,  one  merely  varies  fa,  the  position  of  the  minimum,  until 
one  obtains  a  realizable  solution  with  a  value  of  b  that  is  stationary* 
with  respect  to  variations  in  f„.  We  recognize  that  this  procedure 
introduces  additional  sources  of  error  into  the  estimates  of  the  model 
parameters.  The  errors  in  a  and  b  are  small  because  b  is  near  a 
stationary  value  and  a  is  slowly  varying.  The  error  in  /i,  is  also  small, 

*  Since  6  is  a  monotone  function  of  a/B,  it  is  only  necessary  to  invert  solutions  with 
stationary  vaiusa  of  the  ratio,  a/0. 
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usually  less  than  3  MHz,  but  is  always  in  the  direction  corresponding 
to  moving  the  minimum  nearest  to  the  band  closer.  We  consider  the 
effects  of  these  errors  in  Section  VL 

IV.  MODEL  STATISTICS 

Applying  the  procedures  described  in  Section  III  to  the  scans  in  the 
data  base  results  in  25,000  sets  of  values  of  a,  b,  and  fa.  The  relative 
joint  frequency  of  occurrence  of  these  three  parameters  may  be  de¬ 
scribed  by  the  set  of  distribution  functions  shown  in  Figs.  8  to  12.  The 
distribution  of  the  parameter  b  is  described  in  Fig.  3  in  terms  of  the 
distribution  of  -20  log  (1  -  b),  which  is  approximately  exponentially 
distributed  with  a  mean  of  3.8  dB.  This  distribution  gives  the  time 
that  b  exceeds  the  value  given  by  the  abscissa  as  a  fraction  of  the  time 
in  a  heavy  fading  month  that  the  nns  level  in  the  channel  is  depressed 
by  more  than  15  dB.  For  instance,  we  see  that  40  percent  of  the  time 


,  I  1  MMj 


FREQUENCY  of  model  minimum  measured  from  miobano 

Fig.  7— Locus  of  weighted  fit  error  and  model  parameters  with  f*  at  a  free  variable. 
For  channel  given  by  model  with  r  -  6.31  na,  a  ■  0.04.  b  -  0.7.  -  18.5  x  1.1  MHx. 


NEW  SELECTIVE  FADING  MODEL 

6-14 


-20  Loq  11-0)  IN  DECIBELS 
Fig.  3 — Distribution  of  b. 

when  the  channel  is  depressed  the  value  of  b  exceeds  0.3.  It  exceeds 
0.7  for  4  percent  of  that  time,  and  0.99  about  0.3  percent  of  that  time. 
The  distribution  of  b  can  be  modeled  in  the  form 

20 

P{1- b<X)  •XUU'10 -X™.  (13) 

The  distribution  of  a  is  conditioned  on  b  and  is  approximately 
lognormal  as  shown  in  Figs.  9  and  10.  The  mean  and  standard  deviation 
of  the  distributions  in  Figs.  9  and  10  are  plotted  in  Fig.  11.  From  Figs. 
9  to  11  it  is  apparent  that  a  and  b  are  almost  independent;  however, 
less  shapely  fades  tend  to  occur  at  less  depressed  values.  We  note  that 
shape  occurs  when  the  average  depression  is  20  to  25  dB,*  that  the 


*  Ths  value  of  a  correspond*  to  avsrage  power  level  over  a  large  frequency  span  and 
not  strictly  to  the  avenge  power  in  a  narrowband  channel 
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CURVE  8  >  .  3  v  Time,  mc 


I  6.5  .  903.6  | 

'  8.5  ]  343.6  ; 


A  ■  -20  Log  («) 

.  9— Distribution  of  a  cowiitioood  on  tbo  vaiuo  of  S  for  -20  lot  (1  -  b)  law  than 


average  depression  is  near  25  dB  for  b  greater  than  0.7,  and  that  it  fails 
off  gradually  to  15  dB  for  small  b.  The  distribution  of  A  -  -20  log  a  is 
conditioned  on  b  and  may  be  modeled  as 


P{A  >  Y ) 


■ 

-  1  -P  ■ 


r-A,(5)J 

5 


where  P  is  the  cumulative  distribution  function  of  a  zero  mean,  unit 
variance,  and  Gaussian  random  variable,  and  Ac(b)  is  the  mean  of  A 
for  a  given  value  of  b  as  given  in  Fig.  11.  We  see  from  Fig.  11  that  the 
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standard  deviation  of  A  may  be  taken  as  5  dB  regardless  of  the  value 
of  b\  the  variations  near  -20  log  (1  -  b)  -  20  are  due  to  small  sample 
problems. 

Figure  12  shows  the  time  during  which  scans  had  fa  in  4  x  1.1-MHz 
frequency  intervals.  It  is,  in  effect,  an  estimate  of  the  density  function 
of  the  distribution  of  f,  and  is,  consequently,  quite  noisy.  The  maxims 
near  ±  30  x  1.1  MHz  from  the  center  of  the  band  are  due  in  part  to  the 
movement  of  estimates  of  f,  to  achieve  realizability.  While,  on  physical 
grounds,  one  would  expect  fa  to  have  a  uniform  distribution,  the  fixed 
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delay  model  is  decidedly  not  a  physical  model.  Consider  a  simulated 
set  of  simple  chree-pach  fades  having  a  uniform  distribution  of  f,.  fixed 
values  for  a  and  b,  and  a  delay  r,  fixed  at  a  value  other  than  6.31  ns. 
This  set  of  fades  will  engender  a  nonflat  probability  density  function 
for  the  fn’s  obtained  in  fitting  to  the  6.31-ns  modeL  The  probability 
density  function  is  flat  within  the  band  regardless  of  the  fixed  delay  of 
the  set  of  simulated  fades;  however,  it  will  more  nearly  resemble  that 
shown  in  Fig.  12  if  the  delay  of  the  set  is  greater  than  6.31  ns  than  if  it 
is  less  than  6.31  ns.  In  short,  Fig.  12  is  characteristic  of  a  channel  with 
a  considerable  fraction  of  delay  differences  greater  than  6  ns. 

Based  on  Fig  12,  we  approximate  the  probability  density  function  of 
f„  by  a  two-level  function.  Note  that  fa  is  defined  on  an  interval  of 
length  1/r,  where  r  is  6.3  ns  the  delay  of  the  modeL  Thus,  with  f„ 
measured  from  the  center  of  the  band,  the  probability  density  function 
for  f0  may  be  approximated  by 


PfSfo) 


or 

T 


T 

3 


_1_ 

4r 


<  l^l  <  2r' 


(15) 


Fig.  II — Mno  tad  standard  deviation  of  cha  distribution  of  -20  log  a  is  s  function 
of  -20  log  (1  -  6>. 
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Fig.  12— Tlnw  that  aodai  pannattr,  fm  wm  in  iaunmia  of  width  4  x  1.1  MHz. 


An  extensive  examination  of  various  conditional  distributions  has 
established  that  than  an  no  other  obvious  and  pervasive  dependencies 
among  the  statistics  of  the  parameters. 

V.  CHANNEL  DELAY  DIFFERENCE 

This  section  preeenta  some  results  obtained  in  estimating  the  chan¬ 
nel  delay  difference.  Some  techniques  described  hen  an  used  in  the 
error  analysis  in  Section  VI.  Three  topics  an  considered  in  this  section. 
First  a  simple  method  is  presented  of  estimating  the  average  delay 
spnad  in  the  channel.  A  second  subsection  shows  that  the  distribution 
of  large  delays  (larger  than  10  ns)  can  be  obtained  for  a  simple  three- 
path  fade  model.  The  delay  distribution  is  shown  to  be  consistent  with 
the  estimate  of  average  delay.  A  third  subsection  illustrates  the  prob¬ 
lem  with  an  observed  channel  shape  that  can  be  matched  most 
successfully  using  a  simple  three-path  model  with  a  delay  of  approxi¬ 
mately  26  ns. 

5. 1  Ms an  dMay  OWIerewce  in  (he  cfiatwi 

The  mean  delay  difference  of  a  channel  that  can  be  characterized  by 
a  simpie  three-path  model  is  easily  estimated.  Consider  a  fade  with  a 
delay,  r.  If  f„  the  frequency  of  the  minimum,  is  uniformly  distributed, 
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the  probability  that  such  a  fade  produces  a  minimum  in  a  band  B  Hz 
wide  is  equal  to  the  ratio  of  the  bandwidth  to  the  spacing  of  the 
minima,  or 


-  Bt.  (16) 

1/r 

If  p(r*) At  is  the  fractional  number  of  fades  having  delays  between  (& 
-  Ddr  and  Adr,  then  the  fractional  number  of  fades  having  a  minimum 
in  band  will  be  P^,  where 

B«ui  Br*p(r*)  dr  ■  Br  (17) 


and 


r  -  £  r*p(r*)dr 


J 


rp(r)dr. 


(18) 


It  follows  from  eq.  (17)  that  one  may  estimate  the  mean  delay,  r, 
from  a  knowledge  of  f*nm,  the  fractional  number  of  scans  having  a 
minimum  in  a  band  of  width  B.  Since  any  method  of  determining  Pm m 
is  acceptable,  consider  estimates  of  Pmm  from  the  parameters  estimated 
using  the  fixed  delay  modeL  The  method  of  estimating  the  frequency 
parameter  in  the  model  involved  moving  null  positions  of  some  fades 
that  had  out-of-band  minima.  These  fades  can  be  excluded  by  using 
only  the  central  two-thirds  of  the  band  in  estimating  r.  Of  the  24,920 
scans  in  the  data  base,  3974  had  minima  between  the  4  th  and  ‘20th 
frequencies.  Hence, 


r 


3974 

24920 


1 

16  x  1.1  x  10* 


■  9.1  ns. 


(19) 


One  might  argue  that  the  mean  delay  should  be  estimated  for  a 
more  carefully  screened  set  of  scans.  Table  I  shows  the  mean  delay 
estimates  obtained  from  scan  populations  qualified  by  having  the 
estimate  of  the  model  parameter  a  in  a  given  5-dB  interval  Table  II 


Table  I — Mean  delay  for  scans  selected  by  value  of 


parameter,  a 


-  20  Lof  a,  dB 

Number  of  Scans 

Scam  with  Min.  in 
Band 

Delay,  f.  ns 

0-5 

101 

31 

17.4 

5-10 

725 

235 

18.4 

10-15 

4299 

875 

11.6 

15-20 

6891 

1161 

9.6 

20-25 

7644 

906 

6.7 

25-50 

4184 

806 

8.2 

30-35 

1019 

159 

8.9 

.All 

24920 

3974 

9.1 
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Table  II— Mean  delay  for  scans  selected  by  value  of 
parameter,  b _ 


-20  Log  1-6.  dB  Number  o {  Seen* 

Scene  With  Min.  in 
Band 

Delay,  f.  na 

0-2 

10.442 

1186 

6.5 

2-4 

7040 

1712 

13.8 

3721 

538 

8.2 

a-8 

1474 

191 

7.4 

9-10 

392 

118 

7.5 

10-12 

527 

68 

7.3 

12- 14 

282 

28 

5.6 

14-16 

190 

21 

6.3 

16-18 

146 

46 

17.9 

18-20 

99 

32 

18.4 

All 

24920 

3974 

9.1 

shows  mean  delay  estimates  qualified  by  the  model  parameter  b ,  which 
specifies  the  shapeliness  of  the  fade. 

With  several  exceptions,  the  estimated  delay  spreads  given  in  Tables 
I  and  II  are  reasonably  constant.  One  exception  is  seen  for  large  values 
of  b  (-20  log  1  -  b  greater  than  16).  This  is  consistent  with  a  channel 
for  which  large  differential  attenuation  across  the  channel  is  mors 
likely  to  occur  when  long  delays  are  present-  The  existence  of  3uch  a 
correlation  should  not  be  surprising.  The  other  exception  is  the  large 
delays  estimated  for  small  values  of  b  and  for  values  of  a  between  0 
and  10  dB.  We  provide  strong  evidence  of  the  existence  of  such  a  class 
of  fades  in  the  next  subsection.  The  existence  of  this  subclass  of  fades 
suggests  that  they  have  a  different  physical  source  than  the  other 
fades  in  the  population. 

5.2  Dlmtrtbution  of  defay  dlftoroncm 

To  further  enhance  our  knowledge  of  the  distribution  of  delay  in  the 
channel,  the  data  base  was  processed  to  extract  a  delay  estimate. 
Recall  that,  for  the  fixed  delay  model,  parameter  estimates  are  chosen 
to  minimise  the  weighted  fit  error  [E  in  eq.  (9)]  for  a  given  fixed  r. 
The  present  calculation  was  performed  for  a  set  of  different  values  of 
r  and  the  value  which  produced  the  smallest  weighted  fit  error  and 
corresponded  to  a  realizable  fade  was  designated  as  the  delay  for  that 
scan  if  it  met  certain  qualifications. 

Because  of  the  degeneracy  in  the  simple  three-path  model,  changing 
the  delay  in  the  fixed  delay  model  will  not  appreciably  improve  the  fit 
for  any  3can  that  can  be  well  approximated  by  a  fixed  delay  of  6  ns  or 
less.'  In  performing  the  optimum  delay  calculation,  the  weighted  fit 
error  was  minimized  for  a  predetermined  set  of  delays;  the  differences 
between  adjacent  delay  values  were  chosen  to  be  approximately  15 
percent.  A  given  scan  was  assigned  a  delay  different  from  6.3  ns  only 
if  the  third  best  value  of  the  weighted  fit  error  was  at  least  0.1  dB 
worse  than  the  best  value.  (We  use  the  third  best  value  because  we 
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must  examine  three  values  to  detect  a  minimum.)  This  criterion  sets 
a  threshold  on  the  acceptable  sharpness  of  the  minimum  in  the  fit 
error  with  respect  to  changes  in  delay. 

The  selection  criterion  was  chosen,  after  several  iterations,  to  insure 
regularity  in  the  estimates  derived  from  successive  scans.  With  the 
chosen  criterion,  the  scans  that  were  assigned  a  new  delay  occurred  in 
groups  of  consecutive  scans  and  may  be  said  to  constitute  fading 
events.  During  any  of  these  events,  the  delay  was  consistent  in  that 
indicated  delays  were  within  ±15  percent.  If  we  assume  that  the 
physical  channel  does  not  change  between  scans,  we  can  associate  a 
time  with  each  scan  and  plot  the  distribution  of  the  time  periods 
during  which  the  characterizing  delay  was  greater  than  a  specified 
delay. 

A  series  of  such  plots,  conditioned  on  the  concurrently  estimated 
value  of  b,  is  given  in  Fig.  13.  The  uppermost  curve  contains  the  data 
derived  from  all  scans  which  met  the  selection  criterion;  its  shape  is 


OSLAY  in  NANOSECONDS 

Fig.  13— Distribution  of  optimum  delay  for  simple  three-path  model  ax  qualified 
by  realizability,  the  sharpness  of  minimum,  and  by  several  value*  of  the  model 
parameter  b. 
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dominated  bv  the  627  seconds  during  which  the  channel  was  best 
modeled  by  a  delay  of  43  ns  (the  largest  delay  in  the  test  set)  or  more, 
but  had  littie  shape  ( b  <  0.115).  These  characteristics  contribute  to 
the  large  (18  ns)  mean  u  ’ays  noted  in  the  previous  subsection  for 
small  values  of  b.  They  may  be  due  to  quantization  but  are  apparently 
not  artifacts  of  the  estimation  scheme.  Although  the  origin  of  this  type 
of  channel  defect  is  currently  not  understood,  it  should  not  trouble 
any  existing  radio  system. 

It  is  apparent  from  the  distributions  in  Fig.  13  that  very  few  scans 
qualified  for  a  new  delay  with  delays  less  than  10  ns.  Consequently, 
the  distribution  should  not  be  trusted  for  delays  less  than  12  or  15  ns; 
beyond  15  ns.  it  may  be  interpreted  as  a  lower  bound  to  the  true 
distribution.  The  three  curves  qualified  by  the  parameter  b  correspond 
to  fades  with  peak- to -peak  variability  of  2,  3,  and  6  dB.  (Peak-co-peak 
variability  is  20  log(l  +  b/l  -  b),  as  may  be  seen  in  Fig.  1.)  If  the  delay 
were  exponentially  distributed,  the  distribution  of  delay  would  be  a 
straight  line  on  Fig.  13  and  would  have  the  form 

Pi r  >  x)  -  e',/7.  (20) 

Fitting  a  straight  line  to  the  three  distributions  in  Fig.  13  for  which  b 
>0.115  shows  that  the  average  delay  decreases  with  increasing  b.  The 
corresponding  values  are  5,  5.5,  and  11  ns.  Note  that  this  implies  that 
b  and  r  in  a  simple  three-path  model  are  not  independent. 

5.3  An  axampio  of  a  long  daisy  scan 

To  confirm  the  existence  of  long  delay  scans,  consider  an  event  that 
covered  approximately  10  seconds  on  22  June  1977,  from  23  h,  28  m,  54 
s.  A  representative  scan  from  the  middle  of  this  period  is  shown  with 
the  fit  obtained  with  the  fixed  delay  (6.3  ns)  model  in  Fig.  14a.  To 
emphasize  the  consistency  of  this  channel  condition,  an  average  of  the 
channel  condition  for  the  central  4.2  seconds  (21  scans)  of  this  event 
is  compared  to  the  selected  scan  in  Fig.  14b. 

It  is  apparent  from  Fig.  14a  that  the  6.3  ns  delay  does  not  have 
enough  curvature  (delay  is  too  short)  to  precisely  model  the  channel 
shape.  Figure  15  shows  the  same  scan  modeled  by  three-path  fades 
having  delays  of  22.7,  26,  and  30.3  ns.  The  26-ns  fit  is  the  best;  it  has 
a  weighted  fit  error  0.4  dB  better  than  the  22.7-ns  fit  and  0.8  dB  better 
than  the  30.3  ns  fit.  However,  the  closeness  of  all  three  fits  illustrates 
the  difficulties  in  estimating  channel  delay  differences.  Visually,  one 
would  choose  the  26-ns  model  on  the  basis  that  the  30.3-ns  fit  has  too 
much  curvature  and  the  22.7-ns  fit  too  little. 

VI.  ERROR  ANALYSIS 

To  verify  that  the  model  adequately  represents  the  transmission 
characteristics  of  the  channel,  we  examine  the  errors  between  the 
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Fig.  14 — Scan  from  22  June  1977,  23  h.  28  m.  48.6  s.  (a)  Comparison  with  fixed  delay 
model,  (b)  Comparison  with  average  of  scans  from  23  h.  28  m.  46.4  s  to  23  h.  28  m. 
50.4  s. 


channel  as  observed  and  as  modeled.  In  this  section  we  consider  the 
statistics  of  the  rms  errors  and  the  maximum  errors. 

6.1  RMS  errors 

A  useful  measure  of  the  quality  of  the  fit  of  the  model  to  a  given 
channel  scan  is  the  root-mean-square  value  of  the  decibel  error  at  each 
of  the  sampled  frequencies.  Denoting  this  error  as  Em.,  we  have 

1/2 

(21) 

n#19 

The  model  parameters  were  estimated,  as  described  in  Section  III,  to 
minimize  the  error,  E ,  which  is  a  weighted  sum  of  the  squares  of  the 
power  differences  at  each  frequency  [see  eq.  (9)].  The  weighting  was 
chosen  [eq.  (10)]  30  that  the  error  E  would  approximate  the  error  Em. 
as  given  by  eq.  (21).*  Indeed,  one  may  show  directly  that  the  two 
expressions  are  equivalent  as  long  as 

T 

1 1  -  | «  1  for  all  n.  (22) 


2  (dB  error  at  fn)* 

40  a- 1 


*  Note  that  the  parameter  estimation  problem  cannot  be  solved  in  closed  form  by 
minimizing  Em. 
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As  we  have  seen  in  Fig.  14,  this  inequality  is  not  always  satisfied. 
Consequently,  in  using  Em u  as  a  standard  of  comparison,  we  are 
evaluating  not  only  how  well  the  model  fits  the  observed  channel,  but 
also  how  well  we  have  chosen  the  parameters  to  make  the  match. 

The  error  Erm*  is  a  desirable  quantity  to  work  with  because  we  can 
estimate  its  distribution  under  the  assumption  of  perfect  matching. 
We  observe  that  if  the  decibel  error  were  Gaussian  with  unit  variance 
and  zero  mean,  23  Elm,  would  be  a  x'  variable  with  20  degrees  of 
freedom  (to  account  for  the  three  parameters  estimated  per  scan). 
Observations  of  a  simulated  channel  with  the  transmitter  and  receiver 
back-to-back  indicate  that  the  instrumentation  errors  are  approxi¬ 
mately  Gaussian  with  a  standard  deviation,  Oi,  of  about  0.65  dB. 
Observations  of  the  channel  at  mid-day  with  the  channel  nominally 
flat  and  unfaded  indicate  that  the  standard  deviation  of  the  errors  is 
between  0.68  and  0.73,  varying  frequency  to  frequency  and  day  to  day 
by  a  few  hundredths  of  a  decibel.  Hence,  if  we  enter  a  table  of  the  x' 
distribution,  Q(xJ  1 20),  with 


■> 

X' 


23  EL. 

i 

<J7 


(23) 
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'vo  con  determine  the  distribution  of  Em u  under  the  assumption  of 
perfect  matching.*  This  distribution  is  shown  as  a  reference  on 
Figs.  16  and  17.  It  is  indicated  by  a  solid  curve  labeled  “ideal”  for  a, 
«  0.70  and  by  o’s  for  <r,  »  0.75. 

Figure  16  presents  the  distribution  of  the  rms  error  for  two  scan 
subpopulations  using  the  fixed  delay  (6.3  ns)  model  The  subpopulation 
of  the  distribution  Labeled  “standard”  consists  of  all  scans  that  could 
be  modeled  directly;  the  distribution  labeled  “modified”  shows  the  rms 
error  distribution  for  all  scans  which  required  an  adjustment  of  the 
frequency  of  the  modeled  fade  to  achieve  realizability.  Figure  17  shows 
the  distribution  of  the  rms  error  for  the  composite  of  all  samples  using 
the  fixed  delay  (6.3  ns)  modeL  The  distribution  labeled  simple  three- 
path  model  indicates  the  error  distribution  that  was  obtained  when 
the  scan  fitting  allowed  unqualified  variation  in  model  delay  to  achieve 
the  best  fit.  That  is,  the  calculation  described  in  Section  5.2  was 
performed  and  the  results  were  qualified  only  on  the  basis  of  realiza- 
bility.f 

In  each  case  described  above,  the  mean  value  of  the  rms  error  is 
close  to  the  median  value.  For  the  two  subpopulations  shown  in  Fig. 
16,  the  calculated  mean  fit  errors  correspond  to  <r,  values  of  0.76  and 
0.85  dB,  or  the  errors  are  about  0.09  dB  larger  when  a  realizable  fit  is 
obtained  by  varying  the  frequency  of  the  model  minimum.  Comparing 
the  composite  distributions  in  Fig.  17,  we  find  that  the  mean  error  in 
the  fixed  delay  (6.3  ns)  model  corresponds  to  <j,  —  0.78  dB  or  about 
0.08  dB  higher  than  that  observed  when  the  channel  is  quiescent.  The 
simple  three-path  model  has  a  distribution  of  rms  error  that  very 
nearly  matches  the  ideal  distribution  (with  19  degrees  of  freedom)  for 
a,  -  0.75.  This  is  consistent  with  the  instrumentation  error  imputed  to 
the  standard  distribution  in  Fig.  16  and  is  indicative  of  the  instrumen¬ 
tation  error  in  the  presence  of  multipath  fading.  It  is  exceptionally 
good  considering  that  the  data  are  obtained  from  time  sequential 
measurements  on  a  dynamically  changing  channeL  One  concludes  that 
the  modeling  error  is  negligible  for  the  simple  three-path  model  For 
the  fixed  delay  model  under  the  assumption  that  the  instrumentation 
and  modeling  errors  add  in  quadrature,  the  modeling  error  has  a 
tolerable  value  on  the  order  of  0.2  dB.  That  is, 

[(0.75)1  +  (0.2)2]1/*  -  0.776. 

The  tails  of  the  distributions  in  Figs.  16  and  17  for  large  errors  are 
of  considerable  interest.  The  tails  near  small  values  are  of  little 

*  From  the  central  limit  theorem.  we  know  that  Elm  will  be  approximately  Gaussian, 
as  la  y1,  regardless  of  whether  or  not  the  measurement  errora  are  precisely  Gauaaian. 

t  Note  that  although  one  cannot  alwaya  reliably  localise  the  values  of  the  parameters 
in  fitting  with  the  simple  three-path  model  (see  discussions  in  Section  2.1  and  Ref.  7). 
the  error  in  the  fit  is  always  well  defined. 
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Pig.  IS— Distribution  of  rmt  fit  error  for  two  scan  subpopuistions  with  fixed  daisy 
16.3  ns)  model 

consequence;  they  are  distorted  by  quantization  because  one  cannot 
associate  any  error  with  the  12  flat  fades  included  in  the  data  base. 
The  deviation  of  the  distributions  from  the  ideal  distribution  at  large 
errors  is  significant. 

The  large  deviation  of  the  modified  fits  in  Fig.  16  reflects  the  failure 
of  the  fixed  delay  (6.3  ns)  model  to  accurately  fit  the  long  delay  fades. 
The  tail  deviation  from  ideal  is  modest  down  to  about  the  0.5  percent 
level,  corresponding  to  a  few  tens  of  seconds  per  month.  For  compar- 
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ison,  we  note  that  the  rms  error  of  the  fit  shown  in  Fig.  14a  is  2.3  dB; 
this  was  the  worst  fit  encountered  for  the  fixed  delay  (6.3  ns)  model 
However,  even  in  this  case  the  model  failure  is  hardly  describable  as 
severe.  The  model  of  the  channel  is  depressed  by  40  dB  and  has  9.5  dB 
of  gain  slope;  the  actual  channel  is  depressed  by  39  dB  and  has  11  dB 
of  gain  slope.  Also,  we  note  that  the  6.3- ns  delay  model  has  the 
response  minimum  at  about  the  same  frequency  as  the  best  represen* 
tation,  the  26- ns  delay  model  shown  ip  Fig.  15. 

NEW  SELECTIVE  FADING  MODEL 
6-28 


The  deviation  of  the  tail  of  the  error  distribution  for  the  three- path 
fade  (Fig.  17)  reflects  the  fact  that  there  are  fades  that  even  this  model 
has  difficulty  in  fitting.  An  example  of  such  a  fade  is  shown  in  Fig.  18 
along  with  the  fit  provided  by  the  fixed  delay  (6.3  ns)  model  The  same 
rms  error  (1.6  dB)  is  obtained  for  ail  values  of  model  delay  between 
0.05  and  9  ns;  the  fit  degrades  for  larger  delays.  Either  more  than  three 
rays  are  needed  to  describe  the  channel  shape  in  Fig.  18,  or  the  channel 
is  so  depressed  that  the  amplitudes  in  the  notch  are  distorted  due  to 
closeness  to  the  noise  level  in  the  measuring  equipment.  The  scan 
shown  in  Fig.  18  is  one  of  three  similar  scans  and  has  little  statistical 
significance. 

6.2  Maximum  errors 

Another  type  of  error  that  can  be  used  to  judge  the  quality  of  the  fit 
of  the  model  to  the  channel  is  the  worst-case  error.  That  is,  after  fitting 
to  each  scan,  one  records  the  magnitude  of  the  largest  difference  (in 
decibels)  between  the  observed  channel  shape  and  the  shape  calculated 
from  the  modeL  The  following  paragraphs  consider  the  distribution  of 
these  worst-case  errors. 

As  in  the  preceding  subsection,  we  can  calculate  an  ideal  distribution; 
however,  the  ideal  distribution  is  not  as  realistic  in  this  case  since  it  is 
strongly  dependent  on  the  tails  of  the  distributions  of  the  individual 
measurement  errors.  We  assume  that  each  power  measurement  had 


an  error  in  decibels  that  was  Gaussian,  with  zero  mean,  a  standard 
deviation  of  v  20/23  <r,  to  account  for  the  three  parameters  estimated 
from  the  23  observations  per  scan,  *  and  that  the  errors  are  independent 
frequency  to  frequency.  If  the  probability  of  any  one  measurement 
having  an  error  less  than  x  is  denoted  by  Pi(x),  the  probability  that  all 
23  have  values  less  than  x  is 

PM(x)  -  [Pi(x)]“.  (24) 

This  is  the  probability  that  the  maximum  error  is  less  than  x,  whereas 
we  want  the  probability  that  it  is  greater  than  x  which  we  denote 
Qb(x).  It  follows  immediately  from  eq.  (24)  that 

Qm<x)  -  1  -  [P^x)]43 

-  1  -  [1  -  Qt(x)]a.  (25) 

The  distribution  given  by  (25)  is  used  as  a  reference  in  F:~s.  19  and 
20,  which  show  the  distribution  of  the  maximum  error  L  re  same 
cases  as  in  Figs.  16  and  17.  Since  the  tails  of  these  distributions  are 
well  behaved  for  larger  errors,  the  distribution  of  the  maximum  errors 
is  apparently  dominated  by  the  instrumentation  noise.  That  is,  if  we 
use  for  the  standard  deviation  of  the  measurement  noise  the  value 
obtained  from  the  mean  of  Erm.  for  one  of  these  cases  (as  given  in 
Section  6.1),  the  resulting  worst-case  error  distribution  calculated  with 
eq.  (25)  will  closely  match  the  observed  maximum  error  distribution. 

Vtl.  CONCLUSIONS 

By  analyzing  the  errors  in  fitting  the  observed  channel  characteris¬ 
tics  in  Section  VI,  we  demonstrated  that  the  simple  three-path  fade 
model  is  indistinguishable  from  a  perfect  model  of  a  line-of-sight 
microwave  radio  channel. 

The  simple  three-path  model  was  used  in  Section  V  to  characterize 
the  channel  delay  difference.  By  two  different  methods,  it  was  shown 
that,  when  there  is  3  dB  or  more  shape  present  in  the  channel,  the 
average  delay  difference  is  between  5  and  8  ns.  We  developed  a  lower 
bound  on  the  tails  of  the  distribution  of  delay  difference.  From  these 
results,  which  are  shown  in  Fig.  13,  we  observe  that  a  differential 
channel  attenuation  in-band  of  3  dB  or  more  may  be  due  to  delay 
differences  as  great  as  43  ns.  In  another  dimension,  one  would  expect 
to  see  differential  attenuation  of  3  dB  or  more  in-band  due  to  delays 
greater  than  20  ns  for  at  least  70  seconds  in  a  heavy  fading  month. 
This  is  comparable  to  the  time  the  channel  attenuation  at  a  single 
frequency  exceeds  40  dB. 

*  For  comparison*  with  tht  thrss-psth  modti.  it  is  appropritM  to  us*  v  19/23 a,. 
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Fig.  19 — Distribution  of  mianiman  idfi)  fit  anor  for  two  scan  subpopuianona  with 
filed  daisy  (6.3  ns)  modal. 


From  the  error  analysis  in  Section  VI,  we  also  conclude  that  the 
fixed  delay  (6.3  ns)  modei  is  a  very  good  approximation  to  the  channel 
for  all  observed  conditions.  This  conclusion  is  further  substantiated  by 
Figs.  14  and  18,  which  show  the  scans  for  which  the  fits  with  the  fixed 
delay  model  exhibited  the  largest  rms  fit  error  (2.3  dB)  and  the  largest 
maximum  error  (3.9  dB),  respectively.  The  fixed  delay  model  is  pref¬ 
erable  to  the  three-path  modei  for  channel  modeling  because  it  requires 
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Fig.  20— Distribution  of  maximum  idB)  fit  error  for  compoaita  population  with  fixed 
delay  and  simple  three-path  models. 


only  three  parameters,  and  these  can  always  be  uniquely  determined 
from  a  channel  amplitude  scan. 

The  statistics  of  the  parameters  of  the  fixed  delay  model  as  described 
in  Section  IV  and  shown  in  Figs.  S  to  12  provide  the  means  of 
statistically  generating  all  the  channel  conditions  that  one  expects  to 
see  on  a  nominal  hop  channel  operated  at  6  GHz.  If  one  determines, 
by  laboratory  test,  the  parameter  values  that  will  cause  a  particular 
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error  rate  in  a  digital  radio  system,  one  can  easily  calculate  the  time 
during  a  heavy  fading  month  that  the  error  rate  will  equal  or  exceed 
this  critical  value.  A  companion  paper  describes  the  laboratory  test 
and  the  required  calculations.* 

Future  work  will  be  directed  toward  verifying  the  model  and  model 
statistics  with  additional  fading  data  obtained  both  at  6  GHz  and  at  4 
GHz.  Using  coherent  data  obtained  in  1973.  it  will  be  possible  to 
determine  the  extent  to  which  the  channel  is  actually  a  minimum 
phase  channel. 
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APPENDIX 

Estimation  of  Parameters 

The  problem  of  estimating  the  parameters  a,  0,  and  fa  in  Section  m 
is  equivalent  to  the  problem  of  determining  the  first  three  terms  in  a 
subharmonic  Fourier  series  expansion  of  a  function  in  the  frequency 
domain.  Since  such  expansions  are  not  standard,  we  provide  a  complete 
description  of  the  methodology  here. 

From  eqs.  (4)  and  (9),  we  may  express  the  weighted  mean-square 
error  between  estimated  and  observed  power  as* 

B_SC,(y,-«  +  ^  coslw,  -  uo)r)2 

E - fa. -  (26) 

For  simplicity,  we  use  a  normalized  weighting  function,  defined  by 


so  that 

l  dn  -  1.  (28) 

In  terms  of  the  normalized  weighting  we  may  write  (26)  as 

E  -  £  y,  -  a  +  0  coslw,,  -  wo)r)J  (29) 

*  Throughout  this  appendix,  all  automations  are  taken  over  ail  valuta  of  n  correspond- 
in*  to  ail  frequencies  observed  in  a  scan. 
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or  in  expanded  form  as 

E  -  £  d„ri  +  a2  +  t?  £  AcosJ(w,  -  wo)r 

+  2/3  £  dny/icosiw,  -  uo )r 

-  2a0  y  d*cos(wn  -  <jo)t  -  2«  £  d,y,  .  (30) 

The  error  £  is  a  minimum  when  a,  0,  and  wo  are  chosen  so  chat  the 
partial  derivatives  of  E  with  respect  to  a,  0,  and  wo  are  all  equal  to 
zero.  Setting  the  partial  derivative  of  eq.  (30)  with  respect  to  0  equal 
to  zero  and  solving  for  0  gives 

q  £  dwcosiw,  -  wq)t  —  £  d„Yn coa(w„  -  wq)t 

£  ckcos^w.  -  wo)t  ''*J 

Substituting  (31)  into  eq.  (30),  we  find  Eg,  the  error  minimized  with 
respect  to  0,  as 

■  {a2[£  d«cos2(w«  —  wo)T  —  (£  <f,cos(w«  —  wo)r)-“] 

-  2q  [(£  d»y»)(£  of»cos2(wn  -  w®)t) 

-  (£  ct,coa(w„  -  wo )r)(£  d,]f«coe(w,  -  wo)t)] 

-  (Xd,y»cos(w,-wo)T)2).  (32) 

Minimizing  this  with  respect  to  a  requires  that  we  set  the  partial 
derivative  of  E$  with  respect  to  a  equal  to  zero.  This  gives 

(£  ckT-iXS  d,cos2(w,  -  wo)t)  - 

a  -  (£  <j,cos(w,  ~  wq)t)(£  d»YK cos(w„  -  wo)t) 

£  dKc os^w,  -  wo)t  -  (£  d,cos(w,  —  wo)r)2 

Substituting  (33)  into  (32)  gives  £«*  the  error  minimized  with  respect 
to  both  a  and  8,  as 

£*-£d,n-?1 

(£  dn(Yo  -  ?)coa(w,  -  wo)r)2  (34) 

V  dn cos3(w,  -  wo)r  -  (£  d„coa(w„  -  wo)r)J  ’ 


where 


?-£d,y. 


We  note  that  we  could  have  obtained  this  same  expression  by  first 
minimizing  with  respect  to  a  and  then  with  respect  to  0;  however,  one 
obtains  different  but  equivalent  expressions  for  a  and  0,  depending  on 
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the  order  of  differentiation.  We  develop  the  alternative  expressions  for 
a  and  0  in  the  following  paragraphs. 

Let  us  define  some  new  quantities  to  simplify  these  expressions.  Let 
the  difference  between  the  observed  power  and  the  weighted  mean 
power  in  the  band  be  designated  by  X„;  then 

X„  -  Y.  -  l  d«Y*  -  -  ?.  (36) 

If  we  also  define  the  quantities 

Xe  m  £  d„X,coa  Wnf,  (37) 

X.  *  2  aLXain  w»r,  (38) 

Da  “  £  dnCOtfiUn  —  ta)o)T,  (39) 

Db  m  £  ducoelwp,  —  wo)t,  (40) 

we  may  rewrite  a  and  0  from  eqs.  (31)  and  (33)  as 

[Xcoa  wot  +  Xmd  uor]Db 

•m? - dT^dI -  1411 


0  -  jr  {(a  -  7)Db  -  (.Xrcos  wot  +  Xsin  wot)}  . 

Dm 


Using  (41)  to  eliminate  a  from  (42),  we  obtain 
.  Xcoe  ujbt -*■  X.sin  uot 

* - ~d7-~dI - ■  (43) 

We  may  use  (43)  in  (41)  to  obtain 

a  ■  7  +  0Db  ■  144) 

Equationa  (43)  and  (44)  are  the  estimators  that  would  have  been 
obtained  if  the  order  of  taking  partial  derivatives  in  the  preceding 
development  had  been  reversed.  It  is  apparent  that,  after  one  has 
estimated  wo,  one  may  estimate  a  and  0  by  using  either  eqs.  (41)  and 
(42),  (43)  and  (44).  or  eqs.  (41)  and  (43). 

The  estimate  of  wo  that  minimizes  the  weighted  error  is  obtained  by 
minimizing  £«.  with  respect  to  wo.  Using  eqs.  (35)  to  (40)  in  eq.  (34). 
we  write 

E..  -  £  d.Xl  -  ^  -y.t  145) 

To  see  the  explicit  dependence  of  £*  on  wo,  we  define  the  following 
quantities 

<&-£<£ cosJw„r  -  (£  d,cos  waT)1.  (46) 
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d,  ■  £  <f„sin*M„r  -  ( £  eksin  <j«r)2, 
d„  ■  £  ckcos  w„r  3in  u)„r  —  (£  cf„cos  aj,,r)(£  dn sin  oj„t). 
Substituting  these  into  (45)  gives 

[Xcos  u„r  +  .Y.sin  uof]2 


£*-£  - 


dccoa~  wor  +■  2dctcos  u*yr  sin  uor  +  casin' uof  ' 


(47) 

(45) 

(49) 


Setting  the  partial  derivative  of  £&,,  as  given  by  (49),  equal  to  zero 
gives  the  estimator  for  u>o  as 


(JOT 


Tan- 


dcX.  -  dc.X< 
d.Xc  -  dr.X. 


(50) 


Obviously,  two  values  of  u *>r  in  the  interval  (—■ r,  t]  will  satisfy 
eq.  (50).  One  of  these,  the  principal  value,  lies  in  the  interval  ( — w/2, 
rr/2],  the  other  differs  from  the  first  by  ±ir.  We  shall  show  that  the  two 
solutions  are  equivalent,  but  that  our  chosen  solution  is  unique. 

If  we  replace  uwr  by  u*>r  ±  it  in  eqs.  (39),  (40),  (43),  (44),  and  (45),  we 
see  that  £>»  and  /?  change  sign  and  a  and  E&,  are  unchanged.  Since  we 
want  the  solution  with  0  greater  than  zero,  we  take  the  principal  value 
solution  to  (50)  if  the  resulting  estimate  of  0  is  positive.  Otherwise  we 
add  or  subtract  it  to  obtain  a  positive  value  for  /3  and  a  value  of  u*>t  in 
the  appropriate  interval. 

While  we  could  substitute  the  result  of  eq.  (50)  into  (49)  to  obtain 
the  minimum  error,  £.&«,,  it  is  more  generally  useful  to  evaluate  E„a 
for  the  optimum  uo-  This  is  especially  true  when  we  do  not  use  the 
optimum  u*>,  as  given  by  eq.  (50).  The  simplest  form  for  Ea*  is  obtained 
by  substituting  (43)  into  (45)  to  give 


E» i  -  £  d„Xi  -  [Da  -  Dl)/?.  (51) 


These  equations  were  implemented,  with  the  modifications  de¬ 
scribed  in  Section  3.2,  to  obtain  all  the  fits  described  in  this  paper. 
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Digital  Radio  Outage  Due  to  Selective 
Fading — Observation  vs  Prediction  From 
Laboratory  Simulation 

By  C.  W.  LUNDGREN  and  W.  D.  RUMMLER 
(Manuscript  received  July  7,  1 978) 

A  statistical  model  (introduced  in  a  companion  paper)  of  fading  on 
a  radio  path  is  used  with  laboratory  measurements  on  a  digital  radio 
system  to  estimate  the  outage  due  to  multipath  fading,  where  outage 
is  the  time  that  the  bit  error  rate  (BBR)  exceeds  a  threshold.  Over  the 
range  of  BEX  of  interest  ( 10T*  to  l(T3),  the  calculated  outage  agrees 
favorably  with  the  outage  observed  during  the  period  for  which  the 
fading  model  was  developed.  It  is  further  shown  that  the  calculated 
outage,  when  scaled  to  a  heavy  fading  month  on  the  basis  of  single¬ 
frequency,  time-faded  statistics,  agrees  equally  well  with  the  outage 
observed  on  the  same  path  during  a  heavy  fading  month.  The  agree¬ 
ment  between  measured  and  predicted  outage  substantiates  the  se¬ 
lective  fading  model.  The  prescribed  laboratory  measurements  char¬ 
acterize  the  sensitivity  of  the  radio  system  to  selective  fading.  Thus, 
the  methodology  provides  a  useful  basis  for  comparing  the  outage  of 
alternative  realizations  of  digital  radio  systems. 


I.  INTRODUCTION 

Present  interest  in  using  high-speed  common  carrier  digital  radio1-5 
has  precipitated  a  need  for  estimating  the  performance  of  such  systems 
during  periods  of  selective  (multipath)  fading.  This  paper  describes  a 
method  of  characterizing  a  digital  radio  system  in  the  laboratory  which 
allows  the  outage  to  be  accurately  estimated.  For  a  digital  radio 
system,  outage  requirements  are  stated  in  terms  of  the  number  of 
seconds  in  a  time  period  (usually  a  heavy  fading  month)  during  which 
the  bit  error  rate  (ber)  may  exceed  a  specified'  level;  typically,  10-3  or 
10~*  is  appropriate  to  voice  circuit  application. 

The  method  is  based  upon  a  statistical  channel  model5  developed 
from  measurements  on  an  unprotected  26.4-mile  hop  in  the  6-GHz 
band  in  Palmetto,  Georgia  in  1977  using  a  general  trade  8-psk  digital 
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radio  system  as  a  channel  measuring  probe.  The  modeled  fading 
occurrences  were  scaled  to  the  basis  of  a  heavy  fading  month  using  the 
occurrence  of  time  faded  below  a  level  at  a  single  frequency  as  the 
means  of  calibration.  The  bit  error  rate  performance  of  the  digital 
radio  system  was  measured  during  the  time  period  used  for  channel 
modeling  and  for  an  extended  period  corresponding  to  a  heavy  fading 
month.  This  same  radio  system  was  later  subjected  to  a  measurement 
program  in  the  laboratory  using  a  multipath  simulator  which  provides 
a  circuit  realization  of  the  fading  model.  The  measured  results  are  used 
with  the  channel  model  to  determine  the  occurrence  of  channel  con¬ 
ditions  which  will  cause  the  ber  to  exceed  a  given  threshold.  Compar¬ 
isons  on  the  basis  of  the  modeling  period  and  a  heavy  fading  month 
show  good  agreement  between  calculated  and  observed  outages  for 
bebs  between  10“*  and  10" 3 . 

The  properties  of  the  fixed-delay  channel  model  are  reviewed  briefly 
in  Section  Q  as  a  basis  for  describing  the  measurements  and  for  the 
subsequent  outage  calculations.  This  three-parameter  channel  model 
is  realized  in  the  laboratory  by  an  17  fade  simulator.  The  simulator 
and  its  use  in  obtaining  the  necessary  laboratory  data  are  described  in 
Section  m.  The  procedures  to  be  followed  in  calculating  outage  times 
for  a  given  bek  are  described  in  Section  IV.  Calculated  and  observed 
outage  times  are  compared  in  Section  V.  Conclusions  are  provided  in 
Section  VL 

II.  MODEL  DESCRIPTION— METHODOLOGY 

It  has  been  demonstrated6  that  the  complex  voltage  transfer  function 
of  a  line-of-sight  microwave  radio  channel  is  well  modeled  by  the 
function 


ff{w)  -  a  [1  -  ie-'1--'*1"] 


(1) 


with  r  fixed.  A  S-GHz  channel  (30- MHz  bandwidth)  has  been  charac¬ 
terized  statistically  by  the  model  with  r  =»  6.3  ns.  Such  a  channel  has 
a  power  transfer  function  given  by 

|  H( w)  | 1  —  a2[l  +  b1  —  2b  cosfw  —  wo)r]  (2) 

and  an  envelope  delay  distortion  function,  i.e.,  the  derivative  of  the 
phase  of  H(u)  with  respect  to  u,  given  by 


D(w) 


&t(cos(u)  -  (jq)t  —  b) 

1  +  b1  -  2b  cos (u  —  wo)t' 


(3) 


In  the  following  paragraphs,  we  summarize  the  properties  of  the  model 
the  statistics  of  the  model  parameters,  and  the  measurement  objec¬ 
tives. 


2. 1  Fixed  delay  model 

A  plot  of  the  attenuation  produced  by  the  fixed  delay  model  of  eq. 
(1)  is  shown  in  Fig.  1.  Since  the  delay  is  fixed  at  6.3  ns,  che  spacing 
between  nulls  of  the  response,  156.4  MHz,  is  much  larger  chan  the 
channel  bandwidth.  The  parameters  a  and  b  control  the  depth  and 
shape  of  the  simulated  fade,  respectively.  The  parameter  f0(~wo/2x) 
determines  che  position  of  the  fade  minimum  or  notch.  Both  the  notch 
frequency,  f0 ,  and  the  response  frequency,  f,  are  measured  from  the 
center  of  the  30-MHz  channel  for  convenience. 

The  model  function  of  eq.  (1)  may  be  interpreted  as  the  response  of 
a  channel  which  provides  a  direct  transmission  path  with  amplitude  a 
and  a  second  path  providing  a  relative  amplitude  b  at  a  delay  of  6.3  ns 
and  with  a  phase  of  ojot  +  rr  (independently  controllable)  at  the  center 
frequency  of  the  channel.  This  interpretation  is  represented  in  Fig.  2 
by  a  phasor  diagram  at  w  —  0,  the  center  frequency  of  the  channel. 
Varying  the  frequency,  u,  over  the  channel  bandwidth  (30  MHz)  moves 
the  angle  of  the  interfering  ray  through  an  arc  of  about  60  degrees  (2tt 
x  30  MHz  x  6.3  ns  ~  ir/3),  centered  at  the  position  shown.  This 
diagram  is  useful  for  understanding  the  fade  simulation;  it  also  provides 
an  alternate  means  of  describing  the  position  of  the  notch.  The  notch 
position  may  be  specified  by  its  frequency,  f„  or  by  <*>,  the  angle  of  the 
interfering  path  at  the  center  of  the  channel. 


l-f,,  MHl 

Fig.  1— Attenuation  of  channel  model  function.  Hiu\  •  a(l  -  b  exp <-/(«  -  wnlri],  for 
~  m  6.3  ns,  a  "  0.1,  b  «  0.7. 
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From  Figs.  1  and  2,  or  eq.  (1),  it  may  be  seen  that  varying  a  changes 
the  overall  level  and  varying  b  changes  the  shapeliness  of  the  modeled 
fade.  Furthermore,  if  the  minimum  is  within  the  30-MHz  channel 
bandwidth  ( |$|  <  30°),  the  fixed  delay  model  can  generate  notches 
with  a  wide  range  of  levels  and  notch  widths.  With  the  minimum  out 
of  band,  it  can  generate  a  wide  range  of  combinations  of  levels,  slopes, 
and  curvatures  within  the  channel  bandwidth. 

2.2  Mod*  statistics 

The  statistics  of  the  model  parameters  were  obtained  from  a  selected 
data  base  during  which  heavy  fading  activity  was  observed.0  The 
distribution  of  b  is  best  described  in  terms  of  B  —  —20  log(l  —  b). 
Figure  3  shows  the  distribution  of  B  and  the  least-squares  straight  line 
fit  to  the  distribution  over  the  region  where  it  best  represents  selective 
fading — between  B  values  of  3  and  23  dB.  The  channel  is  described  by 
B  greater  than  23  dB  for  less  than  0.15  percent  of  the  observed  time 
which  makes  the  distribution  less  certain  beyond  this  point.  At  the 
other  extreme,  during  the  periods  of  time  when  there  is  little  or  no 
selective  fading,  the  channel  is  characterized  by  values  of  B  less  than 
3  dB.  Thus,  the  fitted  line  represents  a  lower  bound  on  the  distribution 
for  B  less  than  3  dB.  Since  the  fitted  line  has  an  intercept  of  5400 
seconds,  we  may  model  the  fraction  of  5400  seconds  during  which  B 
exceeded  a  value  X  by  the  probability  distribution 

P{B>X)~  e-xn».  (4) 

Thus  the  probability  of  finding  a  value  of  B  between  X  and  X  +  dX  is 

pB(X)  e"™»  (5) 

0.0 

The  distribution  of  a  is  lognormal  with  a  standard  deviation  of  5  dB 
and  a  mean  that  is  dependent  on  B  (or  6).  Hence,  the  probability  that 
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Fig.  3 —  Distribution  of  B  for  model  data  baw  period. 

A  •  -20  log  a  has  a  value  between  Y  and  Y  +  dY  is  given  by 

pA(  Y)  dY  -  (6) 

The  relationship  between  Ao,  the  mean  of  the  distribution,  and  B  is 
given  in  Fig.  4. 

The  distribution  of  f0  is  found  to  be  independent  of  A  and  B.  It  is 
usually  simpler  to  work  with  <t>  rather  than  f0.  The  two  variables  are 
simply  related  in  that  <t>  is  defined  on  the  interval  (— ir,  it)  and  a  2.5* 
degree  change  in  <t>  corresponds  to  a  1.1-MHz  change  in  f*.  For  the 
fixed  delay  model,  the  variable  $  has  been  found  to  have  a  probability 
density  function  that  can  be  described  as  uniform  at  two  levels,  with 
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Fig.  4 — Mean  and  standard  deviation  of  the  distribution  of  4  «  -20  log  a  as  a 
function  of  3. 


values  less  than  ir/2  being  five  times  more  likely  than  values  greater 
than  v/2.  Thus,  we  have  the  probability  density  function  per  degree 
as: 


pM) 


l 

216 


|*|  <90° 


1 

1080 


90°  s  |  *  |  S  180. 


(7) 


The  functions  in  eqs.  (5)  to  (7)  can  be  used  to  determine  the 
probability  of  finding  a,  b,  and  fa  in  some  region  of  a-b-f0  space.  This 
probability  can  be  converted  to  number  of  seconds  in  the  observation 
period  by  multiplying  by  5400  seconds.  To  convert  this  probability  to 
the  number  of  seconds  in  a  month  requires  scaling  the  data  base.  The 
scaling  may  be  obtained  from  Fig.  5  which  shows;  for  several  frequen¬ 
cies  in  the  band,  the  time  during  the  model  data  base  period  that  the 
channel  was  faded  below  a  given  leveL  Distributions  are  shown  for 
average  power  in  the  band  and  for  power  at  selected  frequencies  at  the 
center  and  near  the  edges  of  the  radio  channel  (Frequencies  indicated 
are  at  if  where  the  center  frequency  is  at  70  MHz.)  For  the  path  used. 
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one  expects  the  received  voltage  (single  frequency)  relative  to  midday 
average  to  be  less  than  L  for  72.5  x  104  L'  seconds  in  a  month/  For  the 
data  base  used,  the  fading  is  best  described  by  48  x  104  L*  hence,  the 
data  represent  %  of  a  fading  month.  To  obtain  outage  on  a  seconds- 
per- heavy -fading- month  basis,  the  probabilities  calculated  with  eqs. 
(5)  to  (7)  must  be  multiplied  by  5400  x  1.5  or  8100  seconds. 

2.3  Outage  intimation 

The  fixed  delay  model  described  above  can  be  simulated  with  an 
equivalent  circuit  laboratory  measurement  to  determine  the  equip¬ 
ment  response  to  multipath  fading.  Conceptually,  one  determines 
critical  values  of  A  and  B  for  which  a  specified  error  rate  is  achieved 
for  each  fade  notch  position.  In  practice,  it  is  difficult  to  maintain  a 
constant  BEK  it  is  more  expedient  to  fix  b  and  vary  the  carrier-to-noise 
ratio  (a)  while  plotting  the  ber.  From  the  resulting  curves,  one  may 
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compute  critical  contours  of  A  and  B  for  each  prescribed  notch  location 
and  beb.  Using  eqs.  (5)  and  (6),  the  probability  that  A  and  B  lie  on  the 
high  error  rate  side  of  a  given  critical  contour  may  be  calculated. 

By  repeating  this  calculation  for  a  uniform  set  of  notch  positions 
and  using  (7)  to  determine  the  probability  weighting  given  to  each  and 
summing,  one  may  estimate  the  probability  of  all  selective  fades  that 
produce  a  bes  exceeding  the  prescribed  one.  Multiplying  this  proba¬ 
bility  by  5400  gives  the  outage  time  expected  over  the  data  base  period; 
multiplying  by  3100  gives  the  expected  outage  time  per  heavy  fading 
month. 

The  following  section  describes  the  laboratory  measurement;  Sec¬ 
tion  IV  describes  the  reduction  of  the  measured  curves  and  parameters 
to  outage  times. 

III.  LABORATORY  MEASUREMENTS 

Figure  6  illustrates  stressing  of  a  digital  radio  system  by  means  of  an 
if  fade  simulator.  The  simulator,  which  is  inserted  after  linear  if 
preamplification  but  before  any  high-gain  amplification,  shapes  both 
the  desired  signal  and  the  effective  received  noise.  It  is  necessary  to 
operate  the  simulator  at  a  sufficiently  large  input  carrier-to-noise  ratio 
that  the  concomitantly  shaped  noise  at  its  output  remains  a  negligible 
contributor  to  degraded  system  performance  throughout  the  operating 
range  of  interest. 

Within  its  restricted  frequency  range  of  operation,  the  if  simulator 
is  adjusted  to  achieve  those  specific  shapes  implied  by  Fig.  2.  Although 
the  measurements  could  have  been  made  using  an  RF  fade  simulator, 
the  choice  of  an  IF  simulator  was  based  primarily  upon  considerations 
of  signal  and  noise  levels,  and  the  repeatability  of  adjustments.  The 
following  section  describes  an  if  shape-stressing  measurement  in  the 
minimum  detail  necessary  to  qualify  the  data  collected. 

3. 1  Rapnaantatlva  IF  two-path  fad a  ttraamnq  maaauramant 

The  block  diagram  of  Fig.  7  illustrates  an  arrangement  employing 
an  if  fade  simulator  and  an  if  flat  noise  source.  A  pseudo-random  test 


Fig.  S — it  fade  iireeeini 
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pattern  modulates  the  6-GHz  radio  transmitter  whose  output  is  nom¬ 
inally  5  watts  ((g),  in  Fig.  7).  The  output  spectrum  is  usually  shaped  by 
a  bandpass  filter  following  the  microwave  power  amplifier  to  comply 
with  FCC  regulations. 

To  enable  back-to-back  operation  of  the  transmitter  and  receiver  of 
a  single  repeater  which  normally  operate  on  different  radio  channels, 
a  radio  test  translator  was  employed.  The  translator  output  power  was 
approximately  -30  dBm  (adjustable,  at  (B>)  to  simulate  the  unfaded 
received  signal  level  (rsl)  observed  typically  in  the  field. 

Assuming  a  linear  RF-IF  conversion  gain  of  20  dB,  the  signal  power 
at  the  input  (©)  to  the  if  fade  simulator  is  -10  dBm.  The  simulator 
incorporates  low-noise  linear  amplification.  A  reference  insertion  loss 
for  the  main  unfaded  ray  is  10  dB,  including  the  output  power  summer. 
Hence  the  mmrimiiin  desired  signal  power  at  the  input  to  the  main  if 
amplifier  (©)  is  —20  dBm. 

Assuming  a  30-MHz  receiver  noise  bandwidth  and  a  current-art 
receiving  system  noise  figure  of  S  dB,  the  total  system  noise  power  is 
approximately  —95  dBm,  referred  to  the  receiver’s  input  port.  This 
results  in  a  flat  receiver  noise  contribution  of  -85  dBm  at  input  ©  to 
the  main  if  amplifier.  Consequently,  the  maximum  attainable  carrier 
to  simulator-shaped  rf  noise  ratio  is  10  log  iCa/Nrf)  “  -20  -  (-85)  » 
65  dB.  The  noise  contributed  by  the  fade  simulator  amplifiers  must 
not  exceed  -100  dBm,  to  be  negligible. 

Flat  if  noise  much  larger  than  the  unwanted  and  shaped  system 
noise  is  added  artificially  at  ©  and  is  adjusted  in  magnitude  by  a 
calibrated  attenuator  ®  to  superpose  thermal  noise  degradations  upon 
the  simulated  selective  fading  degradations  of  the  desired  signal.  One 
would  ideally  measure  the  added  if  noise  power  in  the  final  predetec¬ 
tion  bandwidth  of  the  digital  radio  system,  or  twice  the  Nyquist 
bandwidth.  It  is  more  convenient  in  the  laboratory  to  reference  carrier- 
to-noise  ratios  to  the  output  of  the  main  IF  amplifier  by  using  the 
precalibrated  acc  voltage  (assuming  that  wideband  agc  detection  is 
employed),  to  measure  both  the  unshaped  signal  and  noise  powers. 
The  carrier-to-noise  ratio  at  the  detector  would  be  higher — by  the 
ratio  of  the  system  noise  bandwidths  that  would  be  measured  at  the 
respective  points. 

The  noise  source  output  in  Fig.  7  may  be  adjusted  so  that  an 
attenuator  setting  of  0  dB  3)  results  in  a  noise  power  equivalent  to 
that  of  the  unfaded  signal  power  (the  attenuator  is  then  calibrated 
directly  in  uncorrected  Co/fV,/,  in  dB).  As  the  if  fade  simulator  is 
readjusted  to  achieve  different  prescribed  fade  shapes,  its  mean  inser¬ 
tion  loss  may  also  change.  The  change  in  insertion  loss  is  determined 
by  monitoring  the  change  in  signal  power  at  ©;  the  same  loss  increment 
(dB)  must  be  added  to  the  noise  attenuator  to  reestablish  the  0  dB 
reference. 
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3.2  IF  two-path  tada  simulator 

Figure  3  illustrates  splitting  the  desired  if  signal  into  an  arbitrarily 
phased,  adjustable  “main”  component  and  a  “delayed”  component 
fixed  in  delay  (r  ns)  but  adjustable  in  magnitude.  The  main  component 
is  further  resolved  into  orthogonal  components  (inset  to  Fig.  3)  using 
wideband  networks  exhibiting  flat  gain  and  well-behaved  delay.  A 
particular  sum  vector  is  constructed  by  adjusting  the  orthogonal 
components  to  establish  a  simulated  fade  notch  frequency,  in  practice, 
the  phase  sense  of  0-  and  90-degree  components  are  independently 
reversible,  as  indicated  by  the  switches  in  the  figure,  for  complete 
flexibility  in  notch  frequency  selection. 

The  6.3-ns  fixed  delay  added  to  the  delayed  path  imparts  a  phase 
shift  of  159  degrees  at  the  TO- MHz  if  center  frequency.  This  is  shown 
built  out  to  225  degrees,  relative  to  the  0-degree  transmission  path, 
using  a  66-degree  wideband  network  of  the  same  type.  The  delayed 
vector  is  fixed  in  direction  opposite  the  midrange  position  of  the 
adjustable  main  vector,  corresponding  to  a  channel-centered  fade  (<t> 
—  0  degrees). 

Since  1/r  -  158.4  MHz,  a  change  of  1  degree  in  direction  of  the  main 
vector  corresponds  to  a  frequency  displacement  of  the  fade  notch 
location  of  0.44  MHz.  For  <t>  »  -45  degrees,  the  notch  location  is 
displaced  19.8  MHz  below  the  channel  center  {fo  m  -19.3  MHz).  The 
magnitude  of  the  delayed  component  is  then  adjusted  to  achieve  the 
desired  notch  depth. 


3  -  (-of 


OSLAYtD 

Fi*.  3 — it  fad*  simulator— conceptual  block  diagram. 
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3.3  Digital  radio  performance  atraaaad  by  inland  selectivity  and  thermal 
noiaa 

The  radio  equipment  was  measured  at  uniformly  spaced  notch 
frequencies  separated  by  4.4  MHz  (d<p  -  10  degrees).  To  fully  charac¬ 
terize  a  period  of  variation  in  <t>,  or  one  would  need  to  make  36  sets 
of  measurements.  Ideally,  half  may  be  omitted  because  of  symmetry. 
For  given  values  of  A  and  B,  the  same  error  rate  should  obtain  for  a 
notch  at  a  given  frequency  displacement  above  or  below  the  channel 
band  center.  Variations  in  B  ought  not  to  have  a  significant  effect  for 
|  <t>  |  greater  than  90  degrees.  It  was  determined  that  detailed  measure¬ 
ments  were  required  for  nine  different  values  of  f,  to  characterize  the 
digital  radio  tested. 

Using  a  wideband  rf  fade  simulator  in  the  field,  the  digital  radio 
performance  for  out-of-channel  notch  locations  was  relatively  inde¬ 
pendent  of  whether  minimum  or  noruninimum- phase  fade  simulations 
were  employed.  The  nonminimum  phase  fade  is  modeled  by  eq.  (1) 
with  the  sign  of  the  phase  term  reversed.  This  leaves  the  amplitude 
[eq.  (2)]  unchanged,  but  reverses  the  sign  of  the  envelope  delay 
distortion  [eq.  (3)].  We  conclude  that  the  minimum  phase  channel 
model  is  sufficiently  general  for  use  in  simulating  the  channel  and  in 
estimating  performance. 

The  if  fade  simulator  was  adjusted  for  each  notch  frequency,  and 
the  depth  of  notch  was  varied  by  adjusting  the  magnitude  of  the 
delayed  component.  Then  various  amounts  of  if  thermal  noise  were 
added.  Figure  9  typifies  the  performance  data  collected,  bers  are 
plotted  versus  the  uncorrected  if  carrier-to-noise  ratio  (C/Ntf),  for  a 
constant  fade  notch  offset  from  midchannel  (f0  m  -19.8  MHz,  r  -  6.3 
ns).  Each  curve  corresponds  to  a  different  notch  depth  (B  —  -20 
log(l— b)  dB),  and  hence  a  different  amplitude  and  delay  shape  in  the 
radio  channel.  Each  curve  is  also  identified  with  an  in-band  selectivity, 
defined  as  the  difference  between  the  maximum  and  minimum  atten¬ 
uation  present  in  the  (25.3-MHz)  channel  bandwidth.  The  lower-left 
“baseline”  curve  presents  the  unshaped  signal,  flat  fading  performance 
obtained  by  adding  only  if  thermal  noise.  This  curve  was  verified 
(without  the  added  if  noise)  by  attenuating  the  received  rf  input 
signal  in  the  back-to-back  configuration. 

Consideration  was  given  to  matching  the  order  of  measurements  to 
characteristics  of  the  particular  digital  radio  tested.  For  example, 
considerable  scattering  of  data  at  low  error  rates  can  resuit  from 
synchronizations  involving  different  reference  carrier  phases.  The  au¬ 
thors  elected  to  perform  several  synchronizations  while  observing  the 
ber  for  each  phase,  and  then  chose  that  phase  giving  the  worst 
performance.*  Synchronization  was  accomplished  at  the  low  error  rate 

*  Became  the  phase  information  in  the  measured  system  was  Gray  coded  and  digital 
access  was  on  a  per-raii  basis,  one  rail  had  twice  the  aut  of  the  other  two.  Ail 
measurements  in  the  field  and  in  the  laboratory  were  made  on  this  worst-case  raiL 
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Fig.  9 — High-ipaod  digital  radio  if  dispersive  fad*  simulation*,  r  -  6.3  ns,  f„  »  - 19.8 
MHi.  j 

(bottom)  of  each  curve,  and  this  phase  relationship  was  maintained 
for  all  data  points  obtained  for  each  curve. 

From  the  baseline  curve  of  Fig.  9,  a  BEK  -lx  10**  obtains  for  10 
log (C/N,f)  *  21.5  dB.  For  the  digital  radio  system  installed  on  the 

instrumented  hop  and  reported  in  the  figure,  the  measured  flat  fade  I 

margin  for  a  threshold  BEK  -lx  10"*  was  40.5  dB  This  leads  to  an  I 

unfaded  if  carrier-to-noise  ratio  10  log(C0/N,/)  *  21.5  +  40.5  -  62  dB. 

From  the  baseline  curve  for  a  BEK  -  1  x  10*?,  note  that  insertion  of 
a  fade  whose  notch  depth  is  6.5  dB  results  in  four  orders  of  magnitude 
degradation  in  ber  performance;  equivalently,  an  in-band  selectivity 
of  only  5.7  dB  in  25.3  MHz  results  in  a  BER  >  1  x  10"3. 

The  asymptotic  regions  in  Fig.  9,  corresponding  to  high  values  of  Cl 
N,f,  are  not  normally  presented  in  characterizations  of  this  type: 
however,  system  outage  depends  primarily  upon  the  performance  in 
these  asymptotic  regions.  Thus,  under  typical  fading  conditions,  the 
transmitted  carrier  power  might  be  increased  at  will  without  improving 

the  ber  significantly.  The  effects  of  decreasing  the  carrier  power  are  ' 

discussed  in  Section  4.4.  * 

A  family  of  curves  like  those  shown  in  Fig.  9  was  obtained  (but  are  i 

not  given  here)  for  each  of  nine  uniformly  spaced  frequency  offsets 
below  midchannel  to  characterize  the  digital  radio  system  sufficiently 
for  the  prediction  of  outage.  A  number  of  spot  checks  were  also  made 
using  both  rf  and  if  fade  simulators  at  symmetrical  positive  and 
negative  offset  frequencies,  to  establish  that  acceptable  symmetry 
existed. 
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IV.  CALCULATION  OF  3U7AGE 


This  section  describes  four  methods  of  calculating  outage.  The 
derivation  of  the  critical  curves  of  .4  and  5.  which  provide  the  basis 
for  making  and  understanding  these  calculations,  is  given  in  Section 
4.1.  In  Section  4.2  the  detailed  calculation  of  outage  from  the  critical 
,4-fl  curves  is  described.  It  is  shown  in  Section  4.3  that  for  the  present 
system  this  method  may  be  greatly  simplified  by  calculating  only 
selectivity-caused  outage  (i.e.,  neglecting  thermal  noise).  Section  4.4 
presents  an  approximate  method  of  accounting  for  the  effects  of 
thermal  noise.  Section  4.5  provides  a  basis  for  estimating  the  selectiv¬ 
ity-caused  outage  from  a  single  measurement. 

4. 1  Derivation  of  critical  characteristics 

To  calculate  the  outage  for  a  fixed  bit  error  rate,  one  must  first 
obtain  the  critical  curves  of  .4  and  3  at  each  simulated  value  of  f0,  the 
notch  position.  Thus,  from  Fig.  9  which  corresponds  to  f„  «•  —19.8  MHz 
; or  <t>  —  —45°),  we  obtain  six  points  on  the  critical  curve  of  .4  and  B  for 
a  ber  of  10-J,  one  point  from  each  of  the  six  curves  which  cross  the 
critical  ber.  The  value  of  B  is  obtained  from  the  value  of  b  since 

B  -  -20  log(I  -  b ).  (8) 

For  the  curve  in  Fig.  9  corresponding  to  B  «■  4.4  dB,  we  obtain  the 
corresponding  critical  value  of  .4  for  a  BER  of  10~}  from  the  value  of 
carrier-to-noise  ratio,  which  is  20.2  dB  where  this  curve  crosses  the 
10~3  ber  line.  Since  the  earner- to- noise  ratio  is  62  dB  when  the  channel 
is  unfaded,  the  20.2  dB  value  corresponds  to  a  relative  average  power 
loss  of  41.8  dB, 


L.  -  62  -  20.2  -  41.8  dB.  (9) 

Without  loss  of  generality,  we  assume  that  the  PSK  signal  has  a 
rectangular  spectrum  of  width  fa  consequently,  the  relative  power 
transmitted  by  the  model  is  obtained  from  eq.  (2)  as* 


(10) 


*  The  calculated  result  ia  not  critically  dependent  on  the  flatness  of  the  signal 
spectrum  or  the  spectral  width  chosen.  We  have  used  for  ft  a  value  of  25.3  MHz  aa 
representing  the  effective  width  of  the  signal 
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Defining  a  correction  term  by 

C  —  -10  log'  1  +  b2  -  2b  cos  | ' 


we  may  express  the  signal  loss  as 

L.  -  -10  log  Pan  *■  A  +■  C. 
Thus,  we  obtain  the  critical  value  of  .4  as 


A-L.-C. 


(11) 


(12) 

(13) 


For  B  —  4.4  dB  (6  —  0.4)  and  f0  m  —19.8  MHz,  we  find  C  —  2.06  dB  and 
the  critical  value  of  A  is  41.8  —  2.1  —  39.7  dB. 

Carrying  out  these  calculations  for  the  six  curves  in  Fig.  9,  one  can 
generate  the  critical  curve  of  .4  and  3  for  f0  ”  —19.8  MHz  and  a  ber 
of  10"3.  The  curve  is  shown  in  Fig.  10  along  with  the  critical  curves  for 
several  other  values  of  the  BER.  A  complete  set  of  curves  must  be 
generated  for  all  values  of  /<* 

The  curves  in  Fig.  10  are  typical  of  the  critical  curves  obtained  for 
|  f0 1  3  33  MHz.  The  intercept  with  the  A-axjs  represents  the  flat  fade 
margin  for  the  given  ber;  this  margin  is  independent  of  notch  position. 
The  intercept  of  a  critical  contour  with  the  B-axis  represents  the 
shape,  or  relative  fade  depth,  margin  for  the  given  notch  position.  For 
values  of  B  to  the  right  of  this  intercept,  the  critical  value  of  ber 
cannot  be  obtained  at  any  carrier-to-noise  ratio  for  the  given  notch 
position. 

4.2  Outage  calculation —  detailed  method 

The  probability,  P,,  of  finding  A  and  B  outside  all  critical  contours 
may  be  written  with  eqs.  (3),  (6),  and  (7)  as 


P„  -  |  p.(<t>)PA<t>)  Cfd,  (14) 

where 

Pc(<f>)  - 

and  A*(X)  is  the  functional  relation  of  the  critical  values  of  A  to  B  (or 
X),  for  B  less  than  Bc,  the  B-axis  intercept,  and  for  a  given  ber  and  <t> 
value.  *  Since  measurements  were  made  for  a  uniformly  spaced  set  of 
notch  positions  with  spacing  $4>  “  10°,  we  may  approximate  (14)  by 


*  The  dependence  of  the  function  and  the  asymptote  Br  on  asa  and  a  is  not 
explicitly  denoted  to  keep  notation  simple. 


IT  ■ 

Jo  JA'ixi 


pA(Y)pa(X)  dYdX, 


(15) 
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Fig.  10 — Critical  curves  of  A  and  B  for  f,  -  - 19.8  MHz. 

Po-bt>  z  p*(<t>i)PA<t>i).  (16) 

Ail 

To  illustrate  the  calculation  of  outage  probability  with  eqs.  ( 15)  and 
(16),  we  shall  calculate  the  term  in  the  summation  of  eq.  (16)  corre¬ 
sponding  to  a  BER  of  10"’  and  ■  -45°  (or  fo  m  -19.8  MHz).  From 
Fig.  11,  which  is  taken  from  Fig,  10,  we  note  that  the  double  integral 
in  eq.  ( 15)  may  be  broken  into  integrations  over  two  regions.  Thus 


—LI 


Pa(Y)pb(X)  dYdX 


IT 

JO  JAr 


Pa(Y)pb(X)  dYdX, 


where  the  two  double  integrals  correspond  to  integrations  over  Regions 
1  and  2,  respectively,  in  Fig.  11.  Outage  due  to  the  occurrence  of  A  and 
B  in  Region  1  may  be  described  as  outage  due  only  to  shape  or 
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selectivity.  In  Region  2,  outage  is  due  to  the  combined  effects  of  signal 
loss  and  selectivity. 

Using  eqs.  (5)  and  (6),  the  integral  over  Region  I  is  obtained  as 
The  contribution  due  to  thermal  noise  and  shape  (Region  2)  is 
slightly  more  complicated.  Dividing  the  interval  0  to  Br  in  Fig.  11  into 
N  subintervals,  as  shown  in  Fig.  12,  the  probability  of  being  in  Region 
2  is  the  sum  of  the  probabilities  for  each  subintervaL  Thus  eq.  (17) 
becomes 

PA4>.)  -  e~B'/3a  +  V  [e-"-'3*  -  e-a>'3*]P,  f-— ~A0-i3 - j ,  (18) 

where 

pg(X)  -  -4=  f  dx.  (19) 

V2ff  Jx 


Fig.  11 — Claaeiflcation  of  outage  with  respect  to  critical  curve  for  But  ”  10' '.  fo  - 
19.8  MHz. 
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3  •  -20  Log  U-OI  IN  DECIBELS 


Fig.  12 — Outage  calculation  for  an  incremental  interval 


Evaluating  the  two  components  of  eq.  (IS)  from  Fig.  11,  we  find 

Pc(— 45°)  -  0.181  +  0.003  -  0.184.  (20) 

This  calculation  was  performed  for  10  values  of  from  —5  to  —85 
degrees  in  10-degree  steps.  Using  these  results  in  eq.  (16)  and  multi¬ 
plying  by  two  to  account  for  positive  values  of  <t>,  which  are  assumed  to 
contribute  equally,  we  find  the  probability,  P„  for  a  BER  of  10'3  as 

Pa  -  0.0996. 

The  expected  outage  for  the  data  base  period  is,  then. 

To  m  5400  x  0.0996  »  538  seconds.  (21) 

4.3  Outage  calculation —  selectivity  only 

It  is  apparent  firom  eq.  (20)  that  most  of  the  outage  for  the  system 
under  study  is  caused  by  selectivity,  fades  characterized  by  A  and  B 
values  in  Region  1.  From  eqs.  (14)  and  (IT),  we  may  express  P«,  the 
probability  of  outage  due  to  selectivity,  as 
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(22) 


Pt(<t>)pa(X)  dX  d<t>. 


For  the  system  studied  for  a  ber  of  10-1,  a  finite  Be  is  obtained  only 
for  |  ®i  |  <  90  degrees.  Hence,  we  may  use  eq.  (7)  to  simplify  (22)* 


Pc 


'^Ze- 

216 


■8,14,1/3.8 


(23) 


From  eq.  (23)  we  see  that  the  outage  due  oniy  to  selectivity  depends 
on  the  relationship  between  Bc,  the  asymptote  of  critical  B  values,  and 
the  notch  angle  or  notch  frequency.  Figure  13  shows  the  relationship 
between  B?  and  the  notch  frequency  for  four  values  of  ber.  It  is 
apparent  from  eq.  (22)  that  the  outage  probability  is  the  probability  of 
finding  B  and  f0  values  in  the  region  above  this  curve.  Such  curves, 
therefore,  provide  a  useful  basis  for  evaluating  the  selectivity  outage 
of  a  digital  radio  system. 


4.4  Outage  calculation—  approximata  method 

For  a  radio  system  sensitive  to  both  thermal  noise  and  selectivity, 
the  calculation  of  Section  4.3  is  inadequate  and  that  of  Section  4.2  is 
unduly  cumbersome. 

To  illustrate  a  simpler,  but  generally  applicable,  method  and  at  the 
same  time  to  provide  a  useful  incidental  result,  let  us  evaluate  the 
effect  of  reducing  the  transmitted  power  by  10  dB.  For  the  reduced 
power  system,  the  carrier-to-noise  ratio  would  be  52  dB  for  the  unfaded 
channel,  and  the  critical  curves  of  A  and  B  would  be  shifted  by  10  dB. 
Figure  14  shows  the  critical  curve  of  A  and  B  for  a  10"’  ber  and  f„  » 
— 19.3  MHz  with  an  overplot  of  the  conditional  mean  of  the  distribution 
of  A.  The  dotted  curves  on  Fig.  14  represent  2-sigma  intervals  on 
either  side  of  the  mean.  From  the  properties  of  the  Gaussian  distri¬ 
bution,  one  may  determine  that  more  than  95  percent  of  the  values  of 
A  and  B  will  lie  between  these  two  dotted  curves.  We  designate  as  Am 
and  Bm  the  coordinates  of  the  intersection  of  two  curves:  the  critical  A 
—  B  curve  and  the  conditional  mean  curve.  Then  approximating  the 
critical  curve  of  A  and  B  with  a  straight  line  segment  tangent  at  (Am, 
Bm),  with  slope  s,  we  may  approximate  the  probability  of  outage  by 
integrating  the  probabilities  over  the  region  to  the  right  of  the  tangent 
line.  Using  eqs.  (15)  and  (16),  we  obtain 


Po 


'  Ad>  T 

Ml  * 


P*(<t>i) 


n 

Jo  J* 


pA(Y)p„(X)dY  dX.  (24) 


0  J 


’  The  factor  of  two  is  required  in  eq.  (23)  because  the  indicated  summation  corre¬ 
spond!  to  an  integration  only  over  negative  notch  frequencie*  i«,  <  0). 
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Fig.  13— Asymptotic  performance  curves.  Locus  of  values  of  f>  and  B  that  produce  a 
fixed  ber  at  high  carrier- to- noise  ratio. 


Interchanging  the  order  of  integration  and  ignoring*  the  B  dependence 
of  Ao(£),  this  becomes 

Po-ht>  l  (25) 

All  a, 

Evaluating  eq.  (25)  for  a  10_;l  ber  and  multiplying  the  result  by  5400 
gives  an  outage  estimate  for  the  data  base  period  of  602  seconds. 
Recalculating  the  total  outage  time  at  a  10~J  ber  using  the  method  of 
Section  4.2  [eqs.  (16)  and  (18)]  gives  636  seconds,  which  verifies  the 
accuracy  of  the  approximate  method.  The  estimate  of  636  seconds  was 
calculated  as  an  upper  bound;  the  602  seconds  calculated  using  (25) 
tend  to  be  a  lower  bound.  We  conclude  that  backing  off  transmitted 
power  by  10  dB  would  increase  the  outage  by  about  12  percent  (538  to 
602). 

4.5  A  further  simplification 

In  this  section,  we  show  that  the  outage  due  to  selectivity  can  be 
estimated  approximately  for  a  given  ber  from  a  determination  of  the 
in-band  selectivity  required  (with  the  notch  out-of-band)  to  produce 
that  ber.  Such  a  measurement  may  provide  a  useful  approximation 
for  any  digital  system  using  quadrature  modulation  components;3 
however,  we  provide  a  justification  based  on  the  performance  of  the 
system  at  hand.  In-band  selectivity  is  defined  as  the  difference  between 
the  maximum  and  minimum  attenuation  present  in  the  (25.3-MHz) 
channel  bandwidth. 

Since  the  in-band  selectivity  is  a  constant  for  any  of  the  curves 
shown  in  Fig.  9,  one  can  use  Fig.  9  to  plot  the  asymptotic  ber  against 
in-band  selectivity  for  f0  m  —19.8  MHz.  Such  a  plot  was  generated  for 
each  notch  position  measured  to  produce  the  family  of  curves  shown 
in  Fig.  15.  Note  that,  except  for  notch  positions  near  the  band  center, 
the  BER  is  uniquely  related  to  the  in-band  selectivity.  Neglecting  the 
in-band  notches,  we  find  that  a  10' 1  ber  corresponds  to  an  in-band 
selectivity  of  5.5  dB. 

If  we  use  the  original  model  of  eq.  (2)  to  determine  the  values  of  B 
that  will  produce  an  in-band  selectivity  of  5.5  dB  for  a  number  of 
different  notch  positions,  we  would  generate  Fig.  16.  It  is  apparent  that 
for  this  system  there  is  a  good  correspondence  between  the  curves  of 
asymptotic  performance  (Fig.  13)  and  the  curves  of  constant  in-band 
selectivity  (Fig.  16). 

To  reinforce  this  conclusion,  we  provide  Figs.  17,  18,  and  19.  Figure 
17  shows  the  locus  of  in-band  selectivity  in  a  25.3-MHz  band  corre- 

*  Including  the  effect  of  the  slop*  of  An(B)  at  B  ■  Bm  give*  the  same  symbolic  result 
with  s  interpreted  as  the  algebraic  sum  of  the  slop*  of  the  tangent  and  dAo/dB  evaluated 
at  B  -  Bm. 
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Fig.  15— Measured  asymptotic  bit  error  rate  vs  peak-to-peak  amplitude  difference  in 
a  25.3-MHz  band. 


Fig.  16 — Locus  of  3  and  U  for  modeled  fade*  that  have  fixed  peak-to-peak  amplitude 
in  a  25.3-MHz  band. 


sponding  to  each  of  the  curves  of  constant  ber  in  Fig.  13.  That  is,  for 
each  ber  and  each  value  of  notch  position,  f„,  we  have  plotted  the 
peak-to-peak  amplitude  difference  in  the  band  for  the  corresponding 
value  of  Be,  the  asymptotic  critical  value  of  B.  Figure  IS  shows  a 
similar  set  of  curves  with  the  peak-to-peak  delay  distortion  in  a  25.3- 
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V  DISPLACEMENT  Of  NOTCH  IN  MEGAHERTZ 


Fig.  17 — in-band  selectivity  (in  23.3- M Hi  bandwidth)  corresponding  to  asymptotic 
critical  values  of  notch  depth  (B,)  for  several  values  of  BEJL 


DISPLACEMENT  OP  NOTCH  IN  MEGAHERTZ 


Fig.  18— Peak-to-peak  envelope  delay  distortion  in  23.3-MHt  bandwidth  correspond¬ 
ing  to  asymptotic  critical  values  of  notch  depth  (B.)  for  several  values  of  BEK. 

MHz  band  as  the  ordinate.  Similarly,  Fig.  19  has  as  the  ordinate  the 
“slope,”  or  amplitude  difference  at  a  separation  of  25.3  MHz.  It  is  again 
clear  from  these  three  figures  that  the  in-band  selectivity  is  the 
relevant  channel  impairment  giving  rise  to  errors.  We  see  from  Fig.  IS 
that,  for  out-of-band  notches,  high  bebs  are  obtained  with  very  small 
values  of  peak-to-peak  delay  distortion,  and  from  Fig.  19  that  for  in- 
band  notches  high  bebs  are  obtained  for  very  small  values  (zero  at 
mid-band)  of  slope. 
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Fig.  19 — Amplitude  difference  at  a  25.3-MHz  frequency  separation  corresponding  to 
asymptotic  critical  values  of  notch  depth  (&)  for  several  values  of  bex. 

The  model  data  base  was  analyzed  to  determine  the  time  during 
which  the  in-band  selectivity  in  a  band  of  25.3  MHz  exceeded  a  given 
value.  Figure  20  presents  this  distribution  for  in-band  selectivity  as 
calculated  from  the  modeled  fits.  Figure  20  can  be  used  directly  in 
conjunction  with  Fig.  15  to  calculate  the  outage  times  for  the  model 
data  base.*  For  instance,  from  Fig.  15  we  note  that  5.5  dB  of  selectivity 
corresponds  to  a  10~3  BEK.  We  use  Fig.  20  to  determine  that  5.5  dB 
was  exceeded  for  520  seconds. 

V.  COMPARISONS  OF  CALCULATED  AND  OBSERVED  OUTAGES 

Using  the  methods  of  Sections  4.2  to  4.5,  outage  times  were  calcu¬ 
lated  for  bit  error  rates  of  10~3  to  10"®  for  both  the  model  data  base 
period  and  for  a  heavy  fading  month,  by  multiplying  the  outage 
probabilities  by  5400  and  8100,  respectively. 


5. 1  Modal  data  baaa  partod 

Calculated  and  observed!  outages  for  the  model  data  base  period 
are  shown  in  Table  I.  In  general,  comparing  the  calculated  results  with 
observed  results,  we  see  that  the  outage  is  underestimated  at  high 
beks  and  overestimated  at  low  bees.  Any  estimation  procedure  based 
on  the  current  modeled  state  of  the  channel  will  tend  to  underestimate 


*  In  practice,  on*  would  um  a  single  measurement  of  in-band  selectivity.  For  instance, 
in  Fig.  9  one  would  take  the  5.7-dB  value  corresponding  to  the  curve  asymptotic  at  a 
1C-'  a  xx. 

t  Because  of  quantization,  the  outage  tunes  observed  from  the  field  experiment 
correspond  to  bit  error  ratee  of  1.28  x  10"'.  1.57  x  10'*.  0.981  x  to-',  and  1.19  x  10-” 
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Table  I— Outage  in  modeling  data  base  period  (seconds) 


BCR  -  10" 

10" 

10" 

10" 

Observed 

636 

903 

1191 

1487 

Detailed  calculation  (Section  4.2) 

538 

960 

1430 

1860 

Approximate  calculation  (Section  4.4) 

527 

950 

1420 

1830 

Asymptotic  calculation  i  Section  4.3) 

527 

950 

1420 

1830 

Selectivity  calculation  (Section  4.5) 

510 

900 

1570 

2730 

outage  at  high  bess  because  of  hysteresis  effects  in  the  radio  receiving 
equipment.  That  is,  when  the  channel  condition  becomes  sufficiently 
severe,  the  bit  error  rate  becomes  high  enough  (on  the  order  of  10~3) 
that  the  timing  and/or  phase  of  the  radio  system  breaks  lock.  If  the 
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channel  impairment  becomes  less  severe,  the  ber  will  not  improve 
until  the  system  resynchronizes.  The  hysteresis  is  important  at  the 
10“3  ber,  since  a  significant  fraction  of  the  events  that  cause  10'  *  ber 
will  cause  the  system  to  break  lock. 

One  would  expect  to  overestimate  the  outage  at  low  bers  because 
of  the  method  of  taking  data.  Recall  that,  in  measuring  the  curves  in 
Fig.  9,  it  was  found  that  the  ber  depended  on  the  phase  to  which  the 
system  had  locked.  The  recorded  performance  represented  the  worst- 
phase  condition.  At  a  10~'  ber,  the  best  phase  produces  a  ber  that  is 
about  Hi  that  produced  by  the  worst  phase;  the  difference  in  ber  from 
worst  to  best  phase  at  a  10'*  ber  is  negligible.  Hence,  one  would 
expect  outage  to  be  overestimated  significantly  at  low  bit  error  rates. 

In  comparing  the  outage  calculated  from  in-band  selectivity  (Section 
4.5)  to  the  outage  observed,  we  find  that  the  overestimation  of  outage 
at  low  bers  is  more  severe  than  with  the  other  methods.  This  is  due 
to  the  greater  sensitivity  of  the  differential  selectivity  method  to  the 
bias  induced  by  choosing  the  worst-case  phase.  For  instance,  compar¬ 
ing  calculations  at  a  10'*  ber,  we  find  that  Fig.  20  is  steeper  for 
amplitude  differences  near  2  dB  than  is  Fig.  3  near  B  values  of  3.5  dB. 
(Figure  9  verifies  the  appropriateness  of  this  comparison).  More  gen¬ 
erally,  one  expects  the  method  based  on  in-band  selectivity  to  overes¬ 
timate  the  outage  because  the  method  is  based  on  notches  out  of  band. 
From  Fig.  15,  it  is  apparent  that,  for  a  given  M,  some  scans  will  not 
have  the  BER  specified. 

We  conclude  that,  although  calculation  of  outage  from  sensitivity  to 
in-band  selectivity  provides  quick  estimates,  they  are  less  accurate. 
The  calculation  requires  knowledge  of  the  distribution  of  in-band 
selectivity  over  a  specified  bandwidth.  These  statistics  are  neither 
simple  nor  generally  available.  It  has  been  shown,10  for  instance,  that 
slope  statistics  have  a  nontrivial  dependence  upon  the  measurement 
bandwidth. 

It  is  clear  that  the  calculations  based  on  selectivity  (Sections  4.3  and 
4.5)  agree  for  the  system  studied  here  because  that  system  has  very 
little  outage  due  to  thermal  noise  limitations,  and  because  it  is  sensitive 
primarily  to  in-band  amplitude  excursions.  The  extent  to  which  these 
statements  are  true  for  other  systems  is  currently  unknown. 

5.2  Outag e  on  a  monthly  baaia 

The  results  in  Table  I  may  be  put  on  the  basis  of  a  heavy  fading 
month  by  increasing  them  by  a  factor  of  1.5,  as  discussed  in  Section 
2.2.  The  resulting  outages  (including  the  scaled  observed  outage)  are 
compared  with  the  outage  observed  in  a  one-month  period*  in  Table 
II.  We  see  that  the  outage  times  observed  in  the  total  one- month 
period  agree  well  with  the  values  obtained  by  scaling  the  observed 
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Table  II — Outage  in  a  heavy  fading  month  (seconds) 


BER  *  10" 1 

10'* 

10''' 

to" 

Observed  (Rsf.  8) 

1000 

1320 

2100 

2900 

Scaled  observation  from  Table  1 

955 

1350 

1790 

2230 

Calculation  (Section  4.2-4.4J 

800 

1430 

2140 

2760 

Selectivity  calculation  (Section  4.5) 

770 

1350 

2350 

4100 

outage  for  the  data  base  period  used  in  modeling,  except  for  the  slight 
divergence  appearing  at  low  beks.  This  divergence  should  not  be 
unexpected  for  this  equipment.  As  may  be  seen  in  Fig.  15,  a  10~*  BEK 
may  be  caused  by  differential  amplitude  selectivity  in  band  of  2  dB. 
Such  modest  amounts  of  selectivity  may  be  expected  to  occur  some¬ 
times  in  the  presence  of  very  moderate  selective  fading.  The  modeling 
data  base  was  constructed  by  selecting  only  periods  of  significant 
selective  fading.  This  reinforces  the  comments  made  in  conjunction 
with  Fig.  3,  namely,  that  the  model  distribution  of  B  represents  a  lower 
bound  for  small  values  of  B  which  can  contribute  to  outage  at  the  10~* 
ber  level. 

VI.  CONCLUSIONS 

We  have  demonstrated  the  validity  of  a  technique  for  estimating  the 
unprotected  outage  of  a  digital  radio  system  due  to  selective  fading  on 
a  particular  hop  in  the  6-GHz  common  carrier  band.  The  technique 
required  field  measurements  to  statistically  characterize  the  param¬ 
eters  of  a  model  of  propagation  on  the  hop.  It  also  requires  performance 
data  obtained  in  the  laboratory  for  the  radio  system  by  stressing  it 
with  a  two-path  fade  simulator  with  a  differential  path  delay  of  6.3  ns, 
corresponding  to  the  fixed  delay  channel  model.  Since  the  radio  path 
on  which  these  measurements  were  made  has  a  length  close  to  the 
average  for  the  Bell  System  long  haul  radio  network  and  has  an 
average  incidence  of  fading  activity,  the  channel  model  is  representa¬ 
tive  of  a  typical  path.  At  the  very  least,  the  technique  provides  a  basis 
for  determining  the  relative  merits  of  various  digital  radio  systems 
operating  without  benefit  of  space  diversity. 

For  the  system  under  test,  outage  was  calculated  by  four  different 
methods.  Because  this  system  was  selectivity- limited  rather  than 
noise-limited,  all  four  methods  predicted  approximately  the  same 
outage  as  that  summarized  in  Table  I;  however,  the  method  based  on 
in-band  selectivity  is  more  severely  biased  at  low  bers.  The  method 
based  on  asymptotic  performance  and  that  based  on  in-band  selectivity 
can  only  be  used  to  estimate  outage  due  to  selectivity.  If  the  transmit¬ 
ted  power  of  the  system  under  test  were  reduced  by  10  dB,  both  of  the 
other  two  methods,  the  detailed  and  the  approximate  method,  would 
predict  an  increase  in  outage  time  of  about  12  percent. 
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Transmission  Unavailability  of  Frequency- 
Diversity  Protected  Microwave  FM  Radio 
Systems  Caused  by  Multipath  Fading 

By  A.  VIGANTS  and  M.  V.  PURSLEY 
(Manuscript  received  Apnl  2,  1 979) 

Estimates  of  transmission  unavailability  caused  by  multipath  fad¬ 
ing  are  needed  to  determine  the  adequacy  of  diversity  protection 
arrangements.  A  computer  program  producing  such  estimates  was 
made  available  for  general  Bell  System  use  in  1977.  In  this  paper,  we 
summarize  the  underlying  mathematical  model  of  frequency  diversity 
operation  and  atmospheric  multipath  fading  and  use  it  to  study 
frequency  diversity  behavior. 

I.  INTRODUCTION 

Microwave  fm  radio  systems  utilized  in  the  Bell  System  provide 
protected  broadband  communications  channels  that  generally  have  a 
capacity,  depending  on  the  frequency  band  and  on  the  type  of  equip¬ 
ment,  of  1200  to  2400  message  circuits  per  channel  The  protected 
channels,  with  protection  provided  by  radio  switching  systems,  are 
referred  to  as  “working”  channels,  to  distinguish  these  &om  the  actual 
radio  channels,  which  are  assigned  as  either  “regular"  or  “protection” 
channels.  The  number  of  regular  channels  in  a  radio  system  equals  the 
number  of  working  channels.  Equipment  failures,  atmospheric  fading, 
or  other  impairments  in  a  regular  channel  cause  temporary  transfers 
of  the  communications  traffic  to  a  protection  channel  at  a  different 
radio  frequency,  if  an  unoccupied  and  unimpaired  protection  channel 
is  available.  In  the  case  of  fading,  such  frequency-diversity  protection 
switching  capitalizes  on  the  frequency  selectivity  of  the  multipath 
,  fading  process,  where,  at  an  instant  in  time,  the  strength  of  the  received 

signal  is  a  function  of  frequency.  The  Bell  System  long-haul  microwave 
fm  radio  network  contains  approximately  1600  one-way  frequency- 
.i  diversity  switching  sections  (800  route  segments)  with  an  average  of 

three  radio  hops  per  switching  section. 
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The  communications  traffic  in  a  working  channel  can  be  temporarily 
subjected  to  high  noise  when,  in  a  switching  section,  the  number  of 
unpaired  radio  channels  exceeds  the  number  of  protection  channels. 

Bell  System  practice  is  to  describe  this  in  terms  of  the  amount  of  time  | 

during  which  the  noise  in  a  working  channel  exceeds  55  cLBmcO  in  a  | 

message  circuit  at  the  top  frequency  assignment  in  the  baseband.1  In  1 

this  paper,  this  time  is  referred  to  as  “service  failure  time”  in  accord¬ 
ance  with  designations  on  protection  switching  equipment  and  termi¬ 
nology  used  in  the  operation  of  the  plant.  The  terms  “service  failure” 
and  “oucage”  are  frequently  synonymous,1  but  use  of  the  latter  term 
can  sometimes  result  in  semantic  difficulties,  since  some  service  fail¬ 
ures  may  not  result  in  outages  as  perceived  or  defined  by  some 
communications  users. 

In  modem  microwave  radio  systems,  protection  requirements  are 
frequently  governed  by  the  presence  of  multipath  fading, *'IJ  as  opposed 
to  causes  associated  with  equipment  or  human  intervention.1  Esti¬ 
mates  of  service  failure  time  due  to  multipath  fading  are  therefore 
needed  to  determine  the  adequacy  of  diversity  protection  arrange¬ 
ments  when  new  routes  are  planned,  when  transmission  parameters  of  l 

existing  routes  are  changed,  or  when  alternate  protection  strategies  are 
considered  for  future  use.  In  1977,  a  computer  program  producing  such 
estimates  was  made  available  for  general  Bell  System  application.  It 
is  extensively  used  to  treat  everyday  design  and  planning  tasks  on  an 
individualized  basis  for  microwave  FM  radio  routes  equipped  with 
frequency-diversity  and  space-diversity3  protection. 

A  mathematical  model  of  frequency-diversity  operation  and  multi- 
path  fading  is  central  to  the  computer  program.  In  this  paper,  the 

model  is  summarized  and  utilized  to  study  representative  radio  systems  ' 

on  fully  and  partially  developed  routes.  The  effect  of  radio  channels 

with  reduced  fade  margins,  the  service  failure  time  of  individual 

channels,  and  parameters  that  can  be  used  to  characterize  frequency- 

diversity  systems  are  treated.  A  discussion  is  included  of  the  addition 

of  space-diversity  protection  to  frequency-diversity  protection,  often 

necessary  to  reduce  service  failure  time. 

The  numerical  results  presented  in  this  paper  should  be  considered 
in  the  context  of  transmission  performance  objectives.  A  proposed 
design  or  modification  of  a  frequency-diversity  protection  switching 
section  containing  S  radio  hops  is  satisfactory  when 

S 

I  nsu  (l) 

A-( 

where  T*  is  the  estimated  service  failure  time  of  an  average  radio 
channel  due  to  multipath  fading  on  a  particular  radio  hop.  The  service 
failure  time  due  to  simultaneous  deep  fading  on  different  hops  is  much 
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smaller  and  is  therefore  neglected  in  the  above  expression.  The  design- 
objective  time  Tobi  is  obtained  from  the  reliability  objectives  for 
microwave  fm  radio,  which  specify  an  annual  outage  of  0.02  percent  or 
less  due  to  ail  causes  for  a  two-way  circuit  on  a  4000-mile  long-haul  or 
a  250-mile  short-haul  route.1  In  practice,  this  objective  is  prorated  to 
the  actual  distance  and  applied  to  the  service  failure  time.  One-half 
the  objective  is  allocated  to  multipath  fading,  and  the  other  half  to 
equipment,  maintenance,  and  other  causes.9  For  multipath  fading,  one- 
half  the  0.01-percent  two-way  allocation  is  applied  to  each  direction  of 
transmission.  This  accommodates  customers  that  may  transmit  in  only 
one  direction  and  conforms  to  the  structure  of  the  plant,  where 
independent  protection  switching  systems  and  different  radio  frequen¬ 
cies  are  used  in  the  two  directions  of  transmission.  Furthermore,  the 
service  failure  time  associated  with  a  single  protection  switching  sys- 
cem  can  be  expected  to  be  substantially  larger  than  the  service  failure 
time  occurring  simultaneously  in  the  two  directions,  since  the  latter 
kind  of  failure  requires  simultaneous  deep  fading  on  a  larger  number 
of  radio  channels.  Consequently,  the  one-way  multipath-fading  service 
failure  time  objective  for  a  section  of  a  route  containing  S  hops  is,  in 
seconds  per  year. 

To,,,  -  (1600/ A*)  £  Dh,  (2) 

A- 1 

where  D„f  is  4000  miles  for  long-haui  and  250  miles  for  short-haul 
radio,  and  Dh  is  the  hop  length  in  miles.  In  the  case  of  a  frequency- 
diversity  switching  section  containing  a  single  25-mile  hop,  the  values 
of  Tab,  are  10  seconds  per  year  for  long-haul  and  160  seconds  per  year 
for  short-haul  radio.  In  the  past,  it  was  often  assumed  that  only  half 
the  hops  on  a  long  route  would  experience  significant  fading. 19  This 
assumption  is  no  longer  necessary,  since  fading  estimates  can  now  be 
made  for  every  hop  on -a  route. >a 

II.  THE  MATHEMATICAL  MODEL 

The  service  failure  time  of  an  average  working  channel,  caused  by 
multipath  fading  on  a  particular  radio  hop,  can  be  expressed  as 

Th  -  N"Zh,  (3) 

where  N  is  the  number  of  working  channels  and  Zh,  referred  to  as  the 
facility  service  failure  time,  is  the  sum  of  the  service  failure  times  of 
the  working  channels.  The  facility  service  failure  time  can  be  expressed 
as 

.V 

Zh  -  l  iT'(t),  (4) 

■-I 
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where  T’{i)  is  the  accumulated  time  during  which  exactly  t  of  the  N 
working  channels  experience  service  failure.  The  terms  in  this  sum  are 
associated  with  simultaneous  fading  of  radio  channels.  For  a  frequency* 
diversity  protection  switching  section  with  M  radio  channels  of  which 
u  are  protection  channels, 

M  -  u  +  N,  (5) 

eq.  (4)  becomes 

.V 

Z*  -  l  i T'(u  +  i),  (6) 

<-i 

where  T'(u  +  i)  is  the  accumulated  time  during  which  exactly  u  +  t  of 
the  M  radio  channels  have  failed  simultaneously.  A  failure  in  this 
context  is  an  exceedance  of  the  55-dBmcO  noise  value. 

The  simultaneous  fading  of  radio  channels  is  normally  described  in 
terms  of  the  amount  of  time  during  which  particular  sets  of  a  +  i 
channels  have  failed,  with  the  status  of  the  other  M  —  a  —  i  channels 
not  specified.  With  such  a  time  denoted  by  T>(u  +  t).  the  expression 
for  Zh  becomes  an  alternating  series 

.V  7<u+i> 

Zh  -  i  s  i-ir'csr*  t*(u  +  i),  (7) 

<•1  *•! 

where  J(u  +  t)  is  the  binomial  coefficient  C*  „  and  where  the  subscript 
k  enumerates  the  sets  of  u  +  t  channels  obtainable  from  the  M 
channels.  This  expression  has  been  used  by  other  authors,1  but  its 
general  mathematical  equivalence  to  the  definition  in  eq.  (6)  has  not 
been  demonstrated;  a  proof  of  this  equivalence  is  supplied  in  Appen¬ 
dix  A. 

The  expressions  for  the  simultaneous  failure  time  of  a  number  of 
radio  channels,  denoted  by  TAu  +  i),  are  based  on  a  generalized  model 
of  simultaneous  multipath  fading  in  an  arbitrary  number  of  radio 
channels  obtained  empirically  from  experimental  data.  The  estimation 
of  the  amount  of  fading  and  the  calculation  of  7*(u  + 1)  are  summarized 
in  Appendix  B.  A  simplified  form  of  T*(u  +  t),  without  path  length 
dependence  and  without  accommodation  of  fade  margin  differences 
between  radio  channels,  has  been  utilized  previously  to  demonstrate 
the  feasibility  of  joint  frequency-diversity  protection  for  radio  systems 
in  the  4-GHz  and  6-GHz  bands.1 

Numerical  evaluation  of  Zh  from  the  alternating  series  requires  a 
computer  program  that  efficiently  selects  and  classifies  radio  channel 
combinations  from  a  particular  frequency  plan  and  calculates  the 
corresponding  values  of  7*(u  +  t),  since  che  number  of  terms  can  be 
large.  For  example,  there  are  184,756  ten-channel  combinations  when 
all  radio  channels  in  the  4-GHz  and  6-GHz  bands  are  used  (400 A 
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protection  system  with  2  protection  and  18  regular  channels).  The 
computation  of  the  combinations  is  discussed  in  Appendix  C. 

III.  REPRESENTATIVE  RESULTS  FOR  FULLY  DEVELOPED  ROUTES 

A  fully  developed  4-GHz  route  contains  a  total  of  12  radio  channels 
consisting  of  one  protection  channel  and  11  regular  channels.  The 
notation  1  x  11  is  used  to  describe  this.  On  fully  developed  6-GHz 
routes,  the  total  number  of  radio  channels  is  8  and  the  protection 
scheme  is  1  x  7.  When  both  bands  are  utilized,  the  two  protection 
channels  can  be  placed  in  the  6-GHz  band,  which  results  in  2  x  18 
protection. 

Representative  estimated  values  of  service  failure  times  for  these 
protection  schemes  can  be  obtained  from  calculations  for  a  switching 
section  containing  one  25- mile  hop  with  average  terrain  and  climate. 
Without  diversity  protection,  the  service  failure  time  of  an  average 
channel  in  the  4-GHz  band  is  269  seconds  per  year  for  a  fade  margin 
of  37  dB.  In  the  6-GHz  band,  the  corresponding  service  failure  time  is 
208  seconds  per  year  for  a  fade  margin  of  49  dB. 

When  2  x  18  protection  is  applied,  the  calculated  service  failure  time 
is  6  seconds  per  year,  which  meets  the  long-haul  objective  of  10  seconds 
per  year.  Use  of  one  of  the  protection  channels  for  temporary  service 
in  the  fading  season  increases  the  service  failure  time  to  15  seconds 
per  year  as  a  result  of  the  1  x  19  configuration.  The  fade  margin  affects 
the  service  failure  time.  With  diversity  protection  applied,  a  5-dB  fade 
margin  decrease  in  all  channels  increases  the  service  failure  time  by  a 
factor  of  ten  (Fig.  1). 

For  1  x  11  and  1x7  protection,  the  calculated  service  failure  times 
are  about  15  seconds  per  year  (Fig.  1).  Prior  to  a  conversion  to  one 
protection  channel  per  band,  required  by  spectrum  conservation  mea¬ 
sures,  the  service  failure  times  for  2  x  10  and  2x6  protection  would 
have  been  about  6  seconds  per  year,  meeting  the  long-haul  objective. 

The  relative  positions  of  the  curves  in  Fig.  1  are  significant.  For 
example,  the  service  failure  time  of  an  average  channel  in  the  2  x  18 
scheme  is  between  those  for  average  channels  in  the  2  x  10  and  2x6 
schemes.  This  suggests  that,  because  of  decorrelation  of  fading  between 
the  4-GHz  and  6-GHz  bands,  the  2  X  18  functions  part  of  the  time  as 
a  2  x  12  and  part  of  the  time  as  a  2  x  6. 

The  frequency-diversity  parameters  that  can  be  used  to  scale  the 
results  to  other  path  lengths  and  climatic  conditions  are  discussed  in 
Appendix  O. 

IV.  PARTIALLY  DEVELOPED  ROUTES 

On  partially  developed  routes  where  only  a  part  of  all  available 
frequency  assignments  is  utilized,  the  service  failure  time  is  a  function 
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DEVIATION  F ROM  FADE  MARGIN,  ALL  CHANNELS.  IN  DECIBELS 

Fig.  1 — Service  failure  time  of  average  working  channel  due  to  multipath  fading  on  a 
25-miie  hop  with  average  climate  and  terrain.  Fade  margin  -  37  dB,  4-GHz  band:  fade 
margin  ”  40  dB.  6-4) Hz  band:  average  annual  temperature  •  55°F;  terrain  roughneea 

»  50  ft. 

’  Average  service  failure  time  and  service  failure  time  of  average  working  channel  are 
synonymous  in  the  context  of  this  work. 


of  the  frequency  arrangement  of  the  radio  channels.  A  “best-case” 
arrangement  is  one  where  the  frequency  separation  of  adjacent  chan¬ 
nels  is  maximized.  In  a  “worst-case”  arrangement,  the  channels  are 
crowded  into  the  high  end  of  the  frequency  band.  Such  arrangements 
do  not  necessarily  correspond  to  actual  or  permitted  growth  strategies, 
but  they  do  provide  bounds  for  the  variation  of  service  failure  time  on 
partially  developed  routes. 

In  continuation  of  the  example  in  Fig.  1,  calculated  results  for  the 
best  and  worst  cases  are  shown  in  Fig.  2  as  a  function  of  the  number 
of  working  channels.  The  best-case  results  for  the  4-GHz  band  show 
that  the  10  seconds-per-year  long-haul  objective  cannot  be  met  by 
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Fig.  2— Bounding  cases  for  service  failure  time  of  average  working  channels.  Fre¬ 
quency  assignments  are  listed  in  Appendix  E. 


frequency  diversity  alone  when  the  number  of  working  channels  is 
larger  than  seven.  The  ripple  in  the  worst  case  is  caused  by  the 
alternating  20-MHz  and  60- MHz  increments  as  channels  are  added. 

In  the  6-GHz  band,  the  long-haul  objective  cannot  be  met  by 
frequency  diversity  alone  when  the  number  of  working  channels  is 
larger  than  four.  The  6-GHz  worst  case  is  better  than  the  4-GHz  worst 
case  because  the  6-GHz  fade  margin  is  larger  than  the  4-GHz  fade 
margin  in  this  example. 

In  a  combined  protection  system  for  the  two  frequency  bands, 
introduced  after  the  4-GHz  band  is  fully  developed,  the  arrangement 
of  the  channels  in  the  6-GHz  band  has  only  a  small  effect  on  the 
service  failure  time  of  an  average  channel,  which  varies  between  5  and 
6  seconds  per  year. 

V.  EFFECT  OF  CHANNELS  WITH  REDUCED  FADE  MARGINS 

A  radio  channel  with  a  substandard  fade  margin  experiences  a  failure 
time  that  is  larger  than  normal.  The  presence  of  such  a  channel 
increases  protection  unavailability,  which  results  in  increased  service 
failure  times  for  all  working  channels.  The  increase  in  the  service 
failure  time  of  an  average  working  channel,  caused  by  the  presence  of 
one  or  two  substandard  radio  channels  on  a  fully  developed  route,  is 
shown  in  Fig.  3  as  a  function  of  the  reduction  of  the  fade  margin  from 
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Fig.  3 — Incnaa*  in  mnrica  failura  tim*  of  *v*nLg»  working  channai  whan  fada  marnna 
are  reduced  for  (1)  all  radio  channels;  (2)  two  midband  channela  in  1  x  7  system.  6  GHz; 
(3)  two  midhanri  channela  in  1  x  11  system.  4  GHs  <4>  one  midband  channel  in  1  x  7 
system,  6  GHz;  (5>  one  midband  channai  in  l  x  11  system.  4  GHz. 

its  normal  value.  As  an  example,  the  service  failure  time  of  an  average 
working  channel  in  a  1  x  11  system  is  doubled  when  the  fade  margin 
of  one  midband  radio  channel  is  degraded  by  7  dB;  this  is  equivalent 
to  a  1.5-dB  degradation  of  the  fade  margins  of  all  channels. 

The  service  failure  time  indicated  by  a  protection  switching  monitor 
actuated  by  the  protection  switching  equipment  can  be  larger  than  the 
actual  service  failure  time  if  adjustments  in  the  protection  switching 
system  are  incorrect  so  that,  for  a  particular  channel,  a  transfer  to 
protection  occurs  early,  at  a  noise  value  smaller  than  35  dBmcO.  The 
actual  service  failure  time,  while  smaller  than  that  indicated  (Fig.  4), 
is  larger  than  normal  because  occupancy  of  the  protection  channel  has 
been  increased. 

Transfers  to  protection  are  normally  permitted  only  when  the 
amount  of  noise  after  the  transfer  is  at  least  a  few  decibels  smaller 
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EARLY  TRANSFER  IN  DECIBELS 

Fig.  4—  increase  in  wvica  failure  cinra  of  average  working  channel  when  transfer  to 
protection  occurs  early  for  (1)  ail  channels;  (2)  one  mid  band  channel  in  1  x  7  system.  S 
GHz;  (3)  one  midband  channel  in  1  x  11  system,  4  GHz. 


chan  55  dBmcO.  The  increase  in  service  failure  time  due  to  a  3-dB 
difference  of  this  kind  is  very  roughly  the  same  as  that  due  to  transfers 
to  protection  that  occur  3  dB  early.  From  Fig.  4,  this  is  equivalent  to 
a  fraction  of  a  decibel  decrease  in  the  fade  margins  of  all  channels. 
Consequently,  3uch  differences  associated  with  the  protection  channels 
are  not  included  in  the  mathematical  model  presented  in  this  paper, 
since  their  relatively  small  effects  do  not  justify  the  added  analytical 
and  computational  complexity. 

VI.  SERVICE  FAILURE  TIME  OF  INDIVIDUAL  CHANNELS 

In  addition  to  the  service  failure  time  of  an  average  channel  utilized 
in  engineering  applications,  the  mathematical  modeif  developed  in  this 
work  can  provide  other  information  about  the  operation  of  frequency* 
diversity  protection.  In  particular,  the  service  failure  times  of  individ* 
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ual  working  channels  can  be  obtained  from  a  decomposition  of  the 
facility  service  failure  time. 

As  an  example,  consider  a  1  x  3  protection  system  in  the  4-GHz 
band  with  parameters  as  in  Fig.  2,  utilizing  channels  2,  4,  6  and  3  (80- 
MHz  separation  of  adjacent  channels).  The  simultaneous  failure  times 
for  particular  sets  of  exactly  u  +  i  channels  are  needed  in  the  decom¬ 
position.  These  times,  denoted  by  T*  (u  +  i)  and  listed  in  Table  I,  are 
summations  [as  specified  in  eq.  (14)]  of  the  failure  times  7*(u  +  i) 
calculated  as  outlined  in  Appendix  B.  The  numerical  values  illustrate 
the  frequency  selectivity  and  frequency  dependence  of  multipath  fad¬ 
ing.  For  example,  the  failure  time  of  channels  2  and  4,  but  not  6  and  8, 
which  can  be  caused  by  a  slope  in  the  frequency  response,  is  2.58 
seconds.  The  failure  time  of  channels  6  and  8,  but  not  2  and  4,  is  3.07 
seconds.  This  is  larger  than  the  2.58  seconds  because  the  amount  of 
fading  increases  with  frequency.  The  failure  time  of  channels  4  and  6, 
but  not  2  and  8,  is  2.22  seconds.  This  is  smaller  than  either  the  2.58 
seconds  or  the  3.07  seconds  because  this  failure  results  from  a  minimum 
in  the  frequency  response.  The  facility  service  failure  time  for  this 
example  is,  from  eq.  (6), 

6  4 

Z*  -  £  Ti(2)  +  2  £  T*(3)  +  3  TU4) 

l  1 

-  8.81  +  2  x  2.12  +  3  x  1.00 

■  16.05  channel-seconds  per  year.  (8) 

The  corresponding  service  failure  time  of  an  average  channel  is,  from 
eq.  (3),  one-third  of  this 

Ta  *  5.35  seconds  per  year.  (9) 

About  45  percent  of  this  time  occurs  when  3  or  all  4  of  the  radio 
channels  have  failed  simultaneously. 

The  expressions  for  the  service  failure  time  of  individual  channels 
are  obtained  by  inspection  from  Table  I.  In  the  case  where  channel  2 
is  used  for  protection,  the  service  failure  time  of  the  working  channel 
assigned  to  channel  4  is 

Ta.4  -  r,( 2)  +  <V4)  r«( 2)  +  (V4)  Tl(2) 

+  T'(3)  +  TU3)  +  (%>  T\( 3)  +  7”, (4) 

-  6.484  seconds  per  year,  (10) 

where  the  fractions  are  based  on  statistically  equal  protection  occu¬ 
pancy  in  cases  of  simultaneous  failures. 

The  results  of  such  calculations  (summarized  in  Table  II)  show  that 
the  location  of  the  protection  channel  affects  the  spread  of  the  individ- 
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Table  I — Simultaneous  failures  of 
exactly  u  +  i  radio  channels.  Channel 
numbers:  2,  4,  6,  3:  frequencies,  GHz: 


3.73,  3.81. 3.89,  3,97. 


u  +  * 

k 

Channels 

T/(u  +  <). 

sec/yeir 

2 

i 

2.4 

2.58 

2 

2 

2.6 

0.3? 

2 

3 

2,  8 

0.1? 

o 

4 

4,  6 

2.22 

2 

5 

4,8 

0.40 

2 

6 

6,8 

3.07 

3 

1 

2.  4.  6 

0.77 

3 

2.4.8 

0.23 

3 

3 

2,  6,  8 

0.23 

3 

4 

4.6,8 

0.89 

4 

i 

2.  4.  6,  3 

1.00 

Table  II— Service  failure  time  of  individual  working  channels  for  a 
1  x  3  system  in  the  4-GHz  band.  Channel  numbers:  2.  4,  6,  8; 
frequencies,  GHz:  3.73,  3.81,  3.89,  3.97. 


Working 

Channel* 

Service  Failure  Time,  Second*/ Year 
(percent  of  average  in  parentheses) 

Protection  channel  aaaignment 

2 

4 

6 

8 

2 

5.003 

3.898 

3.618 

(94) 

(73) 

(68) 

4 

6.484 

6.524 

5.434 

(121) 

(122) 

(101) 

8 

5.608 

6.753 

— 

6.998 

(105) 

(126) 

(131) 

8 

3.958 

4.294 

5.628 

— 

(74) 

(80) 

(105) 

Total  seconds 

16.05 

16.05 

16.05 

16.05 

Range  of  percent: 

Minimum 

74 

80 

73 

68 

Maximum 

121 

126 

122 

131 

•  Wortdnf  channel*  in  thii  cable  are  identified  by  the  number  of  the  corresponding 
regular  radio  channel. 


ual  service  failure  times.  The  facility  service  failure  time  and,  therefore, 
the  service  failure  time  of  an  average  channel  are  not  affected  by  the 
location  of  the  protection  channel.  The  largest  spread  of  the  individual 
service  failure  times  (68  to  131  percent  of  average)  occurs  when  the 
channel  at  the  highest  frequency  (channel  3)  is  used  for  protection. 

Protection  channel  assignments  could  be  optimized  based  on  results 
such  as  those  presented  in  Table  H  However,  a  general  treatment  of 
this  is  beyond  the  current  scope  of  this  study. 
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VII.  A00IT10N  OF  SPACE  DIVERSITY  TO  FREQUENCY  DIVERSITY 


Space-diversity  protective  switching  capitalizes  on  the  spati&i  selec¬ 
tivity  of  the  multipath  fading  process,  where,  at  an  instant  in  time  and 
at  a  given  radio  frequency,  the  strength  of  the  received  signal  is  a 
random  function  of  location  along  the  height  of  the  receiving  tower.* 7  9 
For  fm  radio  carrying  analog  message  traffic,  space-diversity  in  most 
cases  in  the  Bell  System  consists  of  two  vertically  separated  receiving 
antennas  and,  for  each  radio  channel,  a  waveguide  switch  that  connects 
a  receiver  to  either  of  the  two  antennas.  The  switch  is  activated  when 
the  receiver  input  decreases  below  a  threshold. 

The  presence  of  digital  traffic  in  the  fm  radio  network  affects  the 
choice  of  timing  in  the  switch.  The  number  of  antenna-to-antenna 
transfers  should  be  as  small  as  possible,  since  these  can  cause  errors  in 
digital  transmission.  For  thresholds  corresponding  to  fades  deeper 
than  about  35  dB,  this  is  accomplished  when,  in  the  case  of  simulta¬ 
neous  fading  on  the  two  antennas,  the  interval  between  repeated 
transfers  is  10  seconds.  Intervals  longer  than  this  cause  significant 
increases  in  the  service  failure  time  of  analog  message  circuits. 

When  added  to  frequency-diversity  to  reduce  service  failure  time, 
space-diversity  is  activated  before  frequency-diversity  to  avoid  unde¬ 
sirable  interactions  of  the  two  protection  systems.  Currently,  fre¬ 
quency-diversity  is  activated  at  a  noise  threshold  of  55  dBmcO,  and  52 
dBmcO  is  the  desired  threshold  value  for  space-diversity  activation 
when  it  is  added  to  frequency-diversity.  Such  an  arrangement  has 
assured  orderly  evolution  and  uniform  maintenance  procedures  of 
protection  switching  systems.  In  the  future,  because  of  increased  use 
of  space-diversity,  changed  threshold  values  may  be  desirable,  ideally 
58  dBmcO  for  frequency-diversity  and  55  dBmcO  for  space-diversity, 
to  further  reduce  service  impairments. 

The  reduced  service  failure  time  resulting  from  an  addition  of  space 
diversity  in  the  manner  described  above  can  be  calculated  from  eq.  (7) 
with  the  values  of  T*(u  +  t)  divided  by  the  appropriate  improvement 
value  for  threshold-switched  space-diversity.9  Such  calculations,  which 
can  be  readily  incorporated  in  a  computer  program,  will  not  be  pursued 
in  this  paper  where  the  emphasis  is  on  an  exploration  of  frequency- 
diversity  behavior. 
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APPENDIX  A 

Ooriv*Hon  at  trio  Attomattng  Soria* 

Generally,  we  prefer  to  describe  fading  in  terms  of  time  rather  than 
in  terms  of  probability,  since  the  latter  description  can  lead  to  inad- 
vertent  confusion  because  of  the  nonstadonarity  of  the  fading  process. 
However,  for  the  purposes  of  combinatorial  manipulation,  Tt>  can  be 
converted  to  a  fraction  designating  a  probability 

P*  -  n/T*  (11) 

where  2»  is  an  observation  period  (time  base)  such  as  a  year.  From 
eq&  (3)  and  (6), 

s 

Pa  -  IV"  £  iP'(u  +  1),  (12) 

where 

P'(u  +  0  -  T'{u  +  fl/T*.  (13) 

To  derive  the  alternating  series  from  eq.  (12),  we  apply  a  theorem 
from  probability  theory  (page  77  in  Ref.  14).  According  to  it,  the 
probability  that  exactly  u  +  t  out  of  M  channels  fail  simultaneously  is 

P'(«  +  0-  i  <-ir~CUS(r),  (14) 


where 

J(n 

Sir)  -  2  (15) 

*-i 

in  which  J (r)  is  the  binomial  coefficient  Cf,  and  where,  in  the  context 
of  this  work, 

Pkir )  -  T„ir)/T>.  (16) 

After  substitution  into  eq.  (12),  and  collection  of  like  terms. 

Pa  -  iv"  2  S  (-irppC^S(u  + »).  <i7) 

!•!  ^1 


Application  of  the  recursion  relationship 

C"  ■  +  2CJT*  +  Cf~‘ 

reveals  that,  because  of  cancellation  of  terms, 


2  (-i)^'pc^  -  err*. 


(18) 

(19) 
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Therefore 


P*  -  N~l  £  (-1  )-lCSt*Slu  +  t).  (20) 

(-1 

After  multiplication  by  NT*  this  becomes  the  alternating  series  expres¬ 
sion  for  Zh  in  eq  (7). 

APPENDIX  B 

Simultaneous  Failure  of  Radio  Channels 

In  the  case  of  a  single  radio  channel,  the  annual  amount  of  time 
during  which  the  received  signal  is  below  a  level,  as  a  result  of 
multipath  fading,  is*-9 

T  -  rTaL \  L  <  0J,  (21) 

where  L  describes  a  normalized  voltage  such  that  the  level  in  decibels 
relative  to  a  nonfaded  received  signal  is  20  log  L.  For  a  fade  margin  of 
F  dB,  where  F  is  a  positive  number,  the  corresponding  voltage  level  is 

U  -  KT"®.  (22) 

When  Lo  is  used  as  value  of  L  in  eq.  (21),  T  becomes  the  time  during 
which  the  channel  has  failed,  where  failure  is  defined  as  noise  in  excess 
of  55  dBmcO.  The  quantity  To  is  a  time  interval  related  to  the  length 
of  the  fading  season  [eq.  (20)  in  Ref.  9j.  The  fade  occurrence  factor  is5  9 

r  -  c(f/4)DnO~\  (23) 

where  D  is  the  path  length  in  miles,  f  is  the  carrier  frequency  in 
gigahertz,  and  c  describes  the  climate  and  terrain  (eqs.  (4)  and  (5)  in 
Ref.  9).  The  value  of  c  is  unity  when  the  climate  and  terrain  are 
average. 

Empirical  expressions  for  simultaneous  fading  of  radio  channels 
have  been  obtained  by  W.  T.  Barnett  from  experimental  data.  From 
these,  expressed  in  a  form  similar  to  that  for  a  single  channel,  the  time 
during  which  u  4-  i  radio  channels  have  failed  simultaneously  is 

Tk(u  4-  i)  ■  r*(u  +  DToLo,  (24) 

where  the  subscript  k  is  used  to  identify  various  sets  of  u  +  i  channels. 
This  expression  is  valid  when  T*<a  4- 1)  is  smaller  than  the  values  of  T 
for  the  channels  in  question.  The  simultaneous  failure  time  varies  with 
fade  margin  more  than  the  single  channel  failure  time  (Lo  versus  Ld 
variation).  The  quantity  Lo  describes  the  actual  fade  margin  when  it 
is  the  same  for  all  channels.  In  the  case  of  channels  with  differing  fade 
margins,  Lo  describes  a  nominal  fade  margin,  and  the  description  of 
differences  from  it  is  contained  in  r*(u  4-  i).  The  occurrence  factor  is 

r„(u  4-  i)  -  (cD*lQ~s/4QQ)fk(u  +  i ),  (25) 
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where 


Ilu+ii 

f„(u  +  0  -  (a  +  0/  £  (U/LxP)t(U/Lv,)t(SP/fi).  (26) 

The  quantities  in  this  expression  are  associated  with  pains  of  channels 
formed  from  the  u  +  i  channels.  The  number  of  such  pairs  is 

Hu  +  i)  -  cr*.  (27) 

The  subscript  p  identifies  the  pairs.  The  fade  margins  of  the  channels 
in  a  pair  are  described  by  Lip  and  Lip.  The  corresponding  carrier 
frequencies  fip  and  ftp  are  expressed  in  gigahertz  ( /ip  <  ftp).  The  average 
frequency  for  a  pair  is 


fp  -  (/ip  +  ftp)  12  (28) 

and  the  fractional  frequency  difference  is 

Sp  -  (ftp  -  fip)/fp.  (29) 

A  value  of  0.05  is  used  for  Sp  when  one  of  the  channels  is  in  the  4-GHz 
band  and  the  other  in  the  6-GHz  band. 

APPENDIX  C 

Computation  of  Combinations  of  Channeia 

The  evaluation  of  Z*  in  eq.  (7)  requires  identification  of  the  combi¬ 
nations  of  r  frequencies  that  can  be  formed  from  the  M  carrier 
frequencies  of  the  radio  channels.  If  the  frequencies  are  tagged  using 
the  integers  I  through  M,  then  the  combinations  required  are  subsets 
of  the  set  (1,  2,  •  •  • ,  M] .  The  integers  appearing  in  the  subsets  serve 
as  pointers  (used  pairwise)  to  precalculated  values  of  5 p/fl  for  fre¬ 
quency  pairs  in  a  subset  [eqs.  (28)  and  (29)]. 

Given  M  and  r,  the  computation  is  initialized  by  defining  (for  i  »  1, 
2.  ...,r) 

a,-  «•  i,  i  r  (30) 

<Xi  -  (r  -  1),  i  m  r  (31) 

fa  -  M  -  r  +  i.  (32) 

The  subsets  are  then  generated  as  follows: 

U)  Starting  with  j  -  r.  find  the  first  j  such  that  aj<  P,. 

(it)  Replace  a,  with  a,  +  1. 

(iii)  For  k  •  (j  +  1),  (j  +  2),  •  •  •,  r  replace  a*  with  a*-i  +  1. 

In  this  procedure,  elements  of  subsets  are  incremented  to  form  new 
subsets,  and  the  elements  of  each  subset  are  ordered  within  the  vector 
a  such  that  or,  <  a,*i.  Repeated  generation  of  the  same  combination  of 
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frequencies  is  thereby  avoided.  The  first  iteration  after  initialization 
produces  the  subset  {1,  2,  •••,  r}.  The  iteration  is  continued  until 
execution  of  step  (i)  fails,  when  all  the  required  C*  combinations  have 
been  produced.  The  potentially  large  number  of  combinations  neces¬ 
sitates  double  precision  computation  in  summations  over  the  subsets. 

APPENDIX  D 

Frequency-Diversity  Parameters 

The  service  failure  time  of  an  average  channel  can  be  expressed  as 

T, i  -  Tua/Ih,  (33) 

where  Tua  is  the  unprotected  service  failure  time  for  the  hop  in  question 
obtained  from  eq.  (21)  after  substitution  of  Lo  from  eq.  (22).  The 
improvement  resulting  from  the  use  of  frequency-diversity  protection 
is 

h  -  qU\  (34) 

where 

<7-100  fo/DG,  (35) 

in  which  D  is  the  path  length  in  miles  and  fo  is  the  frequency  in 
gigahertz  at  which  Tua  is  calculated.  An  average  of  the  carrier  frequen¬ 
cies  of  the  radio  channels  is  a  suitable  value  for  the  reference  frequency 
fo.  Alternatively,  f0  can  be  the  center  frequency  of  a  band,  which 
facilitates  band  identification  when  the  calculations  are  performed  for 
frequency  plans  formulated  by  the  International  Radio  Consultative 
Committee  (ccir).  The  parameter  G  is  determined  by  the  carrier 
frequencies  of  the  radio  channels  and  by  the  channel-to-channel  vari¬ 
ations  of  fade  margins 

W  Jiu+t) 

g  - iV_l  I  l  (-i)‘"lcrrV*(u  +  0,  (36) 

<  —  i  **  i 

where  fk(u  +  i)  and  the  summation  are  defined  in  Appendix  B.  In  the 
simple  case  of  1  x  1  frequency  diversity  with  the  same  fade  margin  for 
both  channels, 

G  -  2  fi/bf,  (37) 

where  fp  is  the  average  and  A  f  is  the  difference  of  the  two  carrier 
frequencies,  expressed  in  GHz.  The  corresponding  value  of  q  in  this  1 
x  1  case, 

q  -  (fo/fp)  (5 Q/fpD)  (Af/fP)  (38) 

has  been  used  previously  with  the  factor  fo/fP  approximated  by  unity.9 
When  the  fade  margins  of  all  radio  channels  are  identical,  the 
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Table  III — Values  of  G  for  fully 
developed  routes.  Fade  margins  are 
identical  for  all  radio  channels  in  a 


given  protection  scheme. 


Frequency 

Band, 

GHz 

Protection 

Scheme 

Value 
of  G 

4 

2  x  10 

1597 

4 

1  X  11 

4682 

6 

2x6 

7380 

6 

1  x  7 

17059 

Frequencies  in  the  4-GHz  band: 

3.670  +  0.040  i,  t  -  1. 3,  S.  .  .  11 

3.660  +  0.040  i.  i  -  2.  4.  6.  ...  12 

Frequencies  in  the  6-GH2  band: 

3.9452,  5.9748,  6.0045,  6.0342, 

6.0638,  6.0935,  6.1231,  6.1528 


parameter  G  depends  only  on  the  frequency  arrangement  of  the 
channels.  Tabulated  values  of  G,  such  as  shown  in  Table  IQ  for  fully 
developed  routes,  can  therefore  be  used  in  eq.  (33)  to  calculate  the 
service  failure  time  for  common  frequency  arrangements. 

In  the  2  x  18  case  in  Fig.  1,  the  choice  of  reference  fade  margin  (—20 
log  La)  becomes  arbitrary  because  there  are  two  fade  margins  (37  and 
40  dB,  in  the  4-  and  6-GHz  bands,  respectively).  The  choice  has  no 
effect  on  Th,  but  it  does  affect  the  various  factors  in  Th  such  as  G. 
When  Tua  and  fa  are  the  averages  for  the  20  radio  channels,  the 
reference  fade  margin  is  38.27  dB  [obtained  from  eq.  (21)].  The  value 
of  G  is  3129  for  this  choice  of  reference  values. 

APPENDIX  E 

Frequency  Assignments  in  Figure  2 

To  simplify  presentation,  the  radio  channels  in  this  work  are  num¬ 
bered  1  to  20  in  order  of  increasing  frequency.  The  frequencies  range 
from  3.71  to  6.1528  GHz  as  listed  in  Table  IQ. 

In  the  4-GHz  worst  case,  channels  12  and  11  are  utilized  in  the  1  x 
1  protection  scheme,  and  the  sequence  in  which  additional  working 
channels  are  assigned  is  10,  9,  8,  •  •  • ,  1.  In  the  best  case,  channels  2  and 
12  (same  polarization)  are  used  for  the  1  x  1  scheme,  and  then  channel 
6  is  added  to  form  a  1  x  2  scheme.  Channels  2,  4,  8,  and  12  are  utilized 
in  the  1x3  scheme,  and  the  sequence  in  which  additional  working 
channels  are  assigned  is  6,  10,  1,  3,  5,  7,  9,  and  11. 

In  the  6-GHz  worst  case,  channels  20  and  19  are  utilized,  in  the  1  x 
1  scheme,  and  the  sequence  in  which  additional  working  channels  are 
assigned  is  18,  17,  •  •  • ,  13.  In  the  best  case,  the  channel  utilization  is 
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m 


I 

X 

1  - 

13, 

20 

I 

X 

2  - 

13, 

16, 

20 

I 

X 

3  - 

13, 

15, 

17, 

20 

1 

X 

4  - 

13, 

14, 

16, 

18, 

20 

1 

X 

5  - 

13, 

14, 

15, 

17, 

19, 

20 

1 

X 

6  - 

13, 

14, 

15, 

16, 

18, 

19, 

20 

1 

X 

7  - 

13, 

14, 

15, 

16, 

17, 

18, 

19,20 

In  the  combined  protection  system  for  the  4-GHz  and  6-GHz  bands, 
all  4-GHz  channels  and  channels  20  and  19  are  used  in  the  worst-case 
2  x  12  scheme.  The  sequence  in  which  additional  working  channels 
are  assigned  is  13, 17,  16,  •  -  •,  13.  In  the  best-case  2  x  12  scheme,  all  4- 
GHz  channels  are  also  used,  but  the  6-GHz  channels  are  13  and  20. 
The  pattern  for  the  utilization  of  additional  channels  is  the  same  as 
that  in  the  6-GHz  best  case. 
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ABSTRACT 

Noise  performance  of  average  as  well  as  individual  working  chan¬ 
nels  in  microwave  1 ine-of-sight  FM  radio  systems  can  be  estimated 
from  a  model  of  multipath  fading  and  frequency-diversity  opera¬ 
tion  utilized  in  the  Bell  System.  The  performance  degradation 
imposed  by  radio  channels  with  substandard  fade  margins  can  also 
be  estimated  from  the  model. 

INTRODUCTION 

Microwave  FM  radio  systems  utilized  in  the  Bell  System  pro¬ 
vide  protected  broadband  communications  channels  that  generally 
have  a  capacity,  depending  or.  the  frequency  band  and  on  the  type 
of  equipment,  of  1200  to  2400  message  circuits  per  channel.  The 
protected  channels,  with  protection  provided  by  radio  switching 
systems,  ace  referred  to  as  "working"  channels,  to  distinguish 
these  from  the  actual  radio  channels,  which  are  assigned  as 
either  "regular”  or  "protection"  channels. 

The  communications  traffic  in  a  working  channel  can  be  tem¬ 
porarily  subjected  to  high  noise  when  the  number  of  impaired 
radio  channels  exceeds  the  number  of  protection  channels.  Bell 
System  practice  is  to  describe  this  effect  in  terms  of  the  amount 
of  time  during  which  the  noise  in  a  working  channel  exceeds  55 
d3rnc0  in  a  message  circuit  at  the  top  frequency  assignment  in 
the  baseband1.  In  this  paper,  this  time  is  referred  to  as  "ser¬ 
vice  failure  time"". 

In  modern  microwave  radio  systems,  protection  requirements 
are  frequently  governed  by  the  presence  of  multipath  fading-'"13, 
as  opposed  to  causes  associated  with  equipment  or  human  interven¬ 
tion*.  Estimates  of  service  failure  time  due  to  multipath  fading 
are  therefore  needed  to  determine  the  adequacy  of  diversity  pro¬ 
tection  arrangements  when  new  routes  are  planned,  when  transmis¬ 
sion  parameters  of  existing  routes  are  changed,  or  when  alternate 
protection  strategies  are  considered  for  future  use. 

A  mathematical  model14  utilized  in  the  Bell  System  to 
obtain  such  estimates  for  frequency-diversity  protected  FM  radio 
systems  is  summarized  and  illustrated  in  this  paper.  The 

CH1514- 9/79/0000-0072:00. 75  £.1979  IEEE 
7-20 


improvement  resulting  from  additions  of  space-diversity  protec¬ 
tion^  to  frequency-diversity  protection  can  be  incorporated  in 
th i s  mod  e 1 . 


PSSFORM^NCS  0?  .V/5RAGS  CHANNELS 

In  FM  radio  systems,  the  dominant  form  of 
multipath  fading  is  thermal  noise,  which  is  pro 
depth.  The  annual  amount  of  time  during  which 
of  noise  is  exceeded  can  therefore  be  described 
amount  of  time  during  which,  as  a  result  of  mul 
received  signals  an  a  radio  hop  are  below  parti 
performance  of  a  frequency-diversity  switching 
a  number  of  hops  is  otained  by  addition  of  the 
values  for  the  various  hops,  since  the  probabil 
aus  deep  fading  on  different  hops  is  much  small 
bility  of  simultaneous  fading  in  radio  channels 


noise  during  deep 
portional  to  fade 
a  particular  value 
in-iterms  of  the 
tipath  fading,  the 
cuiar  levels.  The 
section  containing 
time- be low-level 
ity  of  simultane- 
er  than  the  proba- 
on  the  same  hop. 


In  the  case  of  a  single  unprotected  radio  channel,  the 
jated  annual  tirae-below-Ievel  on  a  particular  radio  hop 


Tu  3  fT0L< 


L  <  0.1 


where  L  describes  a  normalized  voltage,  such  that  the  level  in  dB 
relative  to  a  nonfaded  received  signal  is  20  log  L.  The  time  T^ 
is„related  to  the  length  of  the  fading  season.  Its  value  is  B  x 
10’  seconds  for  an  average  geographic  location  in  the  United 
States.  The  fade  occurrence  fact  is5,9: 


c( f/4)D3l0’5 


where  0  is  the  path  length  in  miles  and  f  is  the  carrier  fre¬ 
quency  in  GHz.  The  value  of  c  is  unity  when  the  climate  and  ter¬ 
rain  are  average.  For  a  fade  margin  of  F  dB  (defined  by  noise 
of  55  dBrncO)  where  F  is  a  positive  number,  the  correspond ing 
voltage  level  is: 


r-*  v  w  ■  m 


WW'- 


io-(P/,0) 


'Vhen  is  asad  as  a  value  of  L  in  aquation  (1),  T  becomes  the 

time  during  which  the  channel  has  failed,  where  failure  is  defind 
as  noise  in  excess  of  55  dBrncO  at  the  top  frequency  assignment 
in  baseband. 

The  equivalent  time-beiow-level  of  an  average  working  chan¬ 
nel  in  a  frequency-diversity  protection  system  (Figure  11  can  be 
expressed  as: 


(DG/l00fo)  rT0L^' 


t  <  r 


where  f-  is  the  frequency,  expressed  in  GHz,  at  which  r  is  c 
lated.  An  average  of  the  carrier  frequencies  of  the  radio  c 
nels  is  a  suitable  value  for  the  reference  frequency  f«.  Th 
parameter  G  is  determined  by  the  carrier  frequencies  or  the 
channels  and  by  the  channel-to-channel  departures  of  fade  ma 

from  wh 


from  the  reference  value  Lg .  The  mathematical  model  from  wh 
equation  (4)  is  obtained,  and  from  which  G  can  be  calculated 
described  in  Appendix  A.  The  expression  for  T  describes  the 
vice  failure  time  occurring  in  switched  frequency-diversity 
terns  in  which  the  values  of  noise  defining  the  fade  margins 
the  initiation  of  switching  are  identical. 


alcu- 

han- 

e 

radio 
rgins 
ich 
,  is 
ser- 
sys- 
and 


Representative  estimated  values  o 
be  otained  from  calculations  for  a  swi 
one  25-mile  hop  with  average  terrain  a 
for  a  fully  developed  4-GHz  route  (one 
working  channels) ,  the  calculated  serv 
average  working  channel  is  approximate 
channels  have  a  fade  margin  of  37  dB  ( 
failure  time  changes  by  a  factor  of  te 
all  channels  change  5  dB.  Similarly, 
fade  margins  increases  the  service  fai 
50  percent. 


f  service  failure  time  can 
tching  section  containing 
nd  climate.  In  this  case, 
protection  channel,  11 
ice  failure  time  of  an 
ly  15  seconds/year  when  all 
Figure  2) .  The  service 
n  when  the  fade  margins  of 
a  one-dB  decrease  in  the 
lure  time  by  approximately 


Transmission  performance  objectives,  as  currently  applied  t 
the  engineering  of  microwave  FM  radio  routes9'1,1  allocate  one- 
quarter  of  a  0.02  percent  annual  objective  to  service  failures 
caused  by  multipath  fading  in  a  one-way  circuit  on  a  4000-mile 
long-haul  or  250-raile  short-haul  route.  This  allocaton  is 
prorated  to  the  actual  length  of  the  switching  section.  In  the 
case  of  a  25-mile  switching  section  in  long-haul  service,  the 


service  failure  time  should  be  less  than  10  seconds/year.  In  the 
example  in  the  previous  paragraph,  where  the  estimated  service 
failure  time  was  15  seconds/year,  larger  fade  margins  or  addition 
of  space-diversity  protection  are  needed  to  meet  design  objec¬ 
tives. 


when  the  fade  margins  of  all  radio  channels  are  identical, 
the  parameter  G  depends  only  on  the  frequency  arrangement  of  the 
channels.  Tabulated  values  of  G,  such  as  shown  in  Table  1  for 
fully  developed  4,-GSz  and  6-GHz  routes,  can  be  used  in  such  cases 
to  simplify  calculations  for  standard  frequency  arrangements. 

The  value  of  G  is  4632  in  the  1  x  111  example  discussed  above 

On  partially  developed  routes  where  only  a  part  of  all 
available  frequency  assignments  is  utilized,  the  service  failure 
tine  is  a  function  of  the  frequency  arrangement  of  the  radio 
channels.  A  “best  case"  arrangement  is  one  where  the  frequency 
separation  of  adjacent  channels  is  maximized.  In  a  "worst  case" 
arrangement,  the  channels  are  crowded  into  the  high  end  of  the 
frequency  band.  Such  arrangements  do  not  necessarily  correspond 
to  actual  or  permitted  growth  strategies,  but  they  do  provide 
bounds  for  the  variaton  of  service  failure  time  on  partially 
developed  routes. 

In  continuation  of  the  axaraple  in  Figure  2,  calculated 
results  for  the  best  and  worst  cases  are  shown  in  Figure  3  as  a 
functon  of  the  number  of  working  channels.  The  best  case  results 
for  the  4-G3z  band  show  that  the  10  sec/year  long-haul  objective 
cannot  be  met  by  frequency  diversity  alone  when  the  number  of 
working  channels  is  larger  than  seven.  The  ripple  in  the  worst 
case  is  caused  by  the  alternating  '20  MHz  and  60  MHz  increments  as 
channels  are  added1,4. 

In  the  6 -GHz  band,  the  long-haul  objective  cannot  be  met  by 
frequency  diversity  alone  when  the  number  of  working  channels  is 
larger  than  four.  The  6-GHz  worst  case  is  better  than  the  4-GHz 
worst  case  because  the  6-G3z  fade  margin  (40  d3)  is  larger  than 
the  4-GHz  fade  margin  (37  dBI  in  this  example. 

In  a  combined'"protaction  system  for  the  two  frequency  bands, 
introduced  after  the  4-GHz  band  is  fully  developed,  the  arrange¬ 
ment  of  the  channels  in  the  6-ghz  band  has  only  a  small  effect  on 
the  service  failure  time  of  an  average  working  channel,  which 
varies  between  5  and  5  sec/year. 

PFHFOHMAMCS  OF  INDIVIDUAL  CHANNELS 

In  addition  to  the  service  failure  time  of  an  average  work¬ 
ing  channel  utilized  in  route  design,  the  mathematical  model 
presented  in  the  work  can  provide  other  information  about  the 
operation  of  frequency-diversity  protection.  In  particular,  the 
service  failure  times  of  individual  working  channels  can  be 
obtained  from  a  summation  of  the  times  during  which  particular 
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sacs  of  radio  channels  have  failed. 

As  an  examole,  consider  a  1  x  3  proCecttion  system  in  the  4- 
GHz  band  with  parameters  as  in  Figure  3,  utilizing  channels  2,  4, 
5  and  3  (30  hHz  separation  of  adjacent  channels)  .  The  simultane¬ 
ous  failure  times  for  particular  sets  of  exactly  u  +  i  channels, 
denoted  by  ?V(u+i)  and  listed  in  Table  2,  have  been  calculated 
as  outlined  in  Appendix  A.  The  numerical  values  illustrate  the 
frequency  selectivity  and  frequency  dependence  of  multipath  fad¬ 
ing."  For  example,  the  failure  time  of  channels  2  and  4,  but  not 
3  and  3,  which  can  be  caused  by  a  slope  in  the  frequency 
rasconse,  is  2.53  seconds.  The  failure  time  of  channels  5  and  3, 
but'not  2  and  4,  is  3.07  seconds.  This  is  larger  than  the  2.53 
seconds  because  the  amount  of  fading  increases  with  frequency. 

The  failure  time  of  channels  4  and1  5,  but  not  2  and  3,  is  2.22 
seconds.  This  is  smaller  than  either  the  2.53  seconds  or  the 
3.07  seconds  because  this  failure  results  from  a  minimum  in  the 
frequency  response,  which  has  a  smaller  probability  of  occurrence 
than  a  slope  in  the  frequancy  response. 

The  service  failure  time  of  an  averaqe  working  channel  in 
this  example  is,  from  the  definition  in  equation  (7)  in  Appendix 


T 


(1/3) 


5 

(.3  T^(2) 


4 

+  25  ?^(3)  +  3  T’(4)) 

.  k«l 


(5) 


»  (1/3)  (3.31  +  2  x  2.12  +  3  :<  1.00) 


*  5.35  seconds/year 


About  45  percent  of  this  time  occurs  when  3  or  all  4  of  the 
radio  channels  have  failed  simultaneously.  This  example  demon¬ 
strates  that  the  contribution  of  higher  order  simultaneities  to 
the  service'  'failure  time  cannot  be  neglected. 

The  expressions  for  the  service  failure  time  of  individual 
channels  are  obtained  by  inspection  from  Table  2.  In  the  case 
where  channel  2  is  used  for  protection,  the  service  failure  time 
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o€  the  working  channel  assigned  to  channel  4  is: 


Tch4  =  r[(2)  *  (1/2)  T4(25  +  (1/2)  Tg(2) 


(<5) 


+T  ^  ( 3 )  +  T^(3)  +  (2/3)  T4(3)  -*-t’(4) 


*  5.434  seconds/year 


where  the  fractions  are  based  on  statistically  equal  protection 
occupancy  in  cases  of  simultaneous  failures. 

The  results  of  such  calculations  (summarized  in  Table  3) 
show  that  the  location  of  the  protection  channel  affects  the 
spread  of  the  individual  service  failure  times.  The  service 
failure  time  of  an  average  channel  is  not  affected  by  the  loca¬ 
tion  of  the  protection  channel.  The  largest  spread  of  the  inui- 
vidual  service  failure  times  (S3  to  131  percent  of  average1) 
occurs  when  che  channel  at  the  highest  frequency  (channel  3)  is 
used  for  protection  (Figure  4)- 

Protection  channel  assignments  could  be  optimized  based  on 
results  such  as  presented  in  Table  3.  However,  a  general  treat¬ 
ment  of  this  is  beyond  the  current  scope  of  this  study. 

THE  EFFECT  OF  SaSST4N0AR0  CHANNELS 

A  radio  channel  with  a  substandard  fade  margin  experiences  a 
failure  time  that  Ts  larger  than  normal.  The  presence  of  such  a 
channel  increases  protection  unavailability,  which  results  in 
increased  service  failure  times  for  all  working  channels.  The 
increase -in  service  failure  time  of  an  average  working  channel, 
caused  by  the  presence  of  one  or  two  substandard  radio  channels 
an  a  fully  developed  route,  is  shown  in  Figure  5  as  a  function  of 
the  reduction  of  the  fade  margin  from  its  normal  value.  4s  an 
example,  the  service  failure  time  of  an  average  working  channel 
in  a  1  x  11  system  is  doubled  when  the  fade  margin  of  one  midband 
radio  channel  is  degraded  by  7  dB;  this  is  equivalent  to  a  1.5  d3 
deqradaton  of  the  fade  margins  of  all  channels. 
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APPENDIX  A 


THE  MATHEMATICAL  MODEL 


The  service  failure  time  of  an  average  working  channel  is 
defined  as 


T  =»  N  Z  iT'  (u+i)  (7) 

i=L 

where  N  is  the  number  of  working  channels  and  u  is  the 
number  of  protection  channels.  The  total  number  of  radio 
channels  is 


M  a  U  +  N  (8) 

The  quantity  T'  (u+i)  is  the  accumulated  time  during  which 
exactly  u+i  of  the  M  radio  channels  have  failed  simultaneously 

The  simultaneous  fading  of  radio  channels  is  normally 
described  in  terms  of  the  amount  of  time  during  which  parti¬ 
cular  sets  of  u+i  channels  have  failed,  with  the  status  of 
the  other  M-u-i  channels  not  specified.  With  such  a  time 
denoted  by  T^Cu+i),  the  expression  for  T  becomes  an  alternat¬ 
ing  series^/  14 


-i  N  JCu+i)  .  . 

T  *  N  Z  Z  (-11  1  C  ,  *  T.  Cu+i)  (9) 

i-1  k»l  u”i  K 

M 

Where  J(u+i)  is  the  binomial  coefficient  ,  and  where  the 

subscript  k  enumerates  the  sets  of  u+i  channels  obtainable 
from  the  M  channels. 

Empirical  expressions  for  simultaneous  fading  of  radio 
channels  have  been  obtained  by  W.  T.  Barnett  from  experimental 
data: 

T^  (u+i)  *  (C  D410”5/400)  fk(u+i)  TQL4  (10) 

wfiera 

X  i)  j  ^  9 

Cu+i)  -  Cu+i)/  £  CL0/L1p)^CL0/L2p)n5p/f‘)  (U) 

P** 


The  quantities  in  this  expression  axe  associated  with  pairs 
of  channels  formed  from  the  u+i  channels .  The  number  of 
such  pairs  is 

ICu+i)  -  (12) 

The  subscript  p* identifies  the  pairs.  The  fade  margins  of 
the  channels  in  a  pair  are  described  by  Lip  and  L2p.  The 

corresponding  carrier  frequencies  fip  and  f20  are  expressed 
in  GHz  Cfip<f2p) •  The  average  frequency  for"  a  pair  is 

f  *  Cf.  +  f-)/2  (13) 

P  IP  2p 

and  the  fractional  frequency  difference  is 

5p  *  «2p  -  V/fp  tl4) 

A  value  of  0.05  is  used  for  S0  when  one  of  the  channels  is 
in  the  4-GHz  band  and  the  other  in  the  6  GHz  band. 


The  Siam 


,  N  J (u+i) 

G  =>  Z  Z 
i=l  k=l 


l-l)i-1  C^i’2  fk(u+i) 


contains  the  effects  of  the  frequency  arrangement  of  the 
radio  channels  and  the  deviations  of  the  fade  margins  from 
a  reference  value.  A  simplified  form  of  T^(u+i)  has  been 
used  previously  to  demonstrate  the  feasibility  of  joint 
frequency-diversity  protection  for  radio  systems  in  the 
4  GHz  and  6  GHz  bands . ^ 

The  alternating  series  in  equation  (9)  can  be  derived 
formally  from  the  definition  in  equation  C7)  when  a  theorem 
from  probability  theory^S  is  used  to  relate  the  primed  and 
unprimed  quantities^4 


T* (u+i)  »  Z 


r»u-*-i 


Cj+i  S(r) 


where 


JCr) 

Z 

k»l 


Tk(r) 


in  which  JCr)  is  the  binomial  coefficient 


7-30 


l  AO-A087  793 

defense  communications 

PROCEEDINGS  OF  SEMINAR 
jAN  80  D  R  SMITH 

ENGINEERING  CENTER  RESTON  VA 

ON  frequency  selective  fading  AND 

- 

F/G  17/2.1 

ITS  EF— ETCIU) 

UNCLASSIFIED 

DCEC-TN-2-80 

SBIE-AD-EIOO  375 

NL 

The  quantity  T'  Cu+i)  can  be  decomposed  into  a  sum  of 
terms,  denoted  by  T£(u+i),  describing  the  failure  of 

particular  sets  of  exactly  u+i  radio  channels.  This  and 
equation  (16)  provide  sets  of  equations  from  which  values 
of  T£(u+i)  can  be  calculated. 


VALUES  OF  G  FOR  FULLY  DEVELOPED  ROUTES 

FAOE  MARGINS  ARE  IDENTICAL  FOR  ALL  RAOIO  CHANNELS  IN  A 
GIVEN  PROTECTION  SCHEME 


FREQUENCY 
BAND,  GHz 

PROTECTION 

SCHEME 

VALUE 

OF  G 

4 

2X10 

1597 

4 

1X11 

4682 

6 

2X  6 

7380 

6 

IX  7 

17059 

FREQUENCIES  IN  THE  4-GHz  SAND: 
3.670  +0.0401 ,  I  si,  3,  5,...  11 
3.650 +0.0  401,  I  s  4,  6, . . .  1 2 

FREQUENCIES  IN  THE  6  GHz  BAND: 
5.9452,  5.9748,  6.0045,  6.0342, 
6.0638,  6.0935,  6.1231,  6.1528 


TABLE  1 
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SIMULTANEOUS  FAILURES  OF  EXACTLY  u+t  RAOIO  CHANNELS 
CHANNEL  NUMBERS:  2,  4,  6,  8. 

FREQUENCIES,  <3Hz:  3.73,  3.81,  3.89,  3.97 


u  + 1 

k 

CHANNELS 

T^(u-H),  SEC/YEAR 

2 

1 

2,  4 

2.58 

2 

2 

2,  6 

0.37 

2 

3 

2,  a 

0.17 

2 

4 

4,8 

2.  22 

2 

3 

4,  8 

0.  40 

2 

6 

6,  8 

3.07 

3 

-1 

2,  4,  8 

0.77 

3 

2 

2,  4,  8 

0.  23 

3 

3 

2,6,  8 

0.23 

3 

4 

4>  8,  8 

0.  89 

4 

1 

2,  4,  8,  8 

1.  00 

TABLE  Z 
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SERVICE  FAILURE  TIME  OF  INDIVIDUAL  WORKING  CHANNELS 
FOR  A  1x3  SYSTEM  IN  THE  4-GHz  BAND 


CHANNEL  NUMBERS:  2,  4,  6,  8 
FREQUENCIES,  GHz:  3.73,  3.81,  3.89,  3.97 


WORKING 

CHANNEL* 

SERVICE  FAILURE  TIME, 

SECONDS/YEAR 

PERCENT  OF  AVERAGE  IN  BRACKETS 

PROTECTION  CHANNEL  ASSIGNMENT 

2 

4 

6 

8 

2 

— 

5.003 

(94) 

3.898 

(73) 

3.618 

(68) 

4 

6.484 

(121) 

— 

6.524 

(122) 

5.434 

(101) 

6 

5.608 

(105) 

6.753 

(126) 

— 

6.998 

(131) 

8 

3.958 

(74) 

4.294 

(80) 

5.628 

(105) 

— 

TOTAL 

SECONOS 

16.05 

IB 

16.05 

16.05 

RANGE  OF 

percent: 

MINIMUM 

MAXIMUM 

74 

121 

80 

126 

73 

122 

68 

131 

*  WORKING  CHANNELS  IN  THIS  TABLE  ARE  IDENTIFIED  BY  THE 
NUMBER  OF  THE  CORRESPONDING  REGULAR  RADIO  CHANNEL. 


TABLE  3 
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FIGURE  1  FREQUENCY  DIVERSITY  PROTECTION 
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BELOW  LEVEL  SECONDS  PER  YEAR 


AVERAGE  TERRAIN  AND  CLIMATE 
25-MILE  PATH 
f0*  3.92  GHz 


■AVERAGE  UNPROTECTED 
RADIO  CHANNEL 


IMPROVEMENT  AVAILABLE 
AT  A  PARTICULAR  LEVEL" 


AVERAGE  WORKING- 
CHANNEL 


LEVEL  RELATIVE  TO  NORMAL,  dB 


FIGURE  2  TIME  BELOW  LEVEL  FOR  A  FULLY  DEVELOPED  (Ix-H)  4-GHi 
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SERVICE  FAILURE  TIME  OF 
AVERAGE  WORKING  CHANNEL,  SECONDS.  /EAR 


>-.L- 


AVERAGE  TERRAIN  AND  CLIMATE 
25 ‘MILE  PATH 

FADE  MARGIN *37dB,  4-GHz  BAND 
*40dB,  6 -GHz  BAND 

4-GHz  WORST  CASE 


5-GHz  WORST 
CASE 


4-GHz  BEST  CASE 


5-GHz 

BEST 


COMBINED  4-GHz  AND  6-GHz 
WORST  CASE 


BEST  CASE 
PROTECTION  CHANNELS 


2  3 


5  6  7  8  9  -10  1-1  i2  -13  -14  15  16  17  18 

NUMBER  OF  WORKING  CHANNELS 


*'GURE  3  80UNDING  CASES  FOR  SERVICE  FAILURE  TIME  OF  AVERAGE 
WORKING  CHANNELS 
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TIME  BELOW  LEVEL,  SECONDS  PER  YEA 


FIGURE  4  EXAMPLE  OF  TIME  BELOW  LEVEL  FOR  INDIVIDUAL  WORKING 
CHANNELS 
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FACTOR  BY  WHICH  AVERAGE  SERVICE  FAILURE  TIME  IS  INCREASED 


8.  THE  EFFECT  OF  SELECTIVE  FADING  ON  DIGITAL  RADIO 


Steve  Barber 

Bell -Northern  Research 

Ottawa,  Ontario,  Canada 


NOTE:  The  materials  provided  in  this  section  are  copies 

of  an  article  from  the  1978  International  Conference 
on  Communications  used  as  the  basis  for  the  extempor¬ 
aneous  remarks  given  by  Mr.  Barber. 


8-1 


THE  EFFECT  Of  SELECTIVE  FACINC  ON  DIGITAL  RADIO 


C.W.  Anderson  S.  Barber  R.  Pacel 

3ell-Morchern  Research,  Ottawa,  Canada 


Abstract 

In  the  coursa  of  developnenc  of  the  DR5-3 
Digital  Radio  System  a  sarlas  of  digital,  propagation 
trials  ware  conducted  over  a  51  km  path  ac  3  GHz  Co 
determine  the  effect  of  multipath  propagation  on  a  40 
MHz  bandwidth  digital  radio  system.  It  was  found  that 
the  affects  of  frequency  selective  fading  result  In 
unaccepcably  high  system  unavailability  unless  adap¬ 
tive  equalization  and  space  diversity  are  employed. 
Specifically  che  results  are: 

L.  Multipath  induced  outage  is  much  higher  chan 
would  be  predicted  from  che  measured  flac  fade  margin 
of  the  equipment.  For  a  non  diversity  system,  the 
probability  of  outage  for  a  single  hop  in  the  worst 
fading  month  is  approximately  1  x  10"'.  This  is  700 
times  the  objective  for  a  long  haul  system  and  corr¬ 
esponds  to  an  effective  fade  margin  of  only  27  dB. 

2.  The  primary  cause  of  oucage  is  Inband  distortion 
caused  by  the  frequency  selectivity  of  che  multipath 
fading  process. 

3.  Phase  adapclve  space  diversity  combining  is  very 
effective  in  reducing  the  amounc  of  fading.  In 
addition  it  affords  soma  increase  in  the  effective 
fade  margin  of  the  radio,  i.a.,  it  reduces  the  sever¬ 
ity  of  Inband  distortion  for  a  given  fade  depth. 

Oucage  for  che  system  with  space  diversity  combining 
wee  2.6  x  10^4  which  is  abouc  13  times  che  long  haul 
objective. 

4.  A  simple  adapclve  linear  amplitude  equalizer  in 
conjunction  with  che  phase  adaptive  space  diversity 
combining  provides  an  additional  improvement  of  appro¬ 
ximately  a  factor  of  20,  reducing  che  multipath  oucage 
to  levels  compatible  with  long  haul  availability 
objectives. 

5.  This  same  equalizer  when  applied  co  a  non  diver¬ 
sity  channel  provides  an  improvement  of  abouc  2.  Thus 
che  preconditioning  effect  of  phase  adaptive  specs 
diversity  combining  is  necessary  to  reduce  distortions 
sufficiently  chac  they  can  be  dealc  with  by  a  relati¬ 
vely  simple  adapclve  equalizer. 

Introduction 

Spurred  by  che  cosc  advantages  of  digital 
multiplex  and  switching  there  is  an  accelerating  move 
cowards  all  digital  necvortcs.  Development  of  digital 
modulation  techniques  with  bsndwldth  efficiency 
greater  chan  2  blcs/hertz  has  made  digital  radio 
transmission  a  practical  way  of  implementing  che 
trunks  in  such  a  network.  In  Canada,  che  3  GHz  band 
has  been  rechannelized  for  40  MHz  channels  co 
facilitate  digital  transmission. 


The  Northern  Telecom  DRS-3  Digital  Radio  System 
was  developed  for  long  haul  transaction  in  this  band. 
Tfie  DRS-8  System  consists  of  the  MR-3A  MULDEM  which 
synchronizes  and  conditions  two  DS3  signals  into  a 
single  91.04  Mb/s  digital  stream  (1344  FCM  voice 
channels)  ,  che  CS-4  1  x  S  Automatic  Protection 
Switching  which  provides  equipment  protection  and  the 
RD-3  Digital  Radio  which  transmits  91.04  Mb/s  in  a  40 
MHz  RF  bandwidth  using  Quadrature  Partial  Response 
Signalling  (QPRS).  To  achieve  che  maximum  economic 
advantage,  che  system  Is  designed  to  be  overbuilt  on 
existing  4  GHz  routes  using  existing  sites,  buildings, 
towers  and  antennas.  This  fixes  the  hop  lengths  ac 
approximately  45-35  km.  The  availability  objective  Is 
99.982  for  a  6560  km  connection.  Of  the  allowable 
.022  outage,  .012  was  allocated  for  equipment  failure, 
.0052  for  rain  and  .OOSZ  for  multipath  fading.  This 
requires  a  per  hop  outage  probability  of  1.43  x  L0"& 
in  the  worst  fading  month.  To  design  a  system  to 
these  stringent  requirements  more  information  was 
required  on  the  effect  of  multipath  propagation  on 
wideband  digital  transmission  over  long  hops. 

Although  information  has  been  published4  on  the 
frequency  selective  nature  of  fading,  little  or  no 
quantitative  information  was  available  on  the  effects 
of  frequency  selective  fading  on  digital  radio  system 
availability.  A  series  of  digital  propagation  trials 
have  been  conducted  over  a  51  km  path  ac  3  GHz  co 
provide  more  information  on  space  diversity,  inband 
distortion,  and  cross  polarlzaclon  discrimination 
during  mulclpath  fading. 

Space  Diversity 

Space  diversity  is  required  co  aeec  che  long 
haul  availability  objective  ac  3  GHz.  Frequency 
diversity  which  is  normally  used  to  combat  fading  in 
the  4  GHz  band  would  be  a  more  economical  solution  but 
it  does  noc  provide  adequace  improvement  for  a  fully 
bullc-ouc  system  because  of  che  increased  fading 
activity  and  smaller  fractional  channel  separation 
at  3  GHz. (I)  Space  diversity  switching,  either  ac 
RF  or  IF,  is  unsuitable  for  a  digital  radio  syscem 
because  of  the  excessive  number  of  interruptions  chac 
such  an  arrangement  would  produce  during  fading. 
Although  post  regenerative  switching  can  be  made 
hitless,  it  is  very  expensive  since  two  complete 
receiver-demodulators  are  required. 

Adapclve  inphase  IF  combining  (Fig.  0  was 
selected  because  it  is  hicless  and  requires  duplica¬ 
tion  of  only  che  mixer-preamplifier  and  channel 
dropping  filter.  The  combiner  consists  of  a  summing 
Junction,  an  electronic  phase  shifter,  and  an  electro¬ 
nic  cootrol  syscem  co  maincaln  che  two  IF  tnpuc  sig¬ 
nals  in  phase.  Because  che  wanted  signals  are  added 
coherently  and  the  noise  from  che  two  receivers  adds 
incoherently,  che  combiner  yields  a  becter  signal  to 
noise  ratio  than  switched  space  diversity  whenever 
the  difference  in  signal  levels  is  less  then  7.7  dB. 


Reprinted  from  1978  International  Conference  on  Communications 
Record,  pp.  33.5.1  -  33.5.6. 
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A  consequence  or  adjusting  :h«  phase  racher  chan  the 
da  lay  of  che  input  signals  is  :hac  cha  bandwidth  of 
:ha  systaa  is  invarsaly  proportional  to  the  absolute 
time  delay  batwean  signals.  If  tha  dynamic  delay 
difference  between  ancannas  during  multipath  fading  is 
large,  significant  signal  distortion  results,  to 
check  this  and  co  confirm  thac  the  response  of  Che 
combiner  was  fast  enough  to  follow  fading  activity,  it 
was  desirable  to  verify  cha  performance  of  the 
combiner  under  accual  multipath  fading  conditions. 

Inband  Distortion 

If  the  only  effect  of  oulcipath  fading  were  a 
reduction  in  the  receive  signal  level  during  fading, 
it  would  be  possible  to  calculate,  using  Barnect's 
empirical  fading  formulas^*),  chat  the  fade  margin 
required  co  achieve  the  oucage  objective  over  an 
average  50  km  pach  wich  a  nondiversicy  system  would  be 
56  dB.  Similarly,  applying  Vigant's  expression^) 
for  space  diversity  improvement,  the  fade  margin 
required  for  a  swicched  space  diversity  syscem  with  a 
9  n  separation  becween  ancennas  and  a  diversity 
antenna  syscem  gain  of  2  dB  less  chan  the  main,  would 
be  37  dB.  However,  ic  is  known  that  mulcipach  fading 
is  accompanied  by  severe  amplitude  and  delay  distor¬ 
tion.  Sabler^4)  presented  decailed  statistics  of 
amplitude  distortion  on  a  20  MHz  channel  at  6  GHz  for 
non  diversity  and  swicched  space  diversity  syscems 
from  which  it  is  apparent  thac  amplitude  slope  could 
be  a  serious  problem  for  boch  types  of  syscems.  A 
limited  amounc  of  information  was  available  for  delay 
distortion  on  a  nondiversity  system^,  6)4  -jha 
published  information  was  not  sufficient  co  predict 
the  availability  of  the  QRS-3  syscem  because  of  the 
unknown  affect  of  inphase  combining,  che  wider  band¬ 
width  of  the  SD-3  system  and  che  limited  information 
available  on  delay  distortion. 

Cross-Pdlarlzatlon  Discrimination 

3arnecc<?)  reported  severe  deterioration  of 
cross-polarization  discrimination  (XPD)  during  mulci¬ 
pach  fading.  This  effect  combined  wich  differential 
up  fading  of  adjacent  channels  could  be  a  problem  in  a 
system  which  relies  on  cross-polarlzacion  discrimin¬ 
ation  co  provide  some  protection  from  adjacent  channel 
interference.  Since  this  effect  will  also  be  modified 
by  the  space  diversity  combiner,  field  measurements 
were  desirable  to  confirm  chat  the  RD-3  objective  for 
tolerance  to  adjacent  channel  interference  (1  dB 
degradation  in  thermal  fade  margin  for  two  adjacent 
channels  at  0  dB  C/1)  was  adequate. 

To  determine  che  Impact  of  mulcipach  propaga¬ 
tion  on  system  availability  and  provide  daca  to  guide 
remedial  accion,  an  extensive  experimental  program  was 
undertaken  during  the  summers  of  1975  co  1977.  The 
remainder  of  this  paper  reports  the  results  of  these 
experiments. 

Ex2erimant_DMC£i£tion 

Experiments  have  been  conducted  at  3  GHz  over 
the  past  3  years  on  the  51  km  pach  from  Kempcville  to 
Avonmore,  Ontario,  which  is  part  of  che  Bell  Canada  4 
GHz  analog  radio  syscem.  The  path  profile  and  the  S 
GHz  first  fresnel  zone  for  this  hop  (Fig.  2),  show 
average  terrain  and  adequate  clearance.  An  early 
prototype  of  the  RD-3  Digital  Sadia  was  installed 
using  the  existing  antennas  and  waveguide.  Slock 
diagrams  of  che  experiment  are  shown  in  Figures  3  and 


In  the  transmit  site  (Fig.  3),  digital  transmit¬ 
ters  transmitted  pseudo  random  bit  streams  on  channel  4 


1,3173.15  MHz)  horizontally  polarized  as  well  as  on  the 
two  adjacent  channels,  3  (8132.41  mHz)  and  5  (3213.39 
MHz),  vertically  polarized.  In  addition,  FM  trans¬ 
mitters,  modulated  by  distinct  single  frequency  tones, 
transmitted  on  the  horizontal  and  vertical  polariza¬ 
tions  on  channel  2  (8091.67  MHz)  for  measurement  of 
cross-polarlzacion  interference. 

In  the  receive  sice  the  signals  were  received 
by  the  standard  KB  horn  ancenna  and  a  space  diversity 
dish  3  m  in  diameter  locaced  9  a  below  the  horn 
ancenna.  On  channel  4,  three  digital  receivers  were 
used  to  demodulate  and  regenerate  signals  from  the 
main  antenna,  space  diversity  ancenna,  and  the  oucpuc 
of  the  space  diversity  combiner.  Error  races  on  the 
three  digital  channels  were  measured  by  bit-by-bit 
comparison  with  locally  synchronized  replicas  of  the 
transmitted  pseudo  random  sequences.  The  level  on 
these  channels  was  monitored  by  the  linearized  AGC 
control  volcage  of  the  IF  amplifiers.  In  addition, 
narrowband  decectors  ac  -13  MHz  from  the  center  of  the 
channel  were  included  to  give  an  indication  of  inband 
amplitude  distortion.  The  adjacenc  channels  were 
monitored  for  level  only.  To  provide  indications  of 
equipment  malfunction,  alarms  uere  broughc  ouc  for 
carrier  recovery  lock,  clock  recovery  lock,  and  the 
primary  power  supply.  The  data  was  recorded  by  a 
PDP/11-E10  data  acquisition  system,  which  sampled  all 
of  che  measured  parameters  every  100  ms  and  wrote  a 
record  on  disk  when  errors  or  significant  changes  in  a 
parameter  occurred.  The  contents  of  the  disk  were 
periodically  transferred  by  a  high  speed  daca  line  co 
che  IBM  370  compucer  ac  che  3ell-Morchern  Research 
Laboratories  in  Ottawa,  where  the  data  was  scored  on 
tape  for  later  analysis. 

For  che  cross-polarization  experiment,  relative 
levels  of  the  tones  ac  che  oucpuc  of  che  FM  demodul¬ 
ators  provided  a  measure  of  the  carrier-co-incerf st¬ 
ance  ratio  for  both  horizontal  into  vertical  and 
vertical  inco  horizontal  interference  on  channel  2. 
Carrier-coincerf erence  ratio  is  che  combined  effect 
of  differencial  fading  between  horizontal  and  vertical 
polarizacion  and  deterioration  of  cross-polarization 
discrimination. 

In  a  supplementary  experiment  co  provide  deta¬ 
iled  in-band  distortion  characterization  during  mulci¬ 
pach  fading,  a  swept  microwave  link  analyzer  with 
sample  and  hold  circuit  at  the  oucpuc  of  che  receiver 
was  used  to  sample  the  amplitude  and  delay  distortion 
ac  5  frequencies  across  a  -15  MHz  band. 

Propagation  Experiment  Results  -  Space  Diversity 

Figure  5  shows  the  probability  distribution  of 
fading  for  che  main,  diversity  and  combined  channels. 
Comparison  of  che  fading  distribution  of  che  main 
channel  to  thac  predicted  for  an  average  hop  using 
Barnett's  formula(2)  shows  che  fading  to  be  typical. 
The  jolne  fading  distribution  for  main  and  space 
diversity  ancennas,  which  represents  the  performance 
of  an  ideal  comparative  space  diversity  switch,  shows 
that  che  combiner  has  slighcly  becter  performance  chan 
the  comparative  switch.  A  mors  exact  evaluation  of 
the  IF  combiner  performance  was  obcained  by  comparing 
its  oucpuc  level  distribution  with  thac  of  the  ideal 
in-phase  combiner  computed  from  che  main  and  space 
diversity  levels,  sample  by  sample.  The  measured  fade 
distribution  was  found  co  accurately  coincide  with 
chat  of  the  ideal  combiner.  To  confirm  thac  ic  was 
accurately  cracking  dynamic  variations  in  the  input 
signals,  internal  voltages  which  give  a  measure  of  a 
phase  error,  were  monicored  over  a  period  of  several 
weeks.  This  cesc  confirmed  thac  the  dynamic  response 
of  the  space  diversity  combiner  was  satisfactory. 


Digital  Performance 

The  performance  of  che  digital  traaaaiaalon 
svstem  can  be  described  by  two  complementary  critarla: 
cha  outaga  probability,  daflnad  as  tha  probability  of 
bit  arror  rata  graatar  chan  10“4  and  tha  affacciva 
fada  margin,  daflnad  as  that  fada  dapch  which  has  cha 
sana  probability  as  cha  obsarvad  probability  of 
oucage.  Tha  affacciva  fada  aargin  is  a  convanlant  way 
of  normalizing  cha  outaga  to  cha  aaounc  of  fading 
obsarvad,  tharaby  providing  a  aaaningful  rafaranca  and 
facilitating  extrapolation  to  ocher  situations.  The 
difference  between  cha  effective  fade  aargin  and  cha 
flac  fada  aargin  of  the  equipaanc  neasurad  with  an 
accenuator  is  a  aaasure  of  che  severity  of  salaetiva 
fading  effects.  These  quantities  are  shown  for  cha 
three  digital  channels  in  Table  I. 

TABLE  1 

RADIO  PERFORMANCE  WITHOUT  ADAPTIVE  EQUALIZATION 

Oucage  Effective 

Probability  Fade  Margin 

Main  1.2  x  10*3  26.5  dB 

Space  Diversity  7.S  x  10*4  23.5  dB 

Coablnad  2.6  x  10“5  30.0  dB 

Joint  1.9  x  10”5 

Non-Diversity  Channel 

The  outage  on  che  non-diversity  channels  exceed 
che  objective  by  a  factor  of  about  700  which  corres¬ 
ponds  co  an  effective  fada  aargin  of  only  27  dB.  The 
measured  flac  fade  aargin  of  cha  radio  was  AS  dB.  This 
dearly  lndicacss  chac  a  design  aachod  us  lag  cha  flac 
fade  aargin  in  conjunction  with  empirical  formulas  for 
fading  will  not  provide  a  good  estimate  of  cha  availa¬ 
bility  of  a  digital  radio  system  over  paths  of  this 
length.  The  narrowband  dataccors  at  +13  MHz  give  an 
indication  of  linear  and  quadratic  amplicude  distor¬ 
tion  during  fading.  In  Figure  7,  che  unconditional 
probability  distribution  of  amplitude  slope  is  com¬ 
pared  with  the  distribution  conditioned  on  oucaga  of 
che  digical  receiver  for  a  non-diversity  system.  This 
illustrates  che  high  degrea  of  correlation  of  outaga 
and  amplitude  slope.  Although  quedracic  distortion 
does  noc  show  a  similar  correlation,  cha  iniclal  uni¬ 
form  slope  of  cha  conditional  distribution  suggests 
chac  ocher  factors  such  as  group  delay  distortion  aay 
be  contributing  causes  of  oucage.  Because  of  an  equip¬ 
ment  malfunction  a  similar  comparison  is  noc  possibla 
for  cha  combined  channal. 

Initially  chare  was  concern  chac  cha  aodulacion 
schema,  QPRS,  was  more  suscaptibla  to  Inband  distor¬ 
tion  chan  ochar  methods;  however  this  proved  not  co  ba 
che  case.  Since  cha  dominant  multipath  affects  are 
odd  ordar  distortions  which  causa  crosstalk  baevaan 
in-phase  and  quadrature  channel,  all  cwo  dimensional 
modulation  systems  such  as  Quadrature  Amplitude 
Modulation  or  Multiphase  PSK  which  have  comparable 
noisa  resistance  are  similarly  affected.  This  was 
confirmed  by  computer  simulation^® )  for  8  PSK  which 
was  che  ochar  major  contending  modulation  system  for 
this  application.  Anochar  question  of  lnceresc  is  why 
che  severe  Inband  distortion  experienced  in  these 
exparlmancs  was  never  a  problem  for  analogue  radio 
transmission  systems.  Four  factors  are  significant, 
first  che  raw  fading  is  only  half  as  severe  at  A  CHz 
compared  co  8  CHz.  Second,  che  larger  frequency 
dlversley  improvement  (8  times  higher  then  ac  8  CHz) 


reduces  che  A  CHz  fading  eceivlty  even  further. 
Thirdly,  Frequency  Modulation,  and  modulation  method 
commonly  used  for  analogue  transmission,  is  inherently 
more  resistant  co  odd  order  inband  distortion  than  che 
digital  modulation  syscems  which  are  suitable  for 
bandwidth  efficient  transmission.  Last  of  all  che 
susceptibility  of  FM  systems  co  inband  distortion 
increases  with  the  loading,  that  is  che  number  of 
voice  channels  carried.  Since  most  frequency  selective 
multipath  fading  occurs  during  che  night  when  che 
network  is  lightly  loaded,  the  effect  is  reduced. 
Digical  aodulacion  systems  on  the  ocher  hand  employ 
scrambling  devices  which  ensure  chac  che  system  is 
fully  loaded  all  the  dme. 

Combined  Channel 

The  outage  on  che  combined  channal  exceeded  che  per 
hop  objective  by  a  factor  of  18.  Note  that  in  addi¬ 
tion  co  providing  an  improvement  in  che  dme  below 
level  for  the  combined  channel,  che  combiner  improves 
che  effective  fade  margin  from  27  co  30  dB.  That  is, 
lc  provides  some  measure  of  relief  from  selective 
fading  effects.  Frequency  diversity  would  provide 
some  improvement  for  an  initial  1  x  1  system  but  che 
improvement  for  a  fully  built  out  system  is  very 
small.  The  evaluate  methods  of  further  improving  che 
availability  che  detailed  distortion  characteristics 
of  the  combined  channel  was  obtained  by  means  of  a 
supplementary  experiment  in  which  che  sampled 
amplitude  and  delay  outputs  of  a  swept  link  analyzer 
ware  recorded.  Figure  3  and  9  show  che  probability 
distributions  of  che  coefficients  of  a  3rd-order 
least-squares  Chebyshev  polynominal  fit  co  cha 
amplicude  and  delay  measured  ac  5  points  ovar  a  +15 
MHz  band  on  the  combined  channel.  Note  chat  che 
linear  amplitude  component  dominates  the  amplitude 
distortion  and  chac  delay  distortions  are  very  small. 
These  results  suggest  chac  a  simple  amplitude  slope 
equalize  would  provide  a  substantial  improvement  in 
availability.  To  confirm  this,  a  linear  analysis 
program^)  was  employed  co  compute  Che  eye  opening 
with  and  without  a  linear  amplicude  equalizer  for  each 
data  record  which  showed  a  significant  distortion  on 
che  combined  channel.  The  probability  distributions. 
Figure  10,  show  that  che  improvement  in  the  0-10Z 
range  of  eye  opening  is  approximately  20,  which 
indicates  chac  a  linear  amplitude  equalizer  in 
conjunction  with  che  space  diversity  combiner  would 
provide  sufficient  improvement  to  meet  che  long  haul 
availability  objective. 

Adaptive  Equalizer  Field  Trial 

An  adaptive  linear  amplitude  equalizer  was 
developed  and  during  che  summer  of  1976,  che  perfor¬ 
mance  of  a  prococype  was  evaluated  on  che  cest  hop. 

In  this  period  there  were  several  perloda  of  anomalous 
fading.  This  fading  was  characterized  by  prolonged 
simultaneous  signal  depressions  on  both  antennas  on 
the  main  digical  channel  and  on  both  polarizations 
of  che  cross-polarization  experiment,  80  MHz  away.  In 
one  instance,  boch  main  and  space  diversity  antennas 
were  continuously  faded  below  20  dB  for  8A  minutes, 

221  times  the  mean  fade  duracion  reported  by 
VlganctO).  Such  periods  are,  of  course, 
characterized  by  poor  space  diversity  improvement  and 
excessive  fading  on  che  combined  channal.  Since  cha 
fading  behaviour  during  these  periods  is  different 
from  anything  encountered  during  the  rest  of  the  two 
years  of  experiments  and  is  different  from  chac  repor¬ 
ted  in  che  extensive  liceratura  on  aulclpach  fading, 
it  is  concluded  chst  a  mechanism  ocher  chan  aulcipach 
fading  was  responsible.  .  Rain  can  be  excluded  since 
che  fading  exhibits  frequency  selectivity.  The 
resulcs  of  an  experiment  conducted  concurrently  on  an 
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antenna  15  dfl  below  grazing  show  cnac  these  events  jtd 
aoc  caused  by  acaospheric  subref activity.  In  vlaw  of 
chaaa  facta,  tha  data  vara  analyzed  separately,  with 
and  without  thasa  anomalous  events. 

With  tha  anoaaloua  fading  included,  che  overall 
performance  of  the  RD-3  radio  bay  with  adaptive 
equalizer  la  satisfactory.  Tha  affective  fada  margin 
of  the  equalized  combined  channel  for  tha  daca  sat  was 
10  dS  which  exceeds  tha  objective  of  37  dB.  ^Iftegnice 
of  this,  the  outage  exceeds  che  objective  by\5.5y 
because  of  extreme  fading  on  che  combined  chahnel  (23 
times  higher  than  in  che  1975  experiment) .  The  cause 
of  the  excessive  fading  was  low  space  diversity  impro¬ 
vement  caused  by  high  correlation  of  che  fading  exper¬ 
ienced  on  che  main  and  space  diversicy  antennas.  These 
observations  reflect  the  presence  in  this  daca  sec  of 
a  disproportionate  amount  of  deep  flat  fading.  Both 
che  combiner  and  adaptive  equalizer  fall  to  show  the 
expected  improvement  during  such  fading  and  che  effec¬ 
tive  fade  margin  appears  better  than  ic  would  be  for 
more  typical  fading.  The  1  GHz  system  also  suffered 
abnormally  long  oucages  on  these  evenings.  Since  anom¬ 
alous  fading  does  not  creace  a  problem  on  the  1  GHz 
network  we  conclude  chat  ic  is  a  sufficiently  rare 
occurance  on  a  network  wide  basis  to  noc  have  a  major 
impact  on  system  availability. 

With  che  anomalous  fading  activity  excluded  che 
fading  on  main  and  space  diversicy  antennas  was  only 
slightly  less  chan  for  che  hypothetical  reference  hop 
down  co  approximately  40  dB.  The  combined  channel 
exhibited  fairly  cypical  fading  down  co  approximately 
30  dB  after  which  tha  probability  falls  off  rapidly. 

No  fades  greater  chan  33  dB  were  experienced.  This  is 
cypical  of  che  original  experimental  results  of 
Vlgancs.  The  effective  fade  margin  of  tha  unequalized 
channel  is  30  dB  which  is  equal  to  the  value  for  the 
1975  experiment.  In  che  selected  data  set  134  seconds 
of  fading  induced  outage  were  recorded  on  che  unequa¬ 
lized  combined  channel  compared  co  0.3  seconds  on  che 
equalized  channel.  The  calculated  improvement  of  137 
is  probably  not  an  accurate  escimate  of  tha  average 
improvement  because  of  che  small  amounc  of  equalized 
channel  outage  on  which  it  is  based.  However,  since  it 
is  an  order  of  magnitude  better  than  required  we  may 
safely  conclude  chac  che  improvement  is  satisfactory. 

A  saparace  experiment  to  evaluate  che  performance  of 
che  linear  amplitude  equalizer  on  a  non-diversity 
channel  yielded  an  effective  fade  margin  of  only  30 
dB.  This  Indicates  that  che  combiner  plays  an  essen- 
dal  pact  in  dealing  with  inband  distortion  and  chac  a 
simple  amplitude  slope  equalizer  is  noc  satisfactory 
for  a  non-diversity  system. 

Cross-Polarization  Discrimination 

Results  from  che  cross-polarlzacion  discrimin¬ 
ation  (XPD)  portion  of  che  experiment  confirmed  che 
severe  deterioration  of  XPD  on  a  non-diversity  system 
reporced  by  Barnecc(Z).  However  the  XPD  performance 
on  a  combined  channel  obtained  numerically  by  procee- 
slng  the  simultaneous  measurements  on  che  two  antennas 
showed  chat  che  XP0  was  considerably  improved  by 
combining.  The  performance  was  good  enough  co  suggest 
chac,  with  soma  further  work,  co-frequency  croes-pola- 
rlzad  operation  would  be  possible  permitting  increa¬ 
sing  the  capacity  of  che  3  GHz  bend  from  5  to  12  chan¬ 
nels. 

Conclusion^ 

An  experimental  investigation  into  che  effects 
of  multipath  propagation  on  a  wideband  digical  radio 
system  over  a  51  km  hop  hae  been  carried  out  at  3  GHz. 

Ic  wee  found  chat  che  affects  of  frequency  selective 
fading  result  in  unacceptably  high  system  unavaila¬ 
bility  unless  adaptive  equalization  and  space 
diverelty  are  employed.  g.5 
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operates  over  long  hops  using  a  two  dimensional 
modulation  system  such  as  QAM  or  multiphase  PSk. 

Specifically  ic  is  conclud'd  chac: 

1.  Multipath  induced  outage  is  much  higher  chan 
predicted  from  che  measured  flat  fade  margin  of  che 
equipment.  For  a  non-diversity  system,  che  probabi¬ 
lity  of  oucage  for  a  single  hop  in  che  worst  fading 
month  is  approximately  1  x  10”J.  This  is  700  times 
che  objective  for  a  long  haul  syscem  and  corresponds 
co  an  effective  fade  margin  of  only  27  dB. 

2.  The  primary  cause  of  outage  is  inband  distortion 
caused  by  the  frequency  selectivity  of  che  mulclpach 
fading  process. 

3.  Phase  adaptive  space  diversicy  combining  is  very 
effective  in  reducing  che  amounc  of  fading.  In  add¬ 
ition  it  affords  some  Increase  in  che  effective  fade 
margin  of  che  radio  i.e.,  it  reduces  che  severity 
inbend  distortion  for  a  given  fade  depth.  Probability 
of  outage  for  the  system  with  space  diversity  com¬ 
bining  was  2.6  x  10”5  which  is  about  18  times  the 
long  haul  objective. 

4.  A  simple  adaptive  linear  amplitude  equalizer  in 
conjunction  with  phase  adaptive  space  diversicy  com¬ 
bining  provides  an  additional  improvement  of  approxi¬ 
mately  20,  reducing  che  mulclpach  oucage  co  levels 
compatible  with  long  haul  avallablllcy  objectives. 

5.  This  same  equalizer  when  applied  co  a  non-diver¬ 
sity  channel  provides  an  improvement  of  about  2. 

the  preconditioning  effect  of  the  phase  adaptive  space 
diversity  combiner  is  necessary  co  reduce  distortion 
sufficiently  chat  they  can  be  dealt  with  by  a  relati¬ 
vely  simple  adapeive  equalizer. 
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AN  ADAPTIVE  EQUALIZER  FOR  CORRECTION  OF  MULTIPATH  DISTORTION  IN  A  90  MB/S  8  PSK  SYSTEM 
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Collins  Transmission  Systems  Oivislon 
Rockwell  International 


ABSTRACT 

Oigital  transmission  over  LOS  microwave  is  now 
a  well  established  technique.  As  the  number 
of  digital  microwave  systems  has  grown,  there 
has  been  an  increasing  need  for  an  adaptive 
equalizer  to  correct  for  the  amplitude  distor¬ 
tion  that  arises  during  multipath  propagation 
conditions. 

This  paper  describes  a  series  of  lab  experiments 
designed  to  verify  the  nature  and  extent  of 
degradation  to  the  8  PSK  signal  arising  from 
a  two-ray  multipath  geometry.  It  also  describes 
the  performance  improvement  that  can  be  obtained 
with  an  adaptive  equalizer. 


INTRODUCTION 

The  use  of  digital  microwave  for  transmission 
of  PCM  voice  signals  is  expanding  at  a  rapid 
rate.  As  the  data  rate  and  bandwidth  efficiency 
of  digital  microwave  systems  has  increased,  the 
distortion  effects  of  multipath  have  become  an 
increasingly  serious  problem.  Adaptive  equalizer 
circuits  have  been  developed  that  correct  the 
more  serious  distortion  effects. 

This  paper  describes  an  adaptive  equalizer, 
a  multipath  model,  effects  of  multipath  on  a 
digital  radio  system,  and  the  system  performance 
improvements  that  can  be  effected  through  the 
use  of  an  adaptive  equalizer. 

MULTIPATH  MODEL 

Multipath  propagation  studies  in  the  past  have, 
in  general,  concentrated  on  the  received  carrier 
level  for  analog  transmission.  In  analog 
transmission  systems,  the  distortion  effects  of 
multipath  were  not  so  important  because  the 
degradation  from  the  distortion  effects  were 
graceful  and  the  loading  conditions  on  the 
analog  system  were  often  light  during  multipath 
propagation. 

In  digital  systems,  however,  the  distortion 
mechanisms  are  more  critical  because  the  system 
has  data  scrambling,  which  means  that  the  loading 
is  constant  at  all  times.  With  scrambling,  the 
transmit  spectrum  is  well  defined  and,  therefore, 
multipath  distortion  condltions'can  be  much 


Figure  1.  8asic  Refraction  Geometry  that 

Results  in  Multipath  Propagation. 


more  readily  measured  and  corrected  than  in 
analog  systems. 

Over  the  years  many  multipath  propagation  models 
have  been  proposed.  These  models  vary  from 
simple  two-  and  three-ray  models  to  elaborate 
multi -ray  models.  The  two-ray  model  is  the  most 
understood  model  for  multipath  propagation  and 
it  accounts  for  most  of  the  distortion  effects 
that  are  seen  in  multipath  propagation.  There¬ 
fore,  the  two-ray  model  was  used  for  this  inves¬ 
tigation.  This  model  is  shown  in  Figure  1.  It 
consists  of  a  direct  ray,  and  a  single  reflected 
or  refracted  secondary  ray.  These  two  rays  add 
at  some  frequencies  and  cancel  at  others.  The 
amplitude  and  group  delay  effects  that  are  caused 
by  the  two-ray  multipath  are  shown  in  Figure  2. 


f 


Figure  2.  Amolitude  and  Delay  Characteristics 
in  Multipath. 
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The  spacing  between  amplitude  nulls  is  inversely 
proportional  to  the  absolute  delay  difference  (t) 
between  the  direct  and  reflected  ray.  The  null 
depth  is  determined  by  the  relative  amplitudes 
of  the  direct  and  reflected  rays. 

EFFECTS  Of  MULTIPATH  ON  DIGITAL  SYSTEMS 


In  digital  radio  systems,  the  bit  error  rate  (BER) 
during  multipath  is  controlled  by  intersymbol 
interference  resulting  from  frequency  dependent 
amplitude  and  group  delay  effects.  To  evaluate 
these  effects  on’digital  transmission,  a  series 
of  experiments  was  conducted  which  measured  the 
threshold  degradation  resulting  from  linear 
amplitude  slopes,  linear  group  delay  slopes  and 
Symnetrical  parabolic  group  delay.  These  measure¬ 
ments  were  conducted  on  a  6  GHz  90  MSS  radio 
system  operating  at  3  bits/Hz.  Figure  3  shows 
the  threshold  degradation  resulting  from  linear 
amplitude  effects  (gain  slope).  Figure  4  shows 
the  effects  of  linear  group  delay  slope  and 
Figure  5  shows  the  effects  of  parabolic  group 
delay. 

To  evaluate  the  importance  of  each  of  these 
effects,  it  is  necessary  to  understand  the  degree 
to  which  each  of  these  effects  is  present  in 
multipath.  Figures  6  and  7  show  the  calculated 
amplitude  and  group  delay  slopes  in  a  30  MHz 
bandwidth  resulting  from  two-ray  10  d8  fades  at 
various  absolute  delays.  It  can  be  seen  from 
Figures  3-7  that  for  a  given  delay  and  fade  depth, 
the  effect  of  amplitude  slope  on  threshold  degra¬ 
dation  is  more  drastic  than  that  due  to  delay 
slope  or  parabolic  delay.  This  implies  that 
amplitude  slope  is  the  primary  source  of  perfor- 


Figure  3.  3ER  Threshold  Degradation  for  Linear 
Amplitude  Slope  on  “DR-6. 


Figure  4.  BER  'hreshold  Degradation  for  Linear 
Group  Delay  on  TOR-6. 
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Figure  5.  BER  Threshold  Degradation  for 

Symmetrica!  Parabolic  Group  Delay 
on  MDR-6. 
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Figure  6.  Amplitude  Slope  vs.  Phase  Angle  of 
Reflected  Signal  for  10-dB  Fade. 


The  amplitude  slope  correction  circuit  can  correct 
for  positive  or  negative  values  of  amplitude 
slope  of  up  to  12  dB  over  a  20  ^Hz  bandwidth.  The 
null  or  bow  correction  circuit  can  correct  for  up 
to  12  dB  of  amplitude  null. 

The  amplitude  slope  correction  circuit  has  essen¬ 
tially  no  group  delay  effects.  The  null  correc¬ 
tion  circuit  provides  the  reouired  group  delay 
characteristics  necessary  to  compensate  for  the 
group  delay  associated  with  the  amplitude  null. 

SYSTEM  PERFORMANCE 


To  evaluate  the  performance  of  an  adaptive 
equalizer  under  multipath  conditions,  a  series 
of  experiments  were  performed  that  simulated 
the  passband  distortion  arising  from  two-ray 
multipath.  The  results  of  these  experiments 
verified  that  amplitude  slope  is  the  most  serious 
distortion  that  arises  in  multipath.  Thus,  the 
most  important  part  of  the  adaptive  equalizer  is 
an  amplitude  slope  correction  circuit.  The 
effectiveness  of  this  circuit  is  shown  in  Figures 
8  and  9.  Figure  3  shows  the  8ER  threshold 
improvement  using  an  adaptive  equalizer  in  the 
presence  of  amplitude  slope.  The  upper  curve  is 
the  result  of  a  severe  multipath  condition  which 
resulted  in  errored  performance  at  all  signal 
levels.  When  the  adaptive  equalizer  is  turned  on, 
the  threshold  performance  returns  to  within  0.5  dB 
of  its  normal  value. 
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Figure  7.  Oelay  Slope  vs.  Phase  Angle  of  Reflected 
Signal  for  10-d8  Fade. 


mance  degradation  during  multipath.  Also,  it  can 
be  shown  that  when  the  amplitude  slope  is  removed, 
the  next  most  serious  problem  is  a  symmetrical 
null  at  the  center  of  the  passband. 

AOAPTIVE  EQUALIZER 

The  adaptive  equalizer  circuit  is  incorporated 
in  the  receiver  70  Mhz  IF  amplifier.  It  consists 
of  an  amplitude  slope  detector  which  controls  an 
amplitude  slope  correction  circuit  and  a  symmetri¬ 
cal  null  detector  which  controls  the  null  correc¬ 
tion  circuit.  The  slope  detector  operates  by 
comparing  the  amplitude  near  the  high  end  of  the 
passband  with  the  amplitude  near  the  low  end  of 
the  passband.  The  linear  amplitude  slope  correc¬ 
tion  circuit  maintains  the  desired  passband  slope. 
The  null  detection  circuit  operates  by  comparing 
the  energy  at  the  center  of  the  passband  with  the 
total  energy  in  the  passband. 


Figure  8.  BER  Threshold  Improvements  Using  An 
Adaptive  Eaualizer  in  the  Presence  of 
Amplitude  Slope. 


Figure  9  shows  the  BER  threshold  improvement  with 
an  adaptive  equalizer  in  the  presence  of  a 
symmetrical  null.  Errored  performance  at  all 
signal  levels  is  again  present  in  the  upper  curve 
and  normal  threshold  performance  is  obtained  when 
the  adaptive  equalizer  is  turned  on. 

The  cumulative  effects  of  amplitude  slope  and 
nulls  are  shown  in  Figure  10.  The  two  outer 
lobes  show  the  performance  without  adaptive 
equalizers.  The  peak  of  these  two  lobes  corres¬ 
pond  to  the  arrival  angles  that  result  in  maximum 
amplitude  slope.  When  the  slope  correction  cir¬ 
cuit  is  applied,  error  free  performance  is 
obtained  at  all  arrival  angles.  As  the  fade  depth 
is  increased,  the  next  region  of  degradation  that 
is  encountered  is  the  point  where  the  amplitude 
null  is  exactly  at  the  center  of  the  passband. 
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Figure  10.  Polar  Plot  of  BER  Improvements  Using  An 
Adaptive  Equalizer. 
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The  importance  of  this  data  is  twofold.  First, 
the  usable  fade  margin  in  multipath  has  been 
increased  significantly.  Second,  the  relative 
phase  angles  of  the  multipath  at  which  transmis¬ 
sion  errors  result  are  greatly  reduced.  8oth 
of  these  factors  help  to  reduce  the  outage  time 
on  digital  systems.  It  is  anticipated  that 
adaptive  equalizers  used  in  conjunction  with 
space  diversity  will  allow  digital  systems  to 
meet  long-haul  multipath  outage  objectives.  This 
particular  adaptive  equalizer  is  now  being 
evaluated  under  actual  field  conditions. 


Figure  9.  BER  Threshold  Improvements  Using  An 
Adaptive  Equalizer  in  the  Presence 
of  a  Symmetrical  Null. 


This  point  is  shown  in  Figure  10  as  the  center 
lobe  centered  on  180s.  Note  that  the  fade  depth 
needed  to  produce  this  effect  is  5  d8  deeper  than 
that  originally  needed  to  produce  the  amplitude 
slope  problem.  When  the  null  correction  circuit 
is  turned  on,  error-free  performance  is  again 
restored  at  all  phase  angles.  Now  the  depth  of 
fade  is  again  increased  until  errors  are  present. 
This  region  is  shown  by  the  two  narrow  spikes 
in  Figure  10  and  correspond  to  an  in-band  ampli¬ 
tude  null  halfway  between  the  center  and  edge  of 
the  passband. 
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APPENDIX  A 


LIST  OF  SEMINAR  ATTENDEES 


NAME 

ORGANIZATION 

PHONE 

Paul  Hartmann 

Rockwell  International 

(214)  996-6359 

Tom  Giuffrida 

Bell  Labs 

(201)  229-6850  x2485 

Stephen  Matsuura 

USACEEIA 

VON  879-2187 

Andrew  E.  Hooper 

NUWES 

(805)  982-8971 

VON  744-4963 

Henry  C.  Merhoff 

PACMISTESTCEN 

VON  351-7881 

Adrian  Eley 

Ft.  Meade 

(301)  796-6835 

Steve  Barber 

Bell -Northern  Research 

(613)  596-5036 

Arvids  Vigants 

Bell  Labs 

(201)  229-6850  x2450 

Larry  Hause 

OOC/NT1A/TT3 

(303)  449-1000  x3945 

Robert  W.  Hubbard 

DOC/NTIA/ITS 

(303)  449-1000  x3414 

James  L.  Weblemoe 

PACMISTESTCEN 

(805)  982-8971 

VON  351-8971 

John  Osterholz 

DCEC/R410 

(703)  437-2453 

Matthew  Blanding 

DCEC/R210 

(703)  437-2266 

Jyoti  S.  Sharma 

Western  Union 

(201)  825-5148 

John  M.  Reardon 

Western  Union 

(703)  790-2248 

Ralph  Morris 

Western  Union 

(703)  790-2350 

Ed  Lyons 

Ft.  Meade 

(301)  796-6835 

John  J.  Knab 

Digital  Comm.  Corp. 

(301)  840-3468 

Pete  Plotkin 

DCEC/R310 

(703)  437-2261 

Dennis  Buck 

1842EEG  Scott  AFB ,  IL 

VON  638-4185 

C.  H.  Blackerby 

DCEC/R71 0 

(703)  437-2237 

J.  A.  Rose 

DCEC/R710 

(703)  2236 

A-l 


NAME 

Walter  Cybrowski 
Victor  Weill 
Harris  A.  Stover 
John  J.  Cortnack 
M.  W.  Horowitz 
Jim  Vest 
Joe  Mensch 
H.  L.  McKinley 
D.  0.  Schultz 
Allen  E.  Post 
John  K.  Webb 
Owen  Cote' 

Tom  Shimabukuro 
Robert  Wallace 
J.  A.  0'8rien 
Gayton  Yancey 
George  H.  Hagn 
Oavid  R.  Smith 


ORGANIZATION 
DCEC/R220 
DCEC/R230 
OCEC/R220 
0CEC/R220 
0CEC/R220 
DCEC/R540 
DCEC/R220 
DCEC/R2Q0 
OCEC/R201 
MITRE 
MITRE 
MITRE 
0CEC/R71 0 
Raytheon 
DCEC/R210 
Raytheon 

SRI  International 
0CEC/R220 


PHONE 

(703)  437-2316 
(703)  437-2416 
(703)  437-2316 
(703)  437-2316 
(703)  437-2316 
(703)  437-2441 
(703)  437-2316 
(703)  437-2466 
(703)  437-2466 
(617)  271-4885 
(617)  271-2596 
VON  478-3237 
(703)  437-2348 
(617)  762-6700 
(703)  437-2266 
(617)  762-6700 
(703)  524-2053 
(703)  437-2316 


